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. The Casimir force for ciiarge-neutral, perfect conductors of non-planar geometric configurations have been 

' investigated. The configurations are: (1) the plate-hemisphere, (2) the hemisphere-hemisphere and (3) the 

I spherical shell. The resulting Casimir forces for these physical arrangements have been found to be attractive. 

fSJ . The repulsive Casimir force found by Boyer for a spherical shell is a special case requiring stringent material 

^jq' property of the sphere, as well as the specific boundary conditions for the wave modes inside and outside of the 

r- { I sphere. The necessary criteria in detecting Boyer's repulsive Casimir force for a sphere are discussed at the end 

, of this investigation. 
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I. INTRODUCTION 

The introduction is divided into three parts: (1) physics, (2) applications, and (3) developments. A brief outline of the physics 
behind the Casimir effect is discussed in item (1). In item (2), major impact of Casimir effect on technology and science is 
outUned. In item (3), the introduction of this investigation is concluded with a brief review of the past developments, followed 
by a brief outline of the organization of this investigation and its contributions to the physics. 

A. Physics 

When two electrically neutral, conducting plates are placed parallel to each other, our understanding from classical electro- 
dynamics tells us that nothing should happen to these plates. The plates are assumed to be that made of perfect conductors for 
simplicity. In 1948, H. B. G. Casimir and D. Polder faced a similar problem in studying forces between polarizable neutral 
molecules in colloidal solutions. Colloidal solutions are viscous materials, such as paint, that contain micron-sized particles in a 
liquid matrix. It had been thought that forces between such polarizable, neutral molecules were governed by the van der Waals 
interaction. The van der Waals interaction is also referred to as the "Lennard- Jones interaction." It is a long range electrostatic 
interaction that acts to attract two nearby polarizable molecules. Casimir and Polder found to their surprise that there existed 
an attractive force which could not be ascribed to the van der Waals theory. Their experimental result could not be correctly 
explained unless the retardation effect was included in van der Waals' theory. This retarded van der Waals interaction or Lienard- 
Wiechert dipole-dipole interaction [1] is now known as the Casimir-Polder interaction [2]. Casimir, following this first work, 
elaborated on the Casimir-Polder interaction in predicting the existence of an attractive force between two electrically neutral, 
parallel plates of perfect conductors separated by a small gap [3]. This alternative derivation of the Casimir force is in terms of 
the difference between the zero-point energy in vacuum and the zero-point energy in the presence of boundaries. This force has 
been confirmed by experiments and the phenomenon is what is now known as the "Casimir Effect." The force responsible for 
the attraction of two uncharged conducting plates is accordingly termed the "Casimir Force." It was shown later that the Casimir 
force could be both attractive or repulsive depending on the geometry and the material property of the conductors [4-6]. 

The Casimir effect is regarded as macroscopic manifestation of the retarded van der Waals interaction between uncharged 
polarizable molecules (or atoms). Microscopically, the Casimir effect is due to interactions between induced multipole moments, 
where the dipole term is the most dominant contributor if it is non- vanishing. Therefore, the dipole interaction is exclusively 

e+ e+ e+ 

AAP - AA/^ + AA/f fUlA + AATf fVV* VaA + 

physical bare V ^ V ^ 

e~ e~ e~ 



Figure 1: The vacuum polarization of a photon. 
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referred to, unless otherwise explicitly stated, throughout this investigation. The induced dipole moments can be qualitatively 
explained by the concept of "vacuum polarization" in quantum electrodynamics (QED). The idea is that a photon, whether real 
or virtual, has a charged particle content. Namely, the internal loop, illustrated in Figure 1 , can be e+e^ , , t~^t~ , tt+tt" 

or qq pairs, etc. Its correctness have been born out from the precision measurements of Lamb shift [7, 8] and photoproduction 
[9-12J experiments over a vast range of energies. For the almost zero energy photons considered in the Casimir effect, these 
pairs last for a time interval At consistent with that given by the Heisenberg uncertainty principle AE ■ At = h, where AE is 
the energy imbalance and h is the Planck constant. These virtual charged particles can induce the requisite polarizabiUty on the 
boundary of the dielectric (or conducting) plates which explained the Casimir effect. However, the dipole strength is left as a 
free parameter in the calculations because it cannot be readily calculated. Its value can be determined from experiments. 

Once this idea is taken for granted, one can then move forward to calculate the effective, temperature averaged, energy due to 
the dipole-dipole interactions with the time retardation effect folded in. The energy between the dielectric (or conducting) media 
is obtained from the allowed modes of electromagnetic waves determined by the Maxwell equations together with the boundary 
conditions. The Casimir force is then obtained by taking the negative gradient of the energy in space. This approach, as opposed 
to full atomistic treatment of the dielectrics (or conductors), is justified as long as the most significant photon wavelengths 
determining the interaction are large when compared with the spacing of the lattice points in the media. The effect of all the 
multiple dipole scattering by atoms in the dielectric (or conducting) media simply enforces the macroscopic reflection laws 
of electromagnetic waves. For instance, in the case of the two parallel plates, the most significant wavelengths are those of 
the order of the plate gap distance. When this wavelength is large compared with the interatomic distances, the macroscopic 
electromagnetic theory can be used with impunity. The geometric configuration can introduce significant complications, which 
is the subject matter this study is going to address. 

In order to handle the dipole-dipole interaction Hamiltonian in this case, the classical electromagnetic fields have to be quan- 
tized into the photon representation first. The photon with non-zero occupation number have energies in units of tuv; where h is 
the Planck constant divided by 27r, and oj, the angular frequency. The lowest energy state of the electromagnetic fields has energy 
Huj/2. They are called the vacuum or the zero point energy state, and they play a major role in the Casimir effect. Throughout 
this investigation, the terminology "photon" is used to represent the entity with energy Hui, or the entity with energy fiu)/2 unless 
explicitly stated otherwise. 

B. Applications 

In order to appreciate the importance of the Casimir effect from industry's point of view, we first examine the theo- 
retical value for the attractive force between two uncharged conducting parallel plates separated by a gap of distance d : 
Fc = —240~^7r^ d~^/ic, where c is the speed of light in vacuum and d is the plate gap distance. To get a sense of the magnitude 
of this force, two mirrors of an area of ^ 1 crri^ separated by a distance of ^ 1 /im would experience an attractive Casimir force 
of roughly ~ 10~^ A'^, which is about the weight of a water droplet of half a milhmeter in diameter. Naturally, the scale of size 
plays a crucial role in the Casimir effect. At a gap separation in the ranges of ~ 10 nm, which is roughly about a hundred times 
the typical size of an atom, the equivalent Casimir force would be in the range of 1 atmospheric pressure. The Casimir force 
have been verified by Steven Lamoreaux [13] in 1996 to within an experimental uncertainty of 5%. An independent verification 
of this force have been done by U. Mohideen and Anushree Roy [14] in 1998 to within an experimental uncertainty of 1%. It is 
however emphasized that these experiments were not done exactly in the same context of Casimir's original configuration due to 
technical difficulties associated with Casimir's idealized perfectly flat surfaces. In 2002, Casimir's parallel plate configuration 
have been examined by G. Bressi, G. Carugno, R. Onofrio, and G. Ruoso [15]. Their force coefficient was measured at the 15% 
precision level. 

The importance of Casimir effect is most significant for the miniaturization of modem electronics. The technology already in 
use that is affected by the Casimir effect is that of the microelectromechanical systems (MEMS). These are devices fabricated 
on the scale of microns and sub-micron sizes. The order of the magnitude of Casimir force at such a small length scale can 
be enormous. It can cause mechanical malfunctions if the Casimir force is not properly taken into account in the design, e.g., 
mechanical parts of a structure could stick together, etc [16]. The Casimir force may someday be put to good use in other fields 
where nonlinearity is important. Such potential applications requiring nonlinear phenomena have been demonstrated [17]. The 
technology of MEMS hold many promising applications in science and engineering. With the MEMS soon to be replaced by the 
next generation of its kind, the nanoelectromechanical systems or NEMS, understanding the phenomenon of the Casimir effect 
become even more crucial. 

Aside from the technology and engineering applications, the Casinnir effect plays a crucial role in accurate force measurements 
at nanometer and micrometer scales [18]. As an example, if one wants to measure the gravitational force at a distance of atomic 
scale, not only the subtraction of the dominant Coulomb force has to be done, but also the Casimir force, assuming that there is 
no effect due to strong and weak interactions. 

Most recently, a new Casimir- like quantum phenomenon have been predicted by Feigel [19]. The contribution of vacuum 
fluctuations to the motion of dielectric liquids in crossed electric and magnetic fields could generate velocities of ~ 50 nm/ s. 
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Unlike the ordinary Casimir effect where its contribution is solely due to low frequency vacuum modes, the new Casimir-like 
phenomenon predicted recently by Feigel is due to the contribution of high frequency vacuum modes. If this phenomenon is 
verified, it could be used in the future as an investigating tool for vacuum fluctuations. Other possible applications of this new 
effect he in fields of microfluidics or precise positioning of micro-objects such as cold atoms or molecules. 

Everything that was said above dealt with only one aspect of the Casimir effect, the attractive Casimir force. In spite of many 
technical challenges in precision Casimir force measurements [13-15], the attractive Casimir force is fairly well established. This 
aspect of the theory is not however what drives most of the researches in the field. The Casimir effect also predicts a repulsive 
force and many researchers in the field today are focusing on this phenomenon yet to be confirmed experimentally. Theoretical 
calculations suggest that for certain geometric configurations, two neutral conductors would exhibit repulsive behavior rather 
than being attractive. The classic result that started all this is due to Boyer's work on the Casimir force calculation for an 
uncharged spherical conducting shell [4]. For a spherical conductor, the net electromagnetic radiation pressure, which constitute 
the Casimir force, has a positive sign, thus being repulsive. This conclusion seems to violate fundamental principle of physics 
for the fields outside of the sphere take on continuum in allowed modes, where as the fields inside the sphere can only assume 
discrete wave modes. However, no one has been able to experimentally confirm this repulsive Casimir force. 

The phenomenon of Casimir effect is too broad, both in theory and in engineering applications, to be completely summarized 
here. I hope this informal brief survey of the phenomenon could motivate people interested in this remarkable area of quantum 
physics. 



C. Developments 



Casimir's result of attractive force between two uncharged, parallel conducting plates is thought to be a remarkable application 
of QED. This attractive force have been confirmed experimentally to a great precision as mentioned earlier [13-15]. 

Casimir's attractive force result between two parallel plates has been unanimously thought to be obvious. Its origin can also 
be attributed to the differences in vacuum-field energies between those inside and outside of the resonator. However, in 1968, 
T. H. Boyer, then at Harvard working on his thesis on Casimir effect for an uncharged spherical shell, had come to a conclusion 
that the Casimir force was repulsive for his configuration, which was contrary to popular belief. His result is the well known 
repulsive Casimir force prediction for an uncharged spherical shell of a perfect conductor [4] . 

The surprising result of Boyer's work has motivated many physicists, both in theory and experiment, to search for its evidence. 
On the theoretical side, people have tried different configurations, such as cylinders, cube, etc., and found many more config- 
urations that can give a repulsive Casimir force [5, 20, 21]. Completely different methodologies were developed in striving to 
correctly explain the Casimir effect. For example, the "Source Theory" was employed by Schwinger for the explanation of the 
Casimir effect [21-24]. In spite of the success in finding many boundary geometries that gave rise to the repulsive Casimir force, 
the experimental evidence of a repulsive Casimir effect is yet to be found. The lack of experimental evidence of a repulsive 
Casimir force has triggered further examination of Boyer's work. 

The physics and the techniques employed in the Casimir force calculations are well established. The Casimir force calculations 
involve summing up of the allowed modes of waves in the given resonator. This turned out to be one of the difficulties in Casimir 
force calculations. For the Casimir's original parallel plate configuration, the calculation was particularly simple due to the fact 
that zeroes of the sinusoidal modes are provided by a simple functional relationship, kd = mr, where k is the wave number, d 
is the plate gap distance and n is a positive integer. This technique can be easily extended to other boundary geometries such 
as sphere, cylinder, cone or a cube, etc. For a sphere, the functional relation that determines the allowed wave modes in the 
resonator is kro = as,i, where Tq is the radius of the sphere; and as,i, the Zth root of the spherical Bessel function jg. The same 
convention is appUed to all other Bessel function solutions. The allowed wave modes of a cylindrical resonator is determined by 
a simple functional relation fcoo = /3s,;, where Oo is the cylinder radius; and (}s,i, the zeroes of cylindrical Bessel functions Jg- 

One of the major difficulties in the Casimir force calculation for nontrivial boundaries such as those considered in this in- 
vestigation is in defining the functional relation that determines the allowed modes in a given resonator. For example, for 
the hemisphere-hemisphere boundary configuration, the radiation originating from one hemisphere would enter the other and 
run through a complex series of reflections before escaping the hemispherical cavity. The allowed vacuum-field modes in the 



resonator is then governed by a functional relation k 



R'2 - R' 



nn, where 



R'2 - R' 



is the distance between two suc- 



cessive reflection points R'l and R'2 of the resonator, as is illustrated in Figure 5 of section (111) on reflection dynamics. As 
will be shown in the subsequent sections, the actual functional form for R'2 — R'l is not simple even though the physics 



behind 



k 



R'' 



rt'2 - 

R'l 



is particularly simple: the application of the law of reflections. The task of obtaining the functional relation 



nn for the hemisphere-hemisphere, the plate-hemisphere, and the sphere configuration formed by bringing in 

two hemispheres together constitutes the major part of this investigation. 

This investigation is not about questioning the theoretical origin of the Casimir effect. Instead, its emphasis is on applying 
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the Casimir effect as already known to determine the sign of Casimir force for the realistic experiments. In spite of a number of 
successes in the theoretical study of repulsive Casimir force, most of the configurations are unrealistic. In order to experimentally 
verify Boyer's repulsive force for a charge-neutral spherical shell made of perfect conductor, one should consider the case where 
the sphere is formed by bringing in two hemispheres together. When the two hemispheres are closed, it mimics that of Boyer's 
sphere. It is, however, shown later in this investigation that a configuration change from hemisphere-hemisphere to a sphere 
induces non- spherically symmetric energy flow that is not present in Boyer's sphere. Because Boyer's sphere gives a repulsive 
Casimir force, once those two closed hemispheres are released, they must repulse if Boyer's prediction were correct. Although 
the two hemisphere configuration have been studied for decades, no one has yet carried out its analytical calculation successfully. 
The analytical solutions on two hemispheres, existing so far, was done by considering the two hemispheres that were separated 
by an infinitesimal distance. In this investigation, the consideration of two hemispheres is not hmited to such infinitesimal 
separations. 

The three physical arrangements being studied in this investigation are: (1) the plate-hemisphere, (2) the hemisphere- 
hemisphere and (3) the sphere formed by brining in two hemispheres together Although there are many other boundary 
configurations that give repulsive Casimir force, the configurations under consideration were chosen mainly because of the 
following reasons: (1) to be able to check experimentally the Boyer's repulsive Casimir force result for a spherical shell, (2) the 
experimental work involving configurations similar to that of the plate-hemisphere configuration is underway [17]; and (3) to 
the best of our knowledge, no detailed analytical study on these three configurations exists to date. 

The motivation behind to mathematically model the plate-hemisphere system came from the experiment done by a group at 
the Bell Laboratory [17] in which they bring in an atomic-force-probe to a flopping plate to observe the Casimir force which 
can affect the motion of the plate. In our derivations of the equations of motion, the configuration is that of the "plate displaced 
on upper side of a bowl (hemisphere)." The Bell Laboratory apparatus can be easily mimicked by simply displacing the plate to 
the under side of the bowl, which we have not done. The motivation behind the hemisphere-hemisphere system actually arose 
from an article by Kenneth and Nussinov [25]. In their paper, they speculate on how the edges of the hemispheres may produce 
effects such that two arbitrarily close hemispheres cannot mimic Boyer's sphere. This led to their heuristic conclusion which 
stated that Boyer's sphere can never be the same as the two arbitrarily close hemispheres. 

To the best of our knowledge, two of the geometrical configurations investigated in this work have not yet been investigated by 
others. They are the plate-hemisphere and the hemisphere-hemisphere configurations. This does not mean that these boundary 
configurations were not known to the researchers in the field, e.g., [25] . For the case of the hemisphere-hemisphere configuration, 
people realized that it could be the best way to test for the existence of a repulsive Casimir force for a sphere as predicted by 
Boyer. The sphere configuration investigated in this work, which is formed by bringing two hemispheres together, contains non- 
spherically symmetric energy flows that are not present in Boyer's sphere. In that regards, the treatment of the sphere geometry 
here is different from that of Boyer. 

The basic layout of this investigation is as follows: (1) Introduction, (2) Theory, (3) Calculations, and (4) Results. The formal 
introduction of the theory is addressed in chapters (1) and (2). The original developments resulting from this investigation are 
contained in chapters (3) and (4). The brief outline of each chapter is the following: In chapter (1), a brief introduction to 
the physics is addressed; and the application importance and major developments in this field are discussed. In chapter (2), 
the formal aspect of the theory is addressed, which includes the detailed outline of the Casimir-Polder interaction and brief 
descriptions of various techniques that are currently used in Casimir force calculations. In chapter (3), the actual Casimir force 
calculations pertaining to the boundary geometries considered in this investigation are derived. The important functional relation 



for 



R'2 - R' 



is developed here. The dynamical aspect of the Casimir effect is also introduced here. Due to the technical 

nature of the derivations, many of the results presented are referred to the detailed derivations contained in the appendices. In 
chapter (4), the results are summarized. Lastly, the appendices have been added in order to accommodate the tedious and lengthy 
derivations to keep the text from losing focus due to mathematical details. 

The goal of this investigation is not to embark so much on the theory side of the Casimir effect. Instead, its emphasis is on 
bringing forth the suggestions that might be useful in detecting the repulsive Casimir effect originally initiated by Boyer on an 
uncharged spherical shell. In concluding this brief outline of the motivation behind this investigation, it must be added that if by 
any chance someone already did these work that have been claimed to be the original developments in this investigation, I was 
not aware of their work at the time of this work was being prepared. And, should that turn out to be the case, I would like to 
express my sincere apology for not referencing their work in this investigation. 



n. CASIMIR EFFECT 



The Casimir effect is divided into two major categories: (1) the electromagnetic Casimir effect and (2) the fermionic Casimir 
effect. As the titles suggest, the electromagnetic Casimir effect is due to the fluctuations in a massless Maxwell bosonic fields, 
whereas the fermionic Casimir effect is due to the fluctuations in a massless Dirac fermionic fields. The primary distinction 
between the two types of Casimir effect is in the boundary conditions among other differences. The boundary conditions 
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appropriate to the Dirac equations are the so called "bag-model" boundary conditions, whereas the electromagnetic Casimir 
effect follows the boundary conditions of the Maxwell equations. The details of the fermionic force can be found in references 
[21,24]. 

In this investigation, only the electromagnetic Casimir effect is considered. As it is inherently an electromagnetic phenomenon, 
we begin with a brief introduction to the Maxwell equations, followed by the quantization of electromagnetic fields. 



A. Quantization of Free Maxwell Field 



There are four Maxwell equations: 



V- E(^R,t^ = 47rp (r, , 
V-B(^R,tj =0, 

I dB [r, 



V xE 



dt 



VxB(R,t^ = ^j(^R,tj + i 



ldE[R,t 



dt 



(1) 
(2) 

(3) 

(4) 



where the Gaussian system of units have been adopted. The electric and the magnetic field are defined respectively by E = 
— V$ — c~^dtA and B = V x A, where $ is the scalar potential and A is the vector potential. Equations (1) through (4) are 
combined to give 



1=1 
+ 



di^+-dtAi + €ijkdjAk 



C & 



di^ + -dtAi 
C 



+ -^eijkdjd^Ah > e; = 



where the Einstein summation convention is assumed for repeated indices. Because the components along basis direction e; are 
independent of each other, the above vector algebraic relation becomes three equations: 



A-TT 



di<i>+^dtAi+eijkdjAk 



47r 



-dtAi 



1 



eijkdjd'tAk = 0, 



(5) 



where Z = 1, 2, 3. 

To understand the full implications of electrodynamics, one has to solve the above set of coupled differential equations. 
Unfortunately, they are in general too compUcated to solve exactly. The need to choose an appropriate gauge to approximately 
solve the above equations is not only an option, it is a must. Also, for what is concerned with the vacuum-fields, that is, the 
radiation from matter when it is in its lowest energy state, information about the charge density p and the current density J must 
be first prescribed. Unfortunately, to describe properly the charge and current densities of matter is a major difficulty in its own. 
Therefore, the charge density p and the current density J are set to be zero for the sake of simpUcity and the Coulomb gauge, 
V • A = 0, is adopted. Under these conditions, equation (5) is simplified to 



0, Z = 1,2,3. 
The steady state monochromatic solution is then of the form 

A (r, t^=a (i) Ao + a* (t) A* (e) 

= a (0) exp {-iivt) Ao (r^ + a* (0) exp {icut) A^ (r^ , 
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where (^Rj is the solution to the Heknholtz equation S/^Aq (^Rj + c '^uj'^Aq (^Rj = 0, and a (t) is the solution of the 
temporal differential equation satisfying a {t) + w^a (t) = 0. With the solution A ^R, tj , the electric and the magnetic fields 



are found to be 



and 



E{RA] = -- 



( 



: {t) Ao + a* {t) 1*0 (i? 



B (^R, tj =a{t)y X Ao (^Rj + a* {t) V x A^ (^Rj . 
The electromagnetic field Hamiltonian becomes: 

Hf = ^ j E(^R,ty E* (^R, tj+B (^R, tj ■ B* {r, tj 

= |^ii»wf, 



dV 



(6) 



An R 



dV = 1. 



where fc is a wave number, and Ao \^Rj have been normalized such that Jy 

We can transform H p into the "normal coordinate representation" through the introduction of "creation" and "annihilation" 
operators, and a. The resulting field Hamiltonian Tip of equation (6) is identical in form to that of the canonically transformed 
simple harmonic oscillator, Hsh + Ksh oh a^a. For the free electromagnetic field Hamiltonian, the canonical 

transformation is to follow the sequence Ksh cx Ija "Hsh (x E^ + B^ under a properly chosen generating function. 

The result is that with the following physical quantities, 



q{t) = ^=[a{t)-a*{t)], p{t) 
c\/47r 



[a{t)+a* (i)], 



the free field Hamiltonian of equation (6) becomes 



Tip = [p' it) + it)] 



(7) 



which is identical to the Hamiltonian of the simple harmonic oscillator. Then, through a direct comparison and observation with 
the usual simple harmonic oscillator Hamiltonian in quantum mechanics, the following replacements are made 



a it) 



a (t) , a* (t) 



CO 



at it).. 



and, the quantized relations for A (^R, tj , E (r, tj and B (r, tj are found, 

(i?, tj = [a (t) Ao (i?) + at it) A* (i?) 



A 



(8) 



E (^R, tj = iV2TThuj a (t) Ao (^Rj - a^ {t) A*o (^Rj 



(9) 



B 



a (i) V X Ao (^Rj + it) V x 1* (i? 



(10) 



where it is understood that A (^R, tj , E ^R, tj and B ^R, tj are now quantum mechanical operators. The associated field 
Hamiltonian operator for the photon is then written 



'Hp^r, 



at (t) a (i) + 



(11) 



where the hat (A) over Hf, mono now denotes an operator. 
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The generalization of the above quantization procedure for a monochromatic field to a multimode field is straightforward. We 
start by redefining the monochromatic mode function, Aq ^Rj , with the multimode counterpart, g ^ (jij , 

\k',x' = exp {ik' ■ r) A' = 1, 2, 

where V is the quantization volume; eg y , the polarization of the field mode, and the subscripts k' and A' denotes the particular 
modes of the wave. Similarly, the replacement is done for monochromatic a (t) : 

a{t) = a (0) exp {-iu)t) ag ^ (t) = ag ^ (0) exp {-icok't) , 

where the angular frequency iv in multimode representation have been replaced by Wk' to denote particular mode of wave. The 
replacements for the Hermitian conjugates, A* (R) and at ,,{t) , are straightforward. With A„ ,-, , , (R) , A*-.^ (R) 

and ag y (t) , at (t) , the monochromatic vector potential A (^R, of equation (8) is replaced by the multimode counterpart, 

the multimode vector potential A^, ^, (R,t\ : 



' 27rftc2 



Ac, X, = 



or substituting in the explicit expressions for A^ g y (^R^ , A* ^, (^R^ and ag y {t) , at [t) . 



^fe' A' 



The linearity of Maxwell's equations then allows us to write for the total vector potential in free space as 

^T{R,t) = Y^Ag y [R,t) 



or 



AT(R,t) = Y^ 



k',\' 



1 27rfe2 
Uk'V 



ag y (O)exp (i 



fe'.A' 



k' ■ R — (jJk't 



)+4',A' (O)exp(-^ 



k' ■ R — ojk't 



(12) 



Similarly, for the total electric and magnetic fields, we find 
ET[R,t^=iY^ 



2TlhuJk 



V 



ag y (O)exp(i 



k' ■ R — Wk't 



fe',A" 



(13) 



fe'.A' 



fe',A' 



'27rfe2 



fe',A' (O)exp(i 



k' ■ R — ujk't 



)+4.A' (0)exp(-^ 



k' ■ R — LOk't 



(14) 



The associated total field Hamiltonian operator for the photon is 

1 



Tip 



Stt 



dV, 



or using the fact that 



f/o,k',x' (r) ■ ^y,x' (^) dV = 4,^<5a',a, 
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along with the commutation relation 



the total field Hamiltonian operator for the photon in free space is written as 



fe',A' 

= E 

fe' A' 



(15) 



where fc' = k' , and iV^-, ^, = at (t) a^, ^, (t) is the occupation number operator. It is understood that in the free space 
limit, i.e., V ^ oo, the 3D Kronecker delta, 6^ - = <5fe^,fex'^fe^,fc^^fci,fc«j is to be replaced by the continuum counterpart, 
6^ (^k' — k^=S (fc^ — kx) S (k'y — ky) 5 {k'^ — kz) ■ Similarly, the total hnear momentum operator for the photon in free space 



IS written as 



or 



PF = ^j^^T (r, t) X Bt dV, 



Pf=Y1 

k',X' 

= E 

fe'.A' 



at (t) ac, w (t) + - 
k',\' ^ ' k',x' ^ ' ' 2 



hk' 



%A' + 



hk'. 



(16) 



The eigenvalue of the electromagnetic field energy is then 



oo 

E {^s;k',\' 




rig 


■,k',X') 


















oo / 

E {ri,;k',>: 




%',A' + I 


hck'^ 


ns]k',\'^ 


n,=0 \ 


k'.X 











2EE 



hck' } , 



where Ug = (ris; k', X 



%',A' 



ns;k\X') . The factor of two here comes from the fact that there are two possible polarizations. 



i.e., (A' = 1,2), for electromagnetic fields. With the following definitions. 



•oo oo oo 



.E[^^K^.)]^ E-EEE' 

> i=l y 711=0 712=0 713=0 



the quantized electromagnetic field energy is written as 



»la=0 



nF = 2hcY, + - E E E 



OO OO OO 



=0n2=0 713=0 \ i=l 



Y.[k'M,L,)\\ 
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Induced dipoles 




Reference origin 



Figure 2: Two interacting molecules through induced dipole interactions. 



where is the quantization length to be determined from the boundary conditions, and rij is the wave mode number for the 
corresponding fc-. The quantized field energy per quantum state 



ns;k',X ) is therefore 



^na,6 = = 



OO OO CX) 

f-^^'^' E E E 



77.1=0 712=0 7l3=0 \ i=l 



(17) 



where the subscript b of Ti.'^ ^ denotes the bounded space; and Qk' = 2, the number of polarizations. 

When the dimensions of boundaries are such that the difference, Ak'^ (rij, Lj) = fc^ (rij + l,Li) — k^ {rii, Li) , is infinitesimally 
small, we can replace the suimnation in equation (17) by integration. 



^ />oo /"OO />oo />oo /"OO /'Oc 

/ / dmdn2dnz^[h{L^)f2{L2)h{Ls)]-^ / / dk'^dk'^dk'^, 

„i=0n2=0«3=0 ^1=0^2=0^3=0 io Jo Jo 

where in the last step the functional definition for fc- = fc^ (n^, L^) = n^/i (L^) have been used to replace drii by dk[/ fi {Li) . In 
free space, the electromagnetic field energy for quantum state ns;k',X') is given by 



[hs + I] HcQk' 



^' h{L,)f2{L2)hiL,)Jo Jo Jo 



OO /■• 'DC /* 'DC' 



(18) 



where the subscript u of H'^^ „ denotes free or unbounded space, and the functional /j (L,) in the denominator is equal to 
Czeron^^ L^^ for a given Lj. Here Czero is the zeroes of the function representing the transversal component of the electric field. 



B. Casimir-Polder Interaction 



The phenomenon referred to as Casimir effect has its root in van der Waals interaction between neutral particles that are 
polarizable. The Casimir force may be regarded as a macroscopic manifestations of the retarded van der Waals force. The 
energy associated with an electric dipole moment pd in a given electric field E is Hd = —Pd ■ E. When the involved dipole 
moment pd, is that of the induced rather than that of the permanent one, the induced dipole interaction energy is reduced by a 
factor of two, TLd = —pd ■ E/2. The role of an external field here is played by the vacuum-field. Since the polarizability is 
hnearly proportional to the external field, the average value leads to a factor of one half in the induced dipole interaction energy. 
Here the medium of the dielectric is assumed to be Unear. Throughout this investigation, the dipole moments induced by vacuum 
polarization are considered as a free parameters. 

The interaction energy between two induced dipoles shown in Figure 2 are given by 



'Hint — 



R2 — Ri 



\Pd,l ■ Pd,2] 



i?2 ~ ^1 



Pd,l ■ (-R2 - ^1) Pd,2 ■ - Rl 



where Ri is the position of ith dipole. For an isolated system, the first order perturbation energy (Ti-l^^f.^ vanishes due to the 
fact that dipoles are randomly oriented, i.e., {pd,i) = 0. The first non-vanishing perturbation energy is that of the second order. 
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Uef f, static = Y<-M ) = Stot^o 1'^'"* I ™) ("^ 1'^*"* I ^) ["^0 - -^m] ^ , which falls off with respect to the separation distance 



like Uef f, static OC 



i?2 ~ ^1 



. This is the classical result obtained by F. London for short distance electrostatic fields. F. 
London employed quantum mechanical perturbation approach to reach his result on a static van der Waals interaction without 
retardation effect in 1930. 

The electromagnetic interaction can only propagate as fast as the speed of Ught in a given medium. This retardation effect 
due to propagation time was included by Casimir and Polder in their consideration. It led to their surprising discovery that the 



interaction between molecules falls off Uke 



i?i — i?2 



. It became the now well known Casimir-Polder potential [2], 



^eff,re 



he 

4-77 



i?2 ~ ^1 
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(1) (2) (1) (2) 



a(i)a(2) (1) (2)1 j 



(i) (i) 

where a^' and a\l represents the electric and magnetic polarizability of ith particle (or molecule). 

To understand the Casimir effect, the physics behind the Casimir-Polder (or retarded van der Waals) interaction is essential. In 
the expression of the induced dipole energy Hd = —Pd ■ E/2, we rewrite pd = a (ui) E^^ for the Fourier component of the dipole 
moment induced by the Fourier component of the field. Here a (w) is the polarizability. The induced dipole field energy 
becomes Hd = —a {u) E^ ■ E^/2. Summing over aU possible modes and polarizations, the field energy due to the induced 
dipole becomes 



fe.A 



where the subscripts (1) and ^1, k, Xj denote that this is the energy associated with the induced dipole moment at location 
^1 as shown in Figure 2. The total electric field E^^^ {^Ri , in mode {k, acting on pd,i is given by 

where E^^^ ^ , is the vacuum-field at location and E^^^ (^R\ , tj is the induced dipole field at due to the neigh- 
boring induced dipole pd,2 located at R2- The effective Hamiltonian becomes 

1 



fe,A 



where 



fe,A 



k,X 



fe.A 



Because only the interaction between the two induced dipoles is relevant to the Casimir effect, the Jip^ ^ ^ term is considered 
solely here. In the language of field operators, the vacuum-field E^^ ^ ^Ri,tj is expressed as a sum: 
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where 



In the above expressions, at and ar , are the creation and annihilation operators respectively; and V, the quantization volume; 



^, the polarization. By convention, e'^^J ^ (jii,tj is called the positive frequency (annihilation) operator and e'^^ J^ 
is called the negative frequency (creation) operator. 

The field operator E^^^ {R\,t^ has the same form as the classical field of an induced electric dipole, 



Pd,2 ■ ^]^} \\pd,2 (t - ^) 



+ 



Pd,2 [t- - 

C 



where r = 



i?2 ~ Ri 



S 



S 



I- ipd,2 

R2 — Ri 



i?2 — Ri 



, Pd,2 = Pd,2 \\Pd,2\\ ^ as shown in Figure 2, and c is the speed 
of fight in vacuum. Because the dipole moment is expressed as pd = a (w) E^^, the appropriate dipole moment in the above 
expression for ^2 it a (^i' to be replaced by 

Pd,2 = K) (^2'^) ' 

fe,A 

where 0:2 {tOk) is now the polarizability of the molecule or atom associated with the induced dipole moment 2 at the location 
R2. With this in place, E^^ ^ (jli , l)j is now a quantum mechanical operator. 
The interaction Hamiltoruan operator ^.p^ 2 written as 

fc,A 

where we have taken into account the fact that e'^^) , ( -R2, i I It^ac) = (vac\ E^~) , ( -R2, i I = 0. It was shown in [24] in great 

o,fe,A \ / o,K,A \ / 

detail that the interaction energy is given by 

U (r) = (Wp-^,i,pd,2) = -^Mfi E k^ojuai {ujk) 0.2 (wfe) exp (-ifcr) exp (^iA; ■ 



1 - 



fe,A 



21 1 

kr 



1 



In the limit of r ^ c |cj„i„ | ^ , where ujmn is the transition frequency between the ground state and the first excited energy state, 
or the resonance frequency, the above result becomes 



rr / \ ~ 3/1^0 2 



^ "'m|pd|0)f. 



This was also the non-retarded van der Waals potential obtained by F. London. Here is the transition frequency, and a is the 
static {uj = 0) polarizability of an atom in the ground state. Once the retardation effect due to light propagation is taken into 
account, the Casimir-Polder potential becomes. 
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What we try to emphasize in this brief derivation is that both retarded and non-retarded van der Waals interaction may be 
regarded as a consequence of the fluctuating vacuum-fields. It arises due to a non-vanishing correlation of the vacuum-fields over 



distance of r 



R2-R1 



o,k,X 



{R1: 



o.k.X 



. The non- vanishing correlation here is defined by ' 

In more physical terms, the vacuum-fields induce fluctuating dipole moments in polarizable media. The correlated dipole-dipole 
interaction is the van der Waals interaction. If the retardation effect is taken into account, it is called the "Casimir-Polder" 
interaction. 

In the Casimir-Polder picture, the Casimir force between two neutral parallel plates of infinite conductivity was found by a 
simple surmnation of the pairwise intermolecular forces. It can be shown that such a procedure yields for the force between two 
parallel plates of infinite conductivity [24] 

f [a, 1j, C) Casimir- Polder 



(19) 



6407r2d4 ' 

where it is understood that the sign of the force is attractive. When this is compared with the force of equation (20) computed 
with Casimir's vacuum-field approach, which will be discussed in the next section, the agreement is within 20% [24]. In other 
words, one can obtain a fairly reasonable estimate of the Casimir effect by simply adding up the pairwise intermolecular forces. 
The recent experimental verification of the Casimir-Polder force can be found in reference [26]. 

The discrepancy of ~ 20% between the two force results of equations (19) and (20) can be attributed to the fact that the force 
expression of equation (19) had been derived under the assumption that the intermolecular forces were additive in the sense that 
the force between two molecules is independent of the presence of a third molecule [24, 27]. The van der Waals forces are not 
however simply additive (see section 8.2 of reference [24]). And, the motivation behind the result of equation (19) is to illustrate 
the intrinsic connection between Casimir-Polder interaction and the Casimir effect, but without any rigor put into the derivation. 

It is this discrepancy between the microscopic theories assuming additive intermolecular forces, and the experimental results 
reported in the early 1950s, that motivated Lifshitz in 1956 to develop a macroscopic theory of the forces between dielectrics 
[28, 29]. Lifshitz theory assumed that the dielectrics are characterized by randomly fluctuating sources. From the assumed delta- 
function correlation of these sources, the correlation functions for the field were calculated, and from these in turn the Maxwell 
stress tensor was determined. The force per unit area acting on the two dielectrics was then calculated as the zz component of 
the stress tensor. In the limiting case of perfect conductors, the Lifshitz theory correctly reduces to the Casimir force of equation 
(20). 



C. Casimir Force Calculation Between Two Neutral Conducting Parallel Plates 

Although the Casimir force may be regarded as a macroscopic manifestation of the retarded van der Waals force between two 
polarizable charge-neutral molecules (or atoms), it is most often alternatively derived by the consideration of the vacuum-field 
energy hjio/2 per mode of frequency lo rather than from the surmnation of the pairwise intermolecular forces. Three different 
methods widely used in Casimir force calculations are presented here. They are: (1) the Euler-Maclaurin sum approach, (2) the 
vacuum pressure approach by Milonni, Cook and Goggin, and lastly, (3) the source theory by Schwinger. The main purpose 
here is to exhibit their different calculational techniques. 



1. Euler-Maclaurin Summation Approach 



For pedagogical reasons and as a brief introduction to the technique, the Casimir's original configuration (two charge-neutral 
infinite parallel conducting plates) shown in Figure 3 is worked out in detail. 

Since the electromagnetic fields are sinusoidal functions, and the tangential component of the electric fields vanish at the 
conducting surfaces, the functions fi (Li) have the form fi (Li) ~ t^L^^- The wave numbers are given by /c- (rii.Li) = 



nifi {Li) = rLiirL^ . The vacuum state is given by 
radiation energy is given by 

tic. 



0; fc', A') . Then, for rig = in equation (17), the ground state 



00 



C50 



E E 

m —0 712—0 n3 



2^2 



\ i=i 



For the arrangement shown in Figure 3, the dimensions are such that Li ^ L3 and L2 ^ L3, where {Li,L2, L3) corresponds 
to {Lx, Ly, Lz) ■ The area of the plates are given by L\X L2- The sunnmation over m and n2 can be replaced by an integration. 



f- foc foo 00 2 2 
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z=0 



z=d 



Figure 3: A cross-sectional view of two infinite parallel conducting plates separated by a gap distance of z = d. The first two lowest wave 
modes are shown. 



For simplicity and without any loss of generality, the designation of Li = L2 = L and Ls = d yields the result 

fc pco 1*00 I 2 2~ 

Here Ti'^ 5 (d) denotes the vacuum electromagnetic field energy for the cavity when plate gap distance is d. In the Umit the gap 
distance becomes arbitrarily large, the sum over ns is also replaced by an integral representation to yield 

K„6 (00) = l^i'Qfc' ^J[K?+[k'yf + [K?dk'Jk'^dk')j . 

This is the electromagnetic field energy inside an infinitely large cavity, i.e., free space. 

The work required to bring in the plates from an infinite separation to a final separation of d is then the potential energy. 



he 



lim [ — 



"3=0 



JO JO 



[k'f + [k'^]+[k'fdk'Jk'ydk 



The result is a grossly divergent function. Nonetheless, with a proper choice of the cutoff function (or regularization function), 
a finite value for U (d) can be obtained. In the polar coordinates representation (r, 6) , we define = [k'^] + [k'yj and 



dk'^dk'y = rdrdO, then 



7r/2 roc 



E 



.2 J 2 



d^ 



rdrdO 



— lim — 

d— ►00 \ TT 



d/'OO /'7r/2 /'Oo 
/ / / / o r, .-,2 



r2 + [k'XrdrdOdk'^ 



where the integration over 9 is done in the range < 6 < 7r/2to ensure fc^ > and k'y > 0. For convenience, the integration 
over is carried out first. 



E 



^ — - — rdr — lim I — 



d2 



d 

i— ►00 \ TT 



00 /"CO 



fe' =0 .'r=0 



rdrdk. 



As mentioned earlier, (7 (d) in current form is grossly divergent. It is regularized through the use of a regularization function in 

the form of / {k') = f (^yjr'^ + [k'^f^ or / (fc') = / ( + ^2^2^-2^ ^ith the condition that / [k') = 1 for k' «; fc^uto// 

and / (fc') = for k' ^ Kutoff- Mathematically speaking, this cutoff function / {k') is able to regularize the above divergent 
function. Physically, introduction of this regularization takes care of the failure at small distance of the assumption that plates are 
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perfectly conducting for short wavelengths. It is a good approximation to assume ^ l/flo, where is the Bohr radius. 

In this sense, one is inherently assuming that Casimir effect is primarily a low-frequency or long wavelength effect. Hence, with 
the regularization function substituted in U {d) above, the potential energy becomes 



uid) = ^L^e,, 



E 



r=0 



r2 + "i'^' 



rf2 



d^ 



rdr 



Um \lr^ + [Kff {\f^+wl) rdrdK 



The summation X]i^=o ^'^'^ integral in the first term on the right hand side can be interchanged. The interchange of 
sums and integrals is justified due to the absolute convergence in the presence of the regularization function. In terms of the new 
definition for the integration variables x = r^d^Tr"^ and k = /c^dTr"^, the above expression for U {d) is rewritten as 



u (d) = ^TT^L^e 



"3=0 



x + nlf [ -\/x + nl ) dx 



lim i -pr / \/x + K'^f ( — \/x + k"^] dxdn ] 

-F(0)+ ^ FM- / F{K)dn 



where 



and 



F (^^3) = ^ J _^ \lx + nlf (^yjx + nl^ dx, 

F(k) = lim ( 4 / Vx + K^f i-,\/x + KA dx] . 
d^oo\d^ J^^o \d J ) 



Then, according to the Euler-Maclaurin summation formula [30, 31], 

\ dn 



roo 

Y^Fim)- f{k) 

n3 = l •^«=0 



2 ^ ' 12 dK 720 dK^ 



for F (oo) 0. Noting that from F (k) = /J {^Vr) dr and dF (k) /dn = -2^2/ (i^) , one can find dF (0) /dn = 0, 
d^F (0) /dn^ = —4, and all higher order derivatives vanish if one assumes that all derivatives of the regularization function 
vanish at k = 0. Finally, the result for the vacuum electromagnetic potential energy U (d) becomes 

This result is finite, and it is independent of the regularization function as it should be. The corresponding Casimir force for the 
two infinite parallel conducting plates is given by 



Fz=d {d; L,c) = - 



dU {d; L, c) Shcw'^ 



dd 



1440^4 



where this is the force on plate at 2: = d due to the presence of plate at = 0. The force on plate at 2: = due to the presence of 
plate sAz = d is, like wise, given by 



Fz=o {d; L, c) 



lUOd^ 



This can be easily verified by replacing 



Rd {z = d)- Ro {z = 0) 



16 



where Rd [z = d) and {z = 0) are the position vectors for the plates located sA z = d and z = Q, respectively. Then for the 
forces -Fz=(i (d; L, c) and F^^^q (d; L, c) , we have 



Fz=d {d; L,c) = - 



dU [d- L, c) 
dRd {z = d)' 
dU {d; L, c) 



dd 



dd 



dRd [z = d) 



and 



-Pj=o {d;L,c) = - — ^ 



dU {d; L, c) 



dRo {z = 0) 
dU {d; L, c) 



dd 



dd 



It is easily shown that 



dd 



dRo {z = 0) 



dd 



dRd {z = d) dRo {z = 0) 



where d ; 



Rd {z = d) — Ro = 0) ; and, therefore the previous result Fz=d (d: L, c) = —Fz=a {d; L, c) . 



Since the electromagnetic wave has two possible polarizations, we have 6fe' =2. The force is therefore written in magnitude 



as 



Fz=did;L,c) = Fz=o{d;L,c) 



240^4^ 



(20) 



where it is understood that the sign of the force is attractive. This is the Casimir force between two uncharged parallel conducting 
plates [3]. 

It is to be noted that the Euler-Maclaurin summation approach discussed here is just one of the many techniques that can 
be used in calculating the Casimir force. One can also employ dimensional regularization to compute the Casimir force. This 
technique can be found in section 2.2 of the reference [21]. 



2. Vacuum Pressure Approach 



The Casimir force between two perfectly conducting plates can also be calculated from the radiation pressure exerted by a 
plane wave incident normally on one of the plates. Here the radiation pressure is due to the vacuum electromagnetic fields. The 
technique discussed here is due to Milonni, Cook and Goggin [32]. 

The Casimir force is regarded as a consequence of the radiation pressure associated with the zero-point energy of huj/2 per 
mode of the field. The main idea behind this techniques is that since the zero-point fields have the momentum p'^ — hk[/2^ 
the pressure exerted by an incident wave normal to the plates is twice the energy TL per unit volume of the incident field. The 
pressure imparted to the plate is twice that of the incident wave for perfect conductors. If the wave has an angle of incidence 
6 inc, the radiation pressure is 

P = FA-^ = e,nc- 

Two factors of cos Omc appear here because (1) the normal component of the hnear momentum imparted to the plate is propor- 
tional to cos 6 inc, and (2) the element of area A is increased by 1/ cos Oinc compared with the case of normal incidence. It can 
be shown then 

P = 2n cos^ 0i„c 

= 2 X - X ^huJ X X COS^ 6inc 

where the factor of half have been inserted because the zero-point field energy of a mode of energy is divided equally 
between waves propagating toward and away from each of the plates. The COS Oinc factor have been rewritten using the fact that 
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k'^ = k' ■ e-z = k' cos 6inc, where e-z is the unit vector normal to the plate on the inside, k' =u)/c and V is the quantization 
volume. 

The successive reflections of the radiation off the plates act to push the plates apart through a pressure P. For large plates 
where fc^, k'y take on a continuum of values and the component along the plate gap is k'z = n-n/d, where n is a positive integer, 
the total outward pressure on each plate over all possible modes can be written as 



Pout — 



2^ 



dk^dky, 



where Qk' is the number of independent polarizations. 

External to the plates, the allowed field modes take on a continuum of values. Therefore, by the replacement of Yl^=i 
n~^d in the above expression, the total inward pressure on each plate over all possible modes is given by 



P - -^fi 



OO pOO 



IK] 



dkLdkLdk' 



'=i=oA;=oAi=o ^[k'J'+lk'^f + ik'^f 



Both Pout and Pin are infinite, but their difference has physical meaning. After some algebraic simplifications, the difference 
can be written as 



out ^in 



'8¥ 



Qk' 



dx 



=0 vx + n^ 



OO />oo 



Jx=0 + 11? 



dxdu 



An apphcation of the Euler-Maclaurin summation formula [30, 31] leads to the Casimir's result 



out ^in 



240(^4' 

where Qk' = 2 for two possible polarizations for zero-point electromagnetic fields. 



3. The Source Theory Approach 

The Casimir effect can also be explained by the source theory of Sch winger [21, 22, 24J. An induced dipole pd in a field E 
has an energy T-Ld = —Pd ■ E/2. The factor of one half comes from the average value of an induced dipole. When there are N 
dipoles per unit volume, the associated polarization is P = Npd and the expectation value of the energy in quantum theory is 
{T^d) = — J (vd ■ d^R. Here the polarizability in pd is left as a free parameter which needs to be determined from the 

experiment. The expectation value of the energy is then 



{Hd) = -^J(pd- + ^(-) • Pd) d^R, 



where E'^^'> (^R, tj = si^^ [R, tj + si''^^ (^R, tj . Here si"^^ is the vacuum-field and E^^'' is the field due to other sources. 
Since Ey'^^ (^R, tj \vac) = {vac\ Ei~^ (^R, tj = 0, the above expectation value of the energy can be written as 

(Hd) = -\j {pd- e''+^) d^'R + c.c, (21) 

where c.c. denotes complex conjugation. From the fact that electric field operator can be written as an expansion in the mode 
functions Aa ^Rj , 

= i ^ ^/2^ [a„ it) (^) - al {t) A^ (p)] , 

Ol 

the Heisenberg equation of motion for {t) and aa,s (t) are obtained as 

fla {t) = -iLO^a^ (t) + J 1* (n) ■ Pd (R, t) d^E, 
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aa,s it) = exp (iw, [t' - t\) dt' J 1* (^) • pd [r, t') d'R, 

where Oq, s {t) is the source contribution part of (t) . The "positive frequency" or the photon annihilation part of Es^^ (^R, tj 
can then be written as 

(i?, = 27rz ^ w^lc (i?) ^ exp [t' - t]) dt' J A*^ (i?') ■ jJ,; (^R', f'^ d^R' 

= 8Tr J G(+) i') • (^', dt'd^R', 

where G'^' t, t'^ is a dyadic Green function 



(22) 



(e, R'; t, t')^^Yl '^"^^ (^) ^« (^') (*'^« - *]) • 
Equations (21) and (22) lead to the result 

(Hd) = -SttRe [ / / G\f (^R,R';t,t'^ {pd,3 (R,t) ■ Pd,z (r! ,t'^) dt' d"" R' d^ R, 

where the summation over repeated indices is understood, and E^; denotes the real part. The above result is the energy of the 
induced dipoles in a medium due to the source fields produced by the dipoles. It can be further shown that for the infinitesimal 
variations in energy, 

{SUd) = -4Mb j^^ r^,- (r, r', (pd,j (r, t) ■ Pd,i (r', t') 

X exp {iuj [t' - t]) dijdt'd^R'd^R, 
where Vij (R,R',(j]is related to g\'^^ (R,R'; t, t'j through the relation 

G(+) (^R, R'; t, t'j = —J Tij i^R, R',Lvj exp {itv [t' - t]) dtu. 



The force per unit area can then be shown to be 



-f" ('^) = ^ y [^2 - £3] (d, d, k±, tSj d'^k±doj, 

where the factor [e2 — £3] ^jj (^d, d, k± , ujj is given by 



(23) 



[£2 - £3] Tjj (d, d,k±,u?j=2 [K3 - K2] + 2K3 



K1+K3 



+ 



e^Ki - eiK3 



K1-K3 
+ e2Kz 



K2 + K3 
K2-K3 



exp {2K3d) - 1 



£314:2 - £2-^:3 



exp {2K3d) - 1 



Here = k\—c '^oj'^e (lu) and is the dielectric constant corresponding to the region i. The plate configuration corresponding 
to the source theory description discussed above is illustrated in Figure 4. 
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z=0 



z=d 



Figure 4: A cross-sectional view of two infinite parallel conducting plates. The plates are separated by a gap distance of z = d. Also, the three 
regions have different dielectric constants Si (ui) . 

The expression of force, equation (23), is derived from the source theory of Schwinger, Milton and DeRaad [21, 22]. It 
reproduces the result of Lifshitz |28, 29], which is a generalization of the Casimir force involving perfectly conducting parallel 
plates to that involving dielectric media. The details of this brief outline of the source theory description can be found in 
references [21, 24]. 



III. REFLECTION DYNAMICS 



Once the idea of physics of vacuum polarization is taken for granted, one can move forward to calculate the effective, 
temperature-averaged energy due to the dipole-dipole interactions with the time retardation effect folded into the van der Waals 
interaction. The energy between the dielectric or conducting media is then obtained from the allowed modes of electromagnetic 
waves determined by the Maxwell equations together with the boundary conditions, granted that the most significant zero-point 
electromagnetic field wavelengths determining the interaction are large when compared with the spacing of the lattice points in 
the media. Under such an assumption, the effect of all the multiple dipole scattering by atoms in the dielectric or conducting 
media is to simply enforce the macroscopic reflection laws of electromagnetic waves; and this allows the macroscopic electro- 
magnetic theory to be used with impunity in the calculation of Casimir force, granted the classical electromagnetic fields have 
been quantized. The Casimir force is then simply obtained by taking the negative gradient of the energy in space. 

In principle, the atomistic approach utilizing the Casimir-Polder interaction explains the Casimir effect observed between any 
system. Unfortunately, the pairwise summation of the intermolecular forces for systems containing large number of atoms can 
become very complicated. H. B. G. Casimir, realizing the linear relationship between the field and the polarization, devised 
an easier approach to the calculation of the Casimir effect for large systems such as two perfectly conducting parallel plates. 
This latter development is the description of the Euler-Maclaurin summation approach shown previously, in which the Casimir 
force can be found by utilizing the field boundary conditions only. The vacuum pressure approach originated by Milonni, Cook 
and Goggin [32] is a simple elaboration of Casimir's latter invention utilizing the boundary conditions. The source theory 
description of Schwinger is an alternate explanation of the Casimir effect which can be inherently traced to the retarded van der 
Waals interaction. 

Because all four approaches which were previously mentioned, (1) the Casimir-Polder interaction, (2) the Euler-Maclaurin 
summation, (3) the vacuum pressure and (4) the source theory, stem from the same physics of vacuum polarization, they are 
equivalent. The preference of one over another mainly depends on the geometry of the boundaries being investigated. For the 
type of physical arrangements of boundary configurations that are being considered in this investigation, the vacuum pressure 
approach provides the most natural route to the Casimir force calculation. The three physical arrangements for the boundary 
configurations considered in this investigation are: (1) the plate-hemisphere, (2) the hemisphere-hemisphere and (3) a sphere 
formed by brining two hemispheres together. Because the geometric configurations of items (2) and (3) are special versions of 
the more general, plate-hemisphere configuration, the basic reflection dynamics needed for the plate-hemisphere case is worked 
out first. The results can then be applied to the hemisphere-hemisphere and the sphere configurations later. 

The vacuum-fields are subject to the appropriate boundary conditions. For boundaries made of perfect conductors, the trans- 
verse components of the electric field are zero at the surface. For this simplification, the skin depth of penetration is considered 
to be zero. The plate-hemisphere under consideration is shown in Figure 5. The solutions to the vacuum-fields are that of the 

Cartesian version of the free Maxwell field vector potential differential equation V^A ^Rj — c~'^dfA ^Rj = 0, where the 
Coulomb gauge V • A = and the absence of the source $ ^p, R ^ = have been imposed. The electric and the magnetic 

field component of the vacuum-field are given hy E = —c^^dtA and B = W x A, where A is the free field vector potential. 
The vanishing of the transversal component of the electric field at the perfect conductor surface implies that the solution for E 
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y Infinity 



Figure 5: The plane of incidence view of plate-liemisphere configuration. The waves that are supported through internal reflections in the 
hemisphere cavity must satisfy the relation A < 2 _R'2 — . 



is in the form of £^ oc sin 



in (2TrX-^ 



L 



where A is the wavelength and 



is the path length between the boundaries. The 



wavelength is restricted by the condition A < 2 



i?'- 



R'l 



2^2, where R'2 and R'l are two immediate reflection points in 

the hemisphere cavity of Figure 5. In order to compute the modes allowed inside the hemisphere resonator, a detailed knowledge 
of the reflections occurring in the hemisphere cavity is needed. This is described in the following section. 



A. Reflection Points on the Surface of a Resonator 

The wave vector directed along an arbitrary direction in Cartesian coordinates is written as 



h' 


ei 


= X, 




62 


= y, 


k' 
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= z. 



Hence, the unit wave vector, k'l 



k'li^i- Define the initial position R'o for the incident wave k'l, 



i = 2 



' O.x' 



(24) 



(25) 
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Here it should be noted that R'q really has only components Tq ^ and Tq ^. But nevertheless, one can always set ^ = 
whenever needed. Since no particular wave vectors with specified wave lengths are prescribed initially, it is desirable to employ 
a parameterization scheme to represent these wave vectors. The line segment traced out by this wave vector k'l is formulated in 
the parametric form 







-1 




k'l 


k' 



(26) 



where the variable is a positive definite parameter. The restriction > is necessary because the direction of the wave 
propagation is set by fc'i . Here R'l is the first reflection point on the hemisphere. In terms of spherical coordinate variables, R'l 
takes the form 



R'i{r'„e[A'i)=r',J2^hA^ < 



i=l 



' A'n = sin 9'i cos (j)'i , 
A'l 2 = sin 6'i sin ^[ , 
- A'i,3 = cos6'i, 



(27) 



where is the hemisphere radius, 9'i and (p'l are the polar and the azimuthal angle respectively of R'l at the first reflection point. 
Notice that subscript i of r ■ denotes "inner radius" not a summation index. 

By combining equations (26) and (27), we can solve for the parameter ^i. It can be shown that 



a = ^i,p = -k'l ■ R'o + 









2 








k'l- 


R'o 


4 


-[r'f- 


R'o 



(28) 



where the positive root for have been chosen due to the restriction > 0. The detailed proof of equation (28) is given in 
Appendix A, where the same equation is designated as equation (All). 

Substituting in equation (26), the flrst reflection point off the inner hemisphere surface is expressed as 



^'i {^i,p;R'o,k'i) =Y^ 

i=l 



k'l 



k'i,i 



(29) 



where ^i.p is from equation (28). 

The incoming wave vector k'i can always be decomposed into parallel and perpendicular components, fc'^ y and k'i^±, 
with respect to the local reflection surface. It is shown in equation (A14) of Appendix A that the reflected wave vector 



k'r has the form fc' 



n' X k'. 



X n' — II n' • k'in', where the quantities \\ and a^-.-L are the reflection co- 



efficients and n' is a unit surface normal. For the perfect reflecting surfaces, a^Jl = ctr.± = 1- In component form. 



k' 



-1 i 



'^'nk'n^'n 



|n'j/c^ n";| '^i' where it is understood that n' is already normalized and Ein- 



stein summation convention is applied to the index n. The second reflection point R'2 is found then by repeating the steps done 



for R'l and by using the expression k'r = k'r/ 



R'2 = R'l + i2,pk'r = R'l + 6,P 



X k'. 



X k'. 



X n' — ar II n' • k'tn' 



xn' — OLr II n' • k'in' 



where ^2,p is the new positive definite parameter for the second reflection point. 

The incidence plane of reflection is determined solely by the incident wave k'i and the local normal n'i of the reflect- 
ing surface. It is important to recognize the fact that the subsequent successive reflections of this incoming wave will be 
confined to this particular incidence plane. This incident plane can be characterized by a unit normal vector. For the sys- 
tem shown in Figure 5, k'i = k'l and n'n'^i = —^i,pk'i — R'o- The unit vector which represents the incidence plane is 



given by n'p^i 



12^=1 ^ijkk'ijr'g i^Si, where the sunmiations over indices j and k are implicit. If the plane 



of incidence is represented by a scalar function f {x',y',z') , then its unit normal vector n'p^i will satisfy the relationship 
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Figure 6: The thick line shown here represents the intersection between hemisphere surface and the plane of incidence. The unit vector normal 
to the plane of incidence is given by n'p,i = — ||n'p,i || X)i=i fiyfcfc'ij^o./cCi- 

n'p^i oc V'/p,i {x', y', z') . It is shown from equation (A43) of Appendix A that 



i=l 



l^v[= x', 
2^u'^ = y', 



(30) 



where —00 < {z/J = .x', v'2 = y' , = z'} < 00. 

The surface of a sphere or hemisphere is defined through the relation fhemi (x' , y' , z') = [r-] — J2i=i I^il ' where r- is the 
radius of sphere and the subscript i denotes the inner surface. The intercept of interest is shown in Figure 6. The intersection 
between the hemisphere surface and the incidence plane 1 ('^[,1^25^) given by fhemi {x' ,y' ,z') — fp^\ {x',y',z') = 0. 
After substitution of fp^i (a;', y' , z') and fuemi {x' , y' , z') , we have 



i=l 



The term [r'^f can be rewritten in the form [r^]'^ = [K i] ^ 

it can be shown from equation (A51) of Appendix A that 







' 1 - 




= x', 


-[r-^f = 0, 


i = < 


2- 




= y', 






. 3- 




= z'. 


, where ^ 


— r' 
' i,x' ' 


<2 







1 ^ -1 ^ ^ ■ 



1/2 



+ Ki 



(31) 



where i = 1, 2, 3 and eijk is the Levi-Civita symbol. The result for i^l shown above provide a set of discrete reflection points 
found by the intercept between the hemisphere and the plane of incidence. 
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Figure 7: The surface of the hemisphere-hemisphere configuration can be described relative to the system origin through R, or relative to the 
hemisphere centers through R'. 







' 1 - 




= r\ sin Q\ cos (^'y , 




i = < 


2- 




= r • sin 0[ sin , 






. 3- 


K,3 


= r- cos6'i, 



Using spherical coordinate representations for the variables r'^ i,r'^2 3 ) '^^e initial reflection point R'l can be expressed 

in terms of the spherical coordinate variables (r-, 0[,(j)[) (equation (A109) of Appendix A), 



(32) 



where r ■ is the hemisphere radius, and 6[, the polar and azimuthal angle, respectively. They are defined in equations (A102), 
(A103), (A107) and (A108) of Appendix A. Similarly, the second reflection point on the inner hemisphere surface is given by 
equation (A15 1) of Appendix A: 

!1 u'2 -y = r -sin cos (/)2, 
2^^2 =r^sin^^^sin0^, (33) 
3 1^2,3 = r'i COS 0'2, 

where the spherical angles 02 and 62 are defined in equations (A143), (A144), (A148) and (A149) of Appendix A. In general, 
leaving the details to Appendix A, the iVth reflection point inside the hemisphere is, from equation (A162) of Appendix A, 



R'n {r'i, 6'^, <f>'^) = J2 ^N,i (r'i, e'^, <f>'^) Bi, 



1 - 

2 • 
3 



ri- sm ( 



C0S( 



r-sin 6'^ sin 0^, 



(34) 



where the spherical angles 6'^ and 4>'j^ are defined in equations (A158), (A159), (A160) and (A161) of Appendix A. The details 
of all the work shown up to this point can be found in Appendix A. 

The previously shown reflection points (-R'l, R'2 and R'n) were described relative to the hemisphere center. In many cases, 
the preferred choice for the system origin, from which the variables are defined, depend on the physical arrangements of the 
system being considered. For a sphere, the natural choice for the origin is its center from which the spherical variables (r'-, 6', <p') 
are prescribed. For more complicated configuration shown in Figure 7, the preferred choice for origin really depends on the 
problem at hand. For this reason, a set of transformation rules between (r^, 9', (j)') and (r^, 9, </>) is sought. Here the primed set 
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is defined relative to the sphere center and the unprimed set is defined relative to the origin of the global configuration. In terms 
of the Cartesian variables, the two vectors R and R' describing an identical point on the hemisphere surface are expressed by 

3 3 

R {ui ,2/2,1/3)=^ i/iii, RI {v[ , f2, ^) = X! ^'i^i-' 



i=l 



i=l 



where [vi, z/2, v-^) (x, y, z) , {l'[,i'2, I's) {x' , y' , z') and (ei, 62, 63) (x, ?}, z) . The vectors R and R' are connected 
through the relation R (i^i, f2, t's) = X^i=i V'T,i + t','] e'i with Rt = Yl^i=i ^T.i^-i which represents the position of hemisphere 
center relative to the system origin. As a result, we have i^i ~ ^T,i — Ci — 0. In terms of the spherical coordinate 
representation for {vi, 1/2, 1/3) and (I'l, ^2, 1^3) , we can solvefor ^ and ^. As shown from equations (BIO) and (B 12) of Appendix 
B, the result is 



= (j) (rj, 0,(1), vt,i,i^t,2) = arctan 



i^T,2 + T^i sm sm ( 
i^T 1 + r', sin 9' cos ( 



(36) 



e = e{r[,6',4>',RT) 

{vT,i + vt,2 + r[ ainO' [cos0' + sin0']} [vt,3 + r[ cosO']~^ 



arctan 



/^arrtan ( ■^T,2+r^ sine'sin0' \\ . / x / i'T.2+r'^ sing' 

I diCLdll I — 7 — : — ^7 77 I I ^ bill I dlL/Uclll 1 — 7 — : — TTT 

y y i^T, 1 H-^"^ sm V COS q> J j \ \ vt, i sin 




(37) 



where the notation ^ and d indicates that ^ and Q are expUcitly expressed in terms of the primed variables, respectively. It is to 
be noticed that for the configuration shown in Figure 7, the hemisphere center is only shifted along y by an amount of vt,2 = a, 
which leads to i'T,i^2 = 0. Nevertheless, the derivation have been done for the case where i'T,i ^ 0,i = 1, 2, 3 for the purpose 
of generalization. 



With the magnitude 



R 



1/2 



{ELi K» + ?--A^]^} ^ , where (0', = sin6''cos(/-', (61', ^i') = sin 6*' sin </>' and 
A3 (9') = cos 0', the vector R ^r-, A, A', Rt^ is given by equation (B13) of Appendix B as 

' Ai {9,^) 

r 3 ^ 3 

R(rl,A,A',RT) = lY,WT,i+r'iA'if\ J^M, A2 f ^, = sin ^ sin (38) 



i=l 



: sm y cos ( 



A2 (f)^ = sin ^ sin (j), 
A3 = cos^. 



The details of this section can be found in Appendices A and B. 



B. Selected Configurations 



Having found all of the wave reflection points in the hemisphere resonator, the net momentum imparted on both the inner 
and outer surfaces by the incident wave is computed for three configurations: (1) the sphere, (2) the hemisphere-hemisphere 
and (3) the plate-hemisphere. The surface element that is being impinged upon by an incident wave would experience the net 

momentum change in an amount proportional to Ak' inner R's,i, R's,o^ on the irmer side, and Ak' outer R's,i + o,R's,ij 
on the outer side of the surface. The quantities Ak' inner and Ak' outer are due to the contribution from a single mode of 
wave traveling in particular direction. The notation ^; R' s,oj of Ak' inner denotes that it is defined in terms of the initial 
reflection point R' s.i on the surface and the initial crossing point -R'^^o of the hemisphere opening (or the sphere cross-section). 
The notation ^; R' s,i + aR'a,ij of Ak' outer implies the outer surface reflection point. The total resultant imparted momentum 
on the hemisphere or sphere is found by sunmiing over all modes of wave, over all directions. 



/. Hollow Spherical Shell 



A sphere formed by bringing in two hemispheres together is shown in Figure 8. The resultant change in wave vector direction 
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Figure 8: Inside the cavity, an incident wave k'i on first impact point R'i induces a series of reflections that propagate throughout the entire 
inner cavity. Similarly, a wave k' i incident on the impact point R'i + aR'i, where a is the thickness of the sphere, induces reflected wave of 
magnitude • The resultant wave direction in the external region is along R'i and the resultant wave direction in the resonator is along 

—R'i due to the fact there is exactly another wave vector traveling in opposite direction in both regions. In both cases, the reflected and incident 
waves have equal magnitude due to the fact that the sphere is assumed to be a perfect conductor. 



upon reflection at the inner surface of the sphere is from the equation (C4) of Appendix CI, 



inner -^'s,! 7 -^'5,0^ — 



An-K COS 9i 



R 



8,2 



< Ginc < 7r/2, 
n = l,2,--- , 



(39) 



where 9inc is from equation (A115); Rs,i (rl,A'g i^ and Rs,2 (j'ii^'s,2^ foUow the generic form shown in the equation (CI) of 
Appendix CI, 



K,N,2 {^'s,N^ (f>'s,N) = ^^^K,N Sin<;^s,jv> 

K,N,3 Kn) =C0S9'^^^. 



(40) 



Here the label s have been attached to denote a sphere and the obvious index changes in the spherical variables 9'^ ^ and (j)'^ ^ 
are understood from the set of equations (A158), (A159), (A160) and (A161). 

Similarly, the resultant change in wave vector direction upon reflection at the outer surface of the sphere is from equation (C5) 
of Appendix CI, 



^k' outer (; R's,l + ai?'s,i) = 4 

The details of this section can be found in Appendix CI. 



ij 



COS 9 incR' s, I, 



< 9inc < 7r/2, 

n = 1,2,--- . 



(41) 
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2. Hemisphere-Hemisphere 



For the hemisphere, the changes in wave vector directions after the reflection at a point R'h,i inside the resonator, or after the 
reflection at location R'h,i + aR'h,i outside the hemisphere, can be found from equations (39) and (41) with obvious subscript 
changes, 

/ - - \ 4n7rcos6l - f < 6'i„c < 7r/2, 



Afc' outer + a-R'A,l) = 4 k'^J COsOincR'h,! 



(43) 



n = l,2,--- ; 

< e,nc < 7r/2, 
n=l,2,--- , 

where the reflection location Rh,N (r'i, ^h,N, ^'h,N^ ^T,h^ follows the generic form as shown in equation (C6) of Appendix C2, 

Rh,N (r'i,Ah,N, ^'h,N, RT,h) = < X] [^T,h,i + 'r'ih'h^N,i\ ^ \ 51 ^^KnA- 



(44) 



In the above equation, the subscript h denotes the hemisphere; and 

A?i,Ar,l (Oh,N,4'h,N^ = sinOh^N COS (j)h,N, 

A/i,jv,2 (Oh,N,4>h,N^ = smOh,N sm^h,N, 

^h,N,3 (^^h,N^ = COS0h,N- 

The expressions for A'^ ^, i = 1,2, 3, are defined identicaUy in form. The angular variables in spherical coordinates, 9h,N and 
4)h,N, can be obtained from equations (36) and (37), where the obvious notational changes are understood. The implicit angular 
variables, ^ and 4>'f^ ^, are the sets defined in Appendix A, equations (A158) and (A159) for 0'^ ^, and the sets from equations 
(A160) and (A161) for ^'^ . 

Unlike the sphere situation, the initial wave could eventually escape the hemisphere resonator after some maximum number 
of reflections. It is shown in the Appendix C2 that this maximum number for internal reflection is given by equation (C8), 

^h,max ~ {^h^maxlQ 5 (45) 

where the greatest integer function [Zh^maxja is defined in equation (C7) of Appendix C2, 



TT — arccos 



R'o 



TT — IVi 

Here ^i,p is given in equation (28) and 6inc is from equation (Al 15). 



+ r 



R'q 



(46) 



The above results of Afc' 



R' 



aR' 



h.l 



have been derived based on the fact that 



and Ak'outer \^,R'hA 

there are multiple internal refiections. For a sphere, the multiple internal refiections are inherent. However, for a hemisphere, it 
is not necessarily true that all incoming waves would result in multiple internal reflections. Naturally, the criteria for multiple 
internal reflections are in order. If the initial direction of the incoming wave vector, fc'i, is given, the internal reflections can be 
either single or multiple depending upon the location of the entry point in the cavity, R'o- As shown in Figure 9, these are two 
reflection dynamics where the dashed vectors represent the single reflection case and the non-dashed vectors represent multiple 
reflections case. Because the whole process occurs in the same plane of incidence, the vector R' f — —XqR'q where Aq > 0. 
The multiple or single internal reflection criteria can be sunmiarized by the relation found in equation (C21) of Appendix C2: 
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(47) 
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Intercept 




Figure 9: The dashed hne vectors represent the situation where only single internal reflection occurs. The dark line vectors represent the 
situation where multiple internal reflections occur 



Finally, because the hemisphere opening has a radius r^, the following criteria are concluded: 

< r-, single internal reflection, 




(48) 



where 



R' 



R'f > ''i) multiple internal reflections, 
is defined in equation (47). The details of this section can be found in Appendix C2. 



3. Plate-Hemisphere 



A surface is represented by a unit vector n'p, which is normal to the surface locally. For the circular plate shown in Figure 10, 
its orthonormal triad {n'p, 9'p, <j)'pj has the form 



n 



^ ^ OA' ■ 

' -VA' e- 0' -V— ^e- 

i=l i=l p 



6' =f ^^e. 

^sin^^a^;'" 

where A^ ^ (6*^, ^'^ = sin d'p cos 0^, A^ 2 (^'p. = sin Q'p sin 0^ and A^ 3 {ff^ = cos 6*^. 

For the plate-hemisphere configuration shown in Figure 1 1, it can be shown that the element Rp on the plane and its velocity 
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Figure 10: The orientation of a disk is given through the surface unit normal n'p. The disk is spanned by the two unit vectors 6'p and cf)' 



dRpI dt are given by (see equation (50) and (C30) in Appendix C3: 



' dO'p 
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(49) 




t^T,p,k + V\ 
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89' 



sin (9; \de'pd4>'p P d(, 




sin 0'^ dcf)'^ 



A, 




(50) 



where ^Ap^i, Ap^2, Ap,3^ is defined in equation (C31) and the angles ((ip and Op are defined in equations (C27) and (C28) of 

Appendix C3. The subscript p of (j)p and 6p indicates that these are spherical variables for the points on the plate of Figure 11, 
not that of the hemisphere. It is also understood that A^ 3 and Ap^s are independent of and (jjp, respectively. Therefore, their 
differentiation with respect to 4)'^ and 4)p respectively vanishes. The quantities O'p and (^'^ are the angular frequencies, and VT.p,i 
is the translation speed of the plate relative to the system origin. The quantities /^'p.e^ and v'p^^i^ are the lattice vibrations along 
the directions 9' p and (()' p respectively. For the static plate without lattice vibrations, v'p,e'^ and J^'p,,/,^ vanishes. 

In the cross-sectional view of the plate-hemisphere system shown in Figure 12, the initial wave vector k'i traveUng toward 
the hemisphere would go through a series of reflections according to the law of reflection and finally exit the cavity. It would 
then continue toward the plate, and depending on the orientation of plate at the time of impact, the wave- vector, now reflecting 
off the plate, would either escape to infinity or re-enter the hemisphere. The process repeats successively. In order to determine 
whether the wave that just escaped from the hemisphere cavity can reflect back from the plate and re-enter the hemisphere or 
escape to infinity, the exact location of reflection on the plate must be known. This reflection point on the plate is found to be. 
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Figure 11: The plate-hemisphere configuration. 



from equation (C54) of Appendix C3, 
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i=l 



(51) 



where the translation parameter z^T,p,i = and the terms (A^, B(^,Cc) , [A^, Bj,Cj) , {Ap, S^, C^) and 70 are defined in 
equations (C46), (C49), (C50) and (C52) of Appendix C3. It is to be noticed that for a situation where the translation parameter 
'^T,p,j = 0, A becomes identical to A' in form. Results for A can be obtained from A' by a simple replacement of primed 
variables with the unprimed ones. 

Leaving the details to the relevant Appendix, the criterion whether the wave reflecting off the plate at location Rp can re-enter 
the hemisphere cavity or simply escape to infinity is found from the result shown in equation (C58) of Appendix C3, 
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Cp^C-^A^A0 + ^oC^^C-^B^Ap + C^^B^Bp^ Aj,,^) ( ^ {a^,^ [kN^,ma.+hK,kn'p,k 



(52) 



where i = 1, 2, 3 and ^„ ^ is the component of the scale vector = Y^i=i ^K.i^i explained in the Appendix C3. 

In the above re-entry criteria, it should be noticed that Rq <r[. This implies Tq ^ <r[, where r[ is the radius of hemisphere. 
It is then concluded that aU waves re-entering the hemisphere cavity would satisfy the condition ^„ 1 = ^^^2 = Ck,3- Oh the 
other hand, those waves that escapes to infinity cannot have all three ^„ ^ equal to a single constant. The re-entry condition 
= ^K,2 = Ck,3 is just another way of stating the existence of a parametric line along the vector kr,Nh,max+i ^^^^ happens 
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Figure 12: The intersection between oscillating plate, hemisphere and the plane of incidence whose normal is n'p,i = 
II ") II ^ X"^*^ u' ' ^ 

to pierce through the hemisphere opening. In case such a line does not exist, the initial wave direction has to be rotated into a 
new direction such that there is a parametric line that pierces through the hemisphere opening. That is why all three values 
cannot be equal to a single constant. The re-entry criteria are summarized here for bookkeeping purpose: 

Cfv 1 = 2 = 3 • wave reenters hemisphere, 

(53) 

ELSE : wave escapes to infinity, 
where ELSE is the case where ^^,1 = ^k,2 = ^k,3 cannot be satisfied. The details of this section can be found in Appendix C3. 

C. Dynamical Casimir Force 

The phenomenon of Casimir effect is inherently a dynamical effect due to the fact that it involves radiation, rather than static 
fields. One of our original objectives in studying the Casimir effect was to investigate the physical implications of vacuum- 
fields on movable boundaries. Consider the two parallel plates configuration of charge-neutral, perfect conductors shown in 
Figure 13. Because there are more wave modes in the outer region of the parallel plate resonator, two loosely restrained (or 
unfixed in position) plates will accelerate inward until they finally meet. The energy conservation would require that the energy 
initially confined in the resonator when the two plates were separated be transformed into the heat energy that acts to raise the 
temperatures of the two plates. 

Davies in 1975 [33], followed by Unruh in 1976 [34], have asked the similar question and came to a conclusion that when 
an observer is moving with a constant acceleration in vacuum, the observer perceives himself to be immersed in a thermal bath 
at the temperature T = HR/ [Znck'] , where R is the acceleration of the observer and k' , the wave number. The details of the 
Unruh-Davies effect can also be found in the reference [24]. The other work that dealt with the concept of dynamical Casimir 
effect is due to Schwinger in his proposals [21, 23] to explain the phenomenon of sonoluminescense. Sonoluminescense is a 
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z=0 



z=d 



Figure 13: Because there are more vacuum-field modes in the external regions, the two charge-neutral conducting plates are accelerated inward 
till the two finally stick. 



phenomenon in which when a small air bubble filled with noble gas is under a strong acoustic-field pressure, the bubble will 
emit an intense flash of light in the optical range. 

Although the name "dynamical Casimir effect" have been introduced by Schwinger, the motivation and derivation behind the 
dynamical Casimir force in this investigation did not stem from that of Schwinger' s work. Therefore, the dynamical Casimir 
force here should not have any resemblance to Schwinger's work to the best of our knowledge. We have only found out of 
Schwinger's proposals on sonoluminescense after the work on dynamical Casimir force have already begun. The terminology 
"dynamical Casimir force" seemed to be appealing enough for its usage at the beginning of this work. After discovering 
Schwinger's work on sonoluminescense, we have found that Schwinger had already introduced the terminology "dynamical 
Casimir effect" in his papers. Our development to the dynamical Casimir force formaUsm is briefly presented in the following 
sections. The details of the derivations pertaining to the dynamical Casimir force can be found in Appendix D. 



/. FormaUsm of Zero-Point Energy and its Force 

For massless fields, the energy-momentum relation is Ti'^ = Exotai = PC, where p is the momentum, c the speed of light, 
and Ti'n^ is the quantized field energy for the harmonic fields of equation (17) for the bounded space, or equation (18) for the 

free space. For the bounded space, the quantized field energy Ti'^ = H'^ ^ of equation (17) is a function of the wave number 
k'^ (rii) , which in turn is a function of the wave mode value rii and the boundary functional fi (Li) , where L,; is the gap distance 

in the direction of Li = R'2 ■ ei — R'l ■ ei Ci. Here R'l and R'2 are the position vectors for the involved boundaries. As an 

illustration with the two plate configuration shown in Figure 13, R'l may represent the plate positioned at 2; = and R'2 may 
correspond to the plate at the position z = d. When the position of these boundaries are changing in time, the quantized field 
energy = H'^ ^ will be modified accordingly because the wave number functional k'^ (rij) is varying in time. 

Here the term proportional to hi refers to the case where the boundaries remain fixed throughout all times but the number of 
wave modes in the resonator are being driven by some active external influence. The term proportional to Li represents the 
changes in the number of wave modes due to the moving boundaries. 

For an isolated system, there are no external influences, hence hi = 0. Then, the dynamical force arising from the fact that 
the time variation of the boundaries is given by equation (D17) of Appendix Dl, 
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(54) 



where Ca,i, Ca,2, Ca,3, Ca,4, Ca,5, Cafi and Ca,7 are defined in equations (D6), (D9), (D14), (D15) and (D16) of Appendix 
Dl. ' ' ' ' ' ' 

The force shown in the above expression vanishes for the one dimensional case. This is an expected result. To understand 
why the force vanishes, we have to refer to the starting point equation (D4) in the Appendix Dl . The summation there obviously 
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X Z axis is out of the page ! 



Figure 14: A one dimensional driven parallel plates configuration. 



runs only once to arrive at the expression, dH'^^ / dk[ = [n^ + i] he. This is a classic situation where the problem has been over 
specified. For the 3D case, equation (D4) is a combination of two constraints, \p'i^ ^'^'^ • Fo"" '^^e one dimensional 

case, there is only one constraint, Ti.'^^ . Therefore, equation (D4) becomes an over specification. In order to avoid the problem 
caused by over specifications in this formulation, the one dimensional force expression can be obtained directly by differentiating 
equation (Dl) instead of using the above formulation for the three dimensional case. The ID dynamical force expression for an 
isolated, non-driven systems then becomes (see equation (D18) of Appendix Dl) 

p^ndldn^^ (55) 
cdL dk' ' ^ ^ 

where ^' is an one dimensional force. Here the subscript aofT'a have been dropped for simplicity, since it is a one dimensional 
force. The details of this section can be found in Appendix Dl . 



2. Equations of Motion for the Driven Parallel Plates 



The Unruh-Davies effect states that heating up of an accelerating conductor plate is proportional to its acceleration through 
the relation T = hR/ [2-Kck'] , where R is the plate acceleration. A one dimensional system of two parallel plates, shown in 
Figure 14, can be used as a simple model to demonstrate the complicated sonoluminescense phenomenon for a bubble subject 
to a strong acoustic field. 

The dynamical force for the ID, Unear coupled system can be expressed with equation (55), 

Ri - rjiRi - r]2R2 = ^rp, R2 - r?3-R2 - r?4-Ri = ^ip, (56) 

where the quantities 771, 772, 773, 774, ^rp, 6pi R2 are defined in equation (D31) of Appendix D2. Here i?i represents the 
center of mass position for the "Right Plate" and R2 represents the center of mass position for the "Left Plate" as illustrated 
in Figure 14. With a slight modification, equation (56) for this linear coupled system can be written in the matrix form, (see 
equations (D33), (D34) and (D35) of Appendix D2): 

J to J to 

and 
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^rp 


R4 


















Mr, 







(57) 



where 

^1 = R3, R2 = Ri, 

< R3= Ri = ^rp + VlRl + ^-^2 = ^rp + TjlRz + r?2-R4, 

_ i?4 = i?2 = S,ip + r73i?2 + = iip + V3R4 + mRz- 
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The matrix equation has the solutions given by equations (D5 1) and (D52) of Appendix D2: 
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(58) 
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(59) 



where the terms A3 and A4 are defined in equation (D37); and it, to) , V'12 it, to) , V'21 it, to) and V'22 it, to) are defined in 
equations (D43) through (D46) in Appendix D2. The quantities Rrp,cm,a and Rip,cm,a are the speed of the center of mass of 
"Right Plate" and the speed of the center of mass of the "Left Plate," respectively, and a defines the particular basis direction. 

The corresponding positions Rrp,cm,a it) and Rip,cm,a it) are found by integrating equations (58) and (59) with respect to 
time, 
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(61) 



The remaining integrations are straightforward and the explicit forms will not be shown here. 

One may argue that for the static case, Rrp,cm,a ito) and Rip,cm,a (to) must be zero because the conductors seem to be 
fixed in position. This argument is flawed, for any wall totally fixed in position upon impact would require an infinite amount 
of energy. One has to consider the conservation of momentum simultaneously. The wall has to have moved by the amount 
^Rwaii = Rwaii^t, where At is the total duration of impact, and Rwaii is calculated from the momentum conservation and it 
is non-zero. The same argument can be appUed to the apparatus shown in Figure 14. For that system 
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For simplicity, assuming that the impact is always only in the normal direction, 
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where the differences under the magnitude symbol imply field energies from different regions counteract the other. The details 
of this section can be found in Appendix D2. 
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Edge of Universe 




Sphere vacuum-field Poynting vector field lines 



Figure 15: Boyer's configuration is such that a sphere is the only matter in the entire universe. His universe extends to the infinity, hence there 
are no boundaries. The sense of vacuum-field energy flow is along the radial vector f, which is defined with respect to the sphere center. 

Edge of Universe 




Poynting vector field line from left hemisphere Sphere vacuum-field Poynting vector field lines 



Figure 16: Manufactured sphere, in which two hemispheres are brought together, results in small non-spherically symmetric vacuum-field 
radiation inside the cavity due to the configuration change. For the hemispheres made of Boyer's material, these fields in the resonator will 
eventually get absorbed by the conductor resulting in heating of the hemispheres. 

IV. RESULTS AND OUTLOOK 

The results for the sign of Casimir force on non-planar geometric configurations considered in this investigation will eventually 
be compared with the classic repulsive result obtained by Boyer decades earher. For this reason, it is worth reviewing Boyer's 
original configuration as shown in Figure 15. 

T. H. Boyer in 1968 obtained a repulsive Casimir force result for his charge-neutral, hollow spherical shell of a perfect 
conductor [4]. For simplicity, his sphere is the only object in the entire universe and, therefore, no external boundaries such 
as laboratory walls, etc., were defined in his problem. Furthermore, the zero-point energy flow is always perpendicular to his 
sphere. Such restriction constitutes a very stringent condition for the material property that a sphere has to meet. For example, 
if one were to look at Boyer's sphere, he would not see the whole sphere; but instead, he would see a small spot on the surface 
of a sphere that happens to be in his line of sight. This happens because the sphere in Boyer's configuration can only radiate in 
a direction normal to the surface. One could equivalently argue that Boyer's sphere only responds to the approaching radiation 
at normal angles of incidence with respect to the surface of the sphere. When the Casimir force is computed for such restricted 
radiation energy flow, the result is repulsive. In Boyer's picture, this may be attributed to the fact that closer to the origin 
of a sphere, the spherically symmetric radiation energy flow becomes more dense due to the inverse length dependence, and 
this density decreases as it gets further away from the sphere center. This argument, however, seems to be flawed because it 
inherently implies existence of the preferred origin for the vacuum fields. As an illustration, Boyer's sphere is shown in Figure 
15. For the rest of this investigation, "Boyer's sphere" would be strictly referred to as the sphere made of material with such a 
property that it only radiates or responds to vacuum-field radiations at normal angle of incidence with respect to its surface. 

The formation of a sphere by bringing together two nearby hemispheres satisfying the material property of Boyer's sphere is 
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A virtual photon along one of the Poynting 
vector field lines from left hemisphere 



Figure 17: The process in which a configuration change from hemisphere-hemisphere to sphere inducing virtual photon in the direction other 
than f is shown. The virtual photon here is referred to as the quanta of energy associated with the zero-point radiation. 

illustrated in Figure 16. Since Boyer's material property only allow radiation in the normal direction to its surface, the radiation 
associated with each hemisphere would necessarily go through the corresponding hemisphere centers. For clarity, let us define 
the unit radial basis vector associated with the left and right hemispheres hy tl and fji , respectively. If the hemispheres are made 
of normal conductors the radiation from one hemisphere entering the other hemisphere cavity would go through a complex series 
of reflections before escaping the cavity. Here, a conductor with Boyer's stringent material property is not considered normal. 
Conductors that are normal also radiate in directions non-normal to their surface, whereas Boyer's conductor can only radiate 
normal to its surface. Due to the fact that Boyer's conducting materials can only respond to radiation impinging at a normal 
angle of incidence with respect to its surface, all of the incoming radiation at oblique angles of incidence with respect to the 
local surface normal is absorbed by the host hemisphere. This suggests that for the hemisphere-hemisphere arrangement made 
of Boyer's material shown in Figure 16, the temperature of the two hemispheres would rise indefinitely over time. This does not 
happen with ordinary conductors. This suggests that Boyer's conducting material, of which his sphere is made, is completely 
hypothetical. Precisely because of this material assumption, Boyer's Casimir force is repulsive. 

For the moment, let us relax the stringent Boyer's material property for the hemispheres to that of ordinary conductors. For 
the hemispheres made of ordinary conducting materials, there would result a series of reflections in one hemisphere cavity due 
to those radiations entering the cavity from nearby hemisphere. For simplicity, the ordinary conducting material referred to 
here is that of perfect conductors without Boyer's hypothetical material property requirement. Furthermore, only the radiation 
emanating normally with respect to its surface is considered. The idea is to illustrate that the "normally emanated radiation" 
from one hemisphere results in elaboration of the effects of "obliquely emanated radiation" on another hemisphere cavity. Here 
the obliquely emanated radiation means those radiation emanating from a surface not along the local normal of the surface. 

When two such hemispheres are brought together to form a sphere, there would exist some radiation trapped in the sphere of 
which the radiation energy flow lines are not spherically symmetric with respect to the sphere center To see how a mere change 
in configuration invokes such non-spherically symmetric energy flow, consider the illustration shown in Figure 17. For clarity, 
only one "normally emanated radiation" energy flow line from the left hemisphere is shown. When one brings together the two 
hemispheres just in time before that quantum of energy escapes the hemisphere cavity to the right, the trapped energy quantum 
would continuously go through series of complex reflections in the cavity obeying the reflection law. But how fast or how slow 
one brings in two hemispheres is irrelevant in invoking such non-spherically symmetric energy flow because the gap S can be 
chosen arbitrarily. Therefore, there would always be a stream of energy quanta crossing the hemisphere opening with ^ ^ as 
shown in Figure 17. In other words, there is always a time interval At within which the hemispheres are separated by an amount 
S before closure. The quanta of vacuum-field radiation energy created within that time interval At would always be satisfying 
the condition ^ 0, and this results in reflections at oblique angle of incidence with respect to the local normal of the walls of 
inner sphere cavity. Only when the two hemispheres are finally closed, would then ^ = and the radiation energy produced in 
the sphere after that moment would be spherically symmetric and the reflections would be normal to the surface. However, those 
trapped quantum of energy that were produced prior to the closure of the two hemispheres would always be reflecting from the 
inner sphere surface at oblique angles of incidence. 

Unlike Boyer's ideal laboratory, realistic laboratories have boundaries made of ordinary material as illustrated in Figure 18. 
One must then take into account, when calculating the Casimir force, the vacuum-field radiation pressure contributions from 
the involved conductors, as well as those contributions from the boundaries such as laboratory walls, etc. We wiU examine the 
physics of placing two hemispheres inside the laboratory. 

For simpficity, the boundaries of the laboratory as shown in Figure 19 are assumed to be simple cubical. Normally, the 
dimension of conductors considered in Casimir force experiment is in the ranges of microns. When this is compared with the 
size of the laboratory boundaries such as the walls, the walls of the laboratory can be treated as a set of infinite parallel plates 
and the vacuum-fields inside the the laboratory can be treated as simple plane waves with impunity. 
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Poynting vector field lines originating 




Figure 18: A realistic laboratory has boundaries, e.g., walls. These boundaries have effect similar to the field modes between two parallel 
plates. In 3D, the effects are similar to that of a cubical laboratory, etc. 




Figure 19: The schematic of sphere manufacturing process in a realistic laboratory. 



The presence of laboratory boundaries induce reflection of energy flow similar to that between the two parallel plate ar- 
rangement. When the two hemisphere arrangement shown in Figure 16 is placed in such a laboratory, the result is to elaborate 
the radiation pressure contributions from obliquely incident radiations on external surfaces of the two hemispheres. If the two 
hemispheres are made of conducting material satisfying Boyer's material property, the vacuum-field radiation impinging on 
hemisphere surfaces at oblique angles of incidence would cause heating of the hemispheres. It means that Boyer's hemispheres 
placed in a realistic laboratory would continue to rise in temperature as a function of time. However, this does not happen with 
ordinary conductors. 

If the two hemispheres are made of ordinary perfect conducting materials, the reflections of radiation at oblique angles of 
incidence from the laboratory boundaries would elaborate on the radiation pressure acting on the external surfaces of two 
hemispheres at oblique angles of incidence. Because Boyer's sphere only radiates in the normal direction to its surface, or only 
responds to impinging radiation at normal incidence with respect to the sphere surface, the extra vacuum-field radiation pressures 
considered here, i.e., the ones involving oblique angles of incidence, are missing in his Casimir force calculation for the sphere. 



A. Results 



T. H. Boyer in 1968 have shown that for a charge-neutral, perfect conductor of hollow spherical shell, the sign of the Casimir 
force is positive, which means the force is repulsive. He reached this conclusion by assuming that all vacuum-field radiation 
energy flows for his sphere are spherically symmetric with respect to its center. In other words, only the wave vectors that 
are perpendicular to his sphere surface were included in the Casimir force calculation. In the following sections, the non- 
perpendicular wave vector contributions to the Casimir force that were not accounted for in Boyer's work are considered. 
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Figure 20: The vacuum-field wave vectors k'i,b and k'ij impart a net momentum of the magnitude ||pnet|| = 'i||fe'i,6 — /2 on 

differential patch of an area dA on a conducting spherical surface. 



/. Hollow Spherical Shell 



As shown in Figure 20, the vacuum-field radiation imparts upon a differential patch of an area dA on the inner wall of the 
conducting spherical cavity a net momentum of the amount 



^Pinner — ^^^^^inner -^'s,! ; -^'5,0^ — 



2m:h cos 9j. 
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where Ak'mner R's,i,R's,o^ is from equation (39). The angle of incidence 6inc is from equation (Al 15); Rs,i ^r-, 1^ 

and Rs,2 (r'i, ^'s,2^ follow the generic form shown in equation (40). 

Similarly, the vacuum-field radiation imparts upon a differential patch of an area dA on the outer surface of the conducting 
spherical shell a net momentum of the amount 
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where /\k' outer R' s.i + is from equation (41). 

The net average force per unit time, per initial wave vector direction, acting on differential element patch of an area dA is 
given by 
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Notice that Ps,avg is called a force per initial wave vector direction because it is computed for k'i^b and k'ij along specific 
initial directions. Here /c'i b denotes a particular initial wave vector k'i entering the resonator at J?s as shown in Figure 8. 
The subscript b for k'i^b denotes the bounded space inside the resonator. The k'ij denotes a particular initial wave vector k'i 
impinging upon the surface of the unbounded region of sphere at point R' s,\ + aR's,i as shown in Figure 8. The subscript / for 
k'ij denotes the free space external to the resonator. 
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Because the wave vector k'ij resides in free or unbounded space, its magnitude 
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Umit is the case where the radius of the sphere becomes very large. Therefore, by designating 
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and summing over all allowed modes, the total average force per unit time, per initial wave vector direction, per unit area is 
given by 
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In the limit — > oo, the second sunnmation to the right can be replaced by an integration, X^^^ — > /q°° dn. Hence, we have 
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the total average force per unit time, per initial wave vector direction, per unit area is written as 
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(62) 



where < Omc < 7r/2 and n = 1,2,3, - ■ ■ . The total average vacuum-field radiation force per unit time acting on the uncharged 
conducting spherical shell is therefore 
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(63) 



where dSgphere is a differential surface element of a sphere and the integration is over the spherical surface. The term R's,o 
is the initial crossing point inside the sphere as defined in equation (25). The notation X^jp , ^ g. ^ gj imply the sunnmation 

over all initial wave vector directions for both inside ^^'^,6^ and outside {k'ij^ of the sphere, over aU crossing points given by 

It is easy to see that Ts,avg of equation (62) is an "unregularized" ID Casimir force expression for the parallel plates (see 
the vacuum pressure approach by Milonni, Cook and Goggin [32]). It becomes more apparent with the substitution At = d/ c. 
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An application of the Euler-Maclaurin summation formula [30, 31] leads to the regularized, finite force expression. The force 



3^s,avg is attractive because 
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where 



is a constant for a given initial wave k'i^b and the initial crossing point R's,q in the 

cross-section of a sphere (or hemisphere). The total average force Fg^ioLai, which is really the sum of Ts.avg o'^^r all R' sfi and 
all initial wave directions, is therefore also attractive. For the sphere configuration of Figure 8, where the energy flow direction is 
not restricted to the direction of local surface normal, the Casimir force problem becomes an extension of infinite set of parallel 
plates of a unit area. 



2. Hemisphere-Hemisphere and Plate-Hemisphere 



Similarly, for the hemisphere-hemisphere and plate-hemisphere configurations, the expression for the total average force per 
unit time, per initial wave vector direction, per unit area is identical to that of the hoUow spherical shell with modifications. 
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where 6inc < 7r/2 and n= 1, 2, 3, • • • . The incidence angle Oinc is from equation (A115); Rh,i (r'i, AJ^ and Rh,2 (j'ij^h,2j 
follow the generic form shown in equation (44). This force is attractive for the same reasons as discussed previously for the 
hollow spherical shell case. The total radiation force averaged over unit time, over all possible initial wave vector directions, 
acting on the uncharged conducting hemisphere-hemisphere (plate-hemisphere) surface is given by 



.total 



E 



,R'h.o} 



2mrh cos 6r 



n=l ll-R, 



M {<,K,2)-Rh,i {<,K,i) 



_2h 

cos fine 

TT 



X lim 



Rh,2 (r;, 4,2) - Rh,i {r'i,Ki) II K 



(65) 



where dShemisphere IS now a differential surface element of a hemisphere and the integration is over the surface of the 

hemisphere. The term R'hfi is the initial crossing point of the hemisphere opening as defined in equation (25). The notation 

Sjfe' b fc' / J?'h 0} ™P^y sunnmation over all initial wave vector directions for both inside {k'i.bj and outside {k'i,f^ of 

the hemisphere-hemisphere (or the plate-hemisphere) resonator, over all crossing points given by R'hfi- 

It should be remarked that for the plate-hemisphere configuration, the total average radiation force remains identical to that of 
the hemisphere-hemisphere configuration only for the case where the gap distance between plate and the center of hemisphere 

is more than the hemisphere radius r-. When the plate is placed closer, the boundary quantization length 

carefully to be either 



must be chosen 



R 



h,2 



{rlA'f,^,)-Rh,i (rU'/.,i) 



or 



Rp (^i:Ap) -Rh,Nk,^^^ (rl,A'h, 



They are illustrated in Figure 12. The proper one to use is the smaller of the two. Here Rp ^r^, A^^ is from equation (C54) of 
Appendix C3 and N^^max is defined in equation (C8) of Appendix C2. 
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B. Interpretation of the Result 

Because only the specification of boundary is needed in Casimir's vacuum-field approach as opposed to the use of a polar- 
izability parameter in Casimir-Polder interaction picture, the Casimir force is sometimes regarded as a configurational force. 
On the other hand, the Casimir effect can be thought of as a macroscopic manifestation of the retarded van der Waals interac- 
tion. And the Casimir force can be equivalently approximated by a summation of the constituent molecular forces employing 
Casimir-Polder interaction. This practice inherently relies on the material properties of the involved conductors through the use 
of polarizability parameters. In this respect, the Casimir force can be regarded as a material dependent force. 

Boyer's material property is such that the atoms in his conducting sphere are arranged in such manner to respond only to the 
impinging radiation at local normal angle of incidence to the sphere surface, and they also radiate only along the direction of 
local normal to its surface. When the Casimir force is calculated for a sphere made of Boyer's fictitious material, the force is 
repulsive. Also, in Boyer's original work, the laboratory boundary did not exist. When Boyer's sphere is placed in a realistic 
laboratory, the net Casimir force acting on his sphere becomes attractive because the majority of the radiation from the laboratory 
boundaries acts to apply inward pressure on the external surface of sphere when the angle of incidence is oblique with respect 
to the local normal. If the sphere is made of ordinary perfect conductors, the impinging radiation at oblique angles of incidence 
would be reflected. In such cases the total radiation pressure applied to the external local-sphere-surface is twice the pressure 
exerted by the incident wave, which is the force found in equation (63) of the previous section. However, Boyer's sphere cannot 
radiate along the direction that is not normal to the local -sphere -surface. Therefore, the total pressure applied to Boyer's sphere 
is half of the force given in equation (63) of the previous section. This peculiar incapability of emission of a Boyer's sphere 
would lead to the absorption of the energy and would cause a rise in the temperature for the sphere. Nonetheless, the extra 
pressure due to the waves of oblique angle of incidence is large enough to change the Casimir force for Boyer's sphere from 
being repulsive to attractive. The presence of the laboratory boundaries only act to enhance the attractive aspect of the Casimir 
force on a sphere. The fact that Boyer's sphere cannot irradiate along the direction that is not normal to the local-sphere-surface, 
whereas ordinary perfect conductors irradiate in all directions, impUes that his sphere is made of extraordinarily hypothetical 
material, and this may be the reason why the repulsive Casimir force have not been experimentally observed to date. 

In conclusion, (1) the Casimir force is both boundary and material property dependent. The particular shape of the conductor, 
e.g. sphere, only introduces the preferred direction for radiation. For example, radiations in direction normal to the local surface 
has bigger magnitude than those radiating in other directions. This preference for the direction of radiation is intrinsically con- 
nected to the preferred directions for the lattice vibrations. And, the characteristic of lattice vibrations is intrinsically connected 
to the property of material. (2) Boyer's sphere is made of extraordinary conducting material, which is why his Casimir force is 
repulsive. (3) When the radiation pressures of all angles of incidence are included in the Casimir force calculation, the force is 
attractive for a charge-neutral sphere made of ordinary perfect conductor. 

C. Suggestions on the Detection of Repulsive Casimir Force for a Sphere 

The first step in detecting the repulsive Casimir force for a spherical configuration is to find a conducting material that most 
closely resembles the Boyer's material to construct two hemispheres. It has been discussed previously that even Boyer's sphere 
can produce attractive Casimir force when the radiation pressures due to obUque incidence waves are included in the calculation. 
Therefore, the geometry of the laboratory boundaries have to be chosen to deflect away as much as possible the oblique incident 
wave as illustrated in Figure 2 1 . Once these conditions are met, the experiment can be conducted in the region labeled "Apparatus 
Region" to observe Boyer's repulsive force. 

D. Outlook 

The Casimir effect has influence in broad range of physics. Here, we list one such phenomenon known as "sonoluminescense," 
and, finally conclude with the Casimir oscillator. 

1. Sonoluminescense 

The phenomenon of sonoluminescense remains a poorly understood subject to date |35, 36]. When a small air bubble of 
radius ^ 10^'' cm is injected into water and subjected to a strong acoustic field of ^ 20 kHz under pressure roughly ^ 1 oim, 
the bubble emits an intense flash of light in the optical range, with total energy of roughly ~ 10'' eV. This emission of light 
occurs at minimum bubble radius of roughly 10~^ cm. The flash duration has been determined to be on the order of 100 ps 
[37-39]. It is to be emphasized that smaU amounts of noble gases are necessary in the bubble for sonoluminescense. 
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Figure 21: To deflect away as much possible the vacuum-field radiation emanating from the laboratory boundaries, the walls, floor and ceiling 
are constructed with some optimal curvature to be determined. The apparatus is then placed within the "Apparatus Region." 



Original Bubble 




Figure 22: The original bubble shape shown in dotted lines and the deformed bubble in solid line under strong acoustic field. 



The bubble in sonoluminescense experiment can be thought of as a deformed sphere under strong acoustic pressure. The 



dynamical Casimir effect arises due to the deformation of the shape; therefore, introducing a modification to Z21 



R2 — Ri 



from that of the original bubble shape. Here L21 is the path length for the reflecting wave in the original bubble shape. In 
general L21 = L21 (t) = R2 (r^ (t) ,9 (t) ,(f> (t)) — Ri {ri (t) ,9 (t) ,(f) (t)) . From the relations found in this work for the 

reflection points .Ri (r^ (t) ,9{t) ,4> {t)) and Rn {ri (t) , 9 (t) , (j) (t)) , together with the dynamical Casimir force expression of 
equation (54), the amount of initial radiation energy converted into heat energy during the deformation process can be found. 
The bubble deformation process shown in Figure 22 is a three dimensional heat generation problem. Current investigation seeks 
to determine if the temperature can be raised sufficiently to cause deuterium-tritium (D-T) fusion to occur, which could provide 
an alternative approach to achieve energy generation by this D-T reaction (threshold ~ 17 KeV) [40]. Its theoretical treatment 
is similar to that discussed on the ID problem shown in Figure 13. 



2. Casimir Oscillator 



If one can create a laboratory as shown in Figure 21, and place in the laboratory hemispheres made of Boyer's material, then 
the hemisphere-hemisphere system will execute an oscillatory motion. When two such hemispheres are separated, the allowed 
wave modes in the hemisphere-hemisphere confinement would no longer follow Boyer's spherical Bessel function restriction. 



R2 — Ri 



, where Ri and R2 are two neighboring reflection 



Instead it will be strictly constrained by the functional relation of 

points. Only when the two hemispheres are closed, would the allowed wave modes obey Boyer's spherical Bessel function 
restriction. 

Assuming that hemispheres are made of Boyer's material and the laboratory envirormient is that shown in Figure 21, the two 
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Figure 23: The vacuum-field radiation energy flows are shown for closed and unclosed hemispheres. For the hemispheres made of Boyer's 
material, the non-radial wave would be absorbed by the hemispheres. 

closed hemispheres would be repulsing because Boyer's Casimir force is repulsive. Once the two hemispheres are separated, 
the allowed wave modes are governed by the internal reflections at oblique angle of incidence. Since the hemispheres made of 
Boyer's material are "infinitely unresponsive" to oblique incidence waves, all these temporary non-spherical symmetric waves 
would be absorbed by the Boyer's hemispheres and the hemispheres would heat up. The two hemispheres would then attract 
each other and the oscillation cycle repeats. Such a mechanical system may have application. 
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APPENDICES ON DEIUVATION DETAILS 

The appendices contain our original derivations and developments that are too tedious and lengthy to be included in the main 
body of the paper. There are four appendices: (1) Appendix A, (2) Appendix B, (3) Appendix C and (4) Appendix D. The 
appendices C and D are further divided into subparts C.l, C.2, C.3, D.l and D.2. The title and the layout of the appendices 
closely follow the main body of the paper. 

Appendix A: REFLECTION POINTS ON THE SURFACE OF A RESONATOR 

In this appendix, the original derivations and developments pertaining to the reflection dynamics used in this paper are in- 
cluded. It is referenced by the text to supply all the details. 



For the conflguration shown in Figure 5, the wave vector directed along an arbitrary direction in Cartesian coordinates is 
written as 

I — \ ^ k-^ 6i = 

k'l {k[,.,k[^y,kl,) = ^/cl,,ei, kl, = 1 1^2^ k^^, = V (Al) 
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The unit wave vector is given by 



k'l = 



(A2) 



The initial crossing position R'q of hemisphere opening for the incident wave k'l is defined as 



* = 1 ^ r' 



, i = 3 — > J,. 



(A3) 



It should be noticed here that -R'o has only two components, Vq ,^ and Tq ., . But nevertheless, one can always set r'g y = whenever 
needed. Since no particular wave with certain wavelength is prescribed initially, it is desirable to employ a parameterization 
scheme to represent these wave vectors. The Une segment traced out by the wave vector k'l is formulated in the parametric form 







-1 


r'o,i + 6 


k'l 


ki,i 



(A4) 



where the real variable is a positive definite parameter. The restriction > is a necessary condition since the direction of the 
wave propagation is set by fc'i . Here R'l denotes the first reflection point on the hemisphere. In terms of spherical coordinates, 
R'l takes the form 



i=l 



' A'li = sin 9[ cos , 
A[ 2 = sin^i sint^i, 

I A'i,3 = cos6'i, 



(A5) 



where r ■ is the hemisphere radius, 6'i and are the polar and azimuthal angles respectively of the first reflection point R'l. The 
subscript i of r ■ denotes "irmer radius" and it is not a summation index. Equations (A4) and (A5) are combined as 



3 r 



E 



r'oA + 6 



k'l 



k'i,i - r'iAi,i 



0. 



(A6) 



Because the basis vectors ii are independent of each other, the above relations are only satisfied when each coefficients of 
vanish independently. 



r'oA + a 



k'l 



k'i.-r'iAi,i = 0, i = 1,2,3. 



The three terms Ai^i=i, Ai_j=2 and Ai^i=3 satisfy an identity 

3 



(A7) 



(A8) 



From equation (A7), Af ^ is computed for each i : 

^l = {rr'\K,f+ei 



k'l 



k'l 



k'.^A, z = 1,2,3. 



Substituting the above result of Af ^ into equation (A8) and after rearrangement, one obtains 



]'-\r'A' = 0. 



(A9) 
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Surface 



Figure 24: A simple reflection of incoming wave fc'i from the surface defined by a local normal n'. 
Further simplifying, it becomes 



=0. 



There are two roots, 



and 
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^ [fc'i • i?'o 


+ [r'f- 


R'o 


y \k\ ■ R'o 


2 „ 


R'o 



ii,b = -k'i- R'o 

The root to be used should have a positive value. For the wave reflected within the hemisphere, R'o < r'^, 
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k'l ■ R'o 



where the equaUty 



k'l ■ R'o 



(AlO) 



= —k'l - R'o happens when k'l ■ R'o < 0. Therefore, ^i,a < and ^i^ > 0; the positive root ^i,6 



should be selected. For bookkeeping purposes, ^i^ is designated as ^i,p : 



= -k'l ■ R'o 



k'l ■ R'o 



R'o 



Using this positive root, the first reflection point of the inner hemisphere is found to be 



R'i{^i,p;R'o,k'i) =Y1 
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(All) 



(A12) 



The incident wave k'i, shown in Figure 24 and where i here stands for incident wave, can always be decomposed into 
components parallel and perpendicular to the vector n' normal to the reflecting surface. 



k'i = k'i.n + k'iA = i-^n' + 



If the local normal n' is already normalized to unity, the above expression reduces to 







n' X k'i 


X n' 







k'i = n' ■ k'iTi' 



n' X fc'. 



(A13) 



Here the angle between k'i and n' is it — 6i. The action of reflection only modifies fc'j^n in the reflected wave. The reflected wave 
part of k'i in equation (A13) is 



k'r = (Xr,±k'i^± — ar^\\k'i^\\ = ar,i_ n' X k'i x n' — ar^\\n' ■ k'iii' 



(A14) 
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where fc'^ y have been rotated by 180" on the plane of incidence. The new quantities ar^\\ and ar,± are the reflection coefficients. 
For a perfect reflecting surfaces, a^.n = ctr.± = 1- Because of the frequent usage of the component for k'r, equation (A14) is 



also written in component form. The component of the double cross product 



n' X k'i 



X n' is computed first, 



I ^n' X k'i X n' I = e; 



X fc'. 



= [SnqSlr - SnrSlq] u'^k^^^u'^ = 6nq5lrn'gk[ ,^n'^ - 6nr5lqn'qk[ ,^n'^ 



or 



n' X k'i 



X w 



1=1 



where the summation over the index n is implicit. In component form, k'r is hence expressed as 

3 



(A15) 



(A16) 



where it is understood n' is already normalized. 

The second reflection point R'2 is found by repeating the steps done for R'l, 



R'2 = R'l + ^2,pk'r = R'l + 6,i 





n' X k'i 


X n' — Q!r, 
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i'iU' 




n' X k'i 


X n' — a^, 


\n'-) 





where ^2,p is the new positive parameter corresponding to the second reflection point. The procedure can be repeated for any 
reflection point. Although this technique is sound, it can be seen immediately that the technique suffers from the lack of elegance. 
For this reason, the scalar field technique will be exclusively used in studying the reflection dynamics. For a simple plane, the 
scalar field function can be inferred rather intuitively. However, for more complex surfaces, one has to work it out to get the 
corresponding scalar field. For the purpose of generalization of the technique to any arbitrary surfaces, we derive the scalar field 
functional for the plane in great detail. 

In simple reflection dynamics, there exists a plane of incidence in which all reflections occur. The plane of incidence is 
determined by the incident wave k'i and the local surface normal n',. For the system shown in Figure 5, k'i and n'i are given by 

k'i = k'x, n'„j,i = -R'l ^'0, fc'i) = -^i,pk'i - R'o- 
The normal to the incidence plane is characterized by the cross product, 

3 

n'p,i = k'l X n'n'.,i = k'l x -^i,pk'i - R'o = -k'l x R'o = - ^i3kK,f'o,k^u 



where the summations over the indices j and k are imphcit. The normal to the incidence plane is normaUzed as 

.. , 3 



n'p,i = - 



(AI7) 



In order to take advantage of the information given above, the concept of scalar fields in mathematical sense is in order A 
functional / (a;', y', z') is a scalar field if to each point (x' , y' , z') of a region in space, there corresponds a number A. The study 
of a scalar field is a study of scalar valued functions of three variables. Scalar fields are connected to its normals, e.g., equation 
(A17), through the relation 



(A18) 



2'p,i a V7p,i {x',y',z') ^^e^—fp^i {x',y',z') , 
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Introducing a constant proportionality factor (3p^i, equation (A18) becomes 



d 



(A19) 



The proportionality factor /3p,i is intrinsically connected to the normalization of V'/p,i- Because the vector n'p,i is a unit vector, 
its magnitude squared is 



d 



jfp,! {x',y',z') 



= 1 - Pp,,=±{Y, 



_d 



jfpA ix',y',z') 



-1/2 



In equation (A19), the directions for vectors n'p^i and V'/p,i are intrinsically built in. Therefore, the proportionality factor Pp^i 
has to be a positive quantity. 



_d 



jfp,! {x',y',z') 



-1/2 



(A20) 



The exact form of the proportionality coefficient Pp^i requires the knowledge of /p,i, which is yet to be determined. However, 
we can use it formally for now until the solution for fp^i is found. 

Substituting the gradient function V'/p,i into equation (A19), and using equations (A17) and (A19), one arrives at 
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(A21) 



Because the basis vectors are linearly independent, the equation for each component is obtained as 



Integrating both sides of equation (A22) over the variable v'^ = a, 
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(A22) 



ra Q r" -1 

J an J an 



where the dummy variable a' is introduced for integration purpose. The terms ea'jkk'ijVQ ^, /3p,i and 
of the dummy variable a', and they can be moved out of the integrand, 

r" d 



n'p,i 



J ao 

Because the total differential of /p,i is given by 



n'p,i 



pa 

^ajkk'ijr'o^k / da'. 

Jan 



are independent 



(A23) 



dfp,^ = ^da' + ^dl3+^dj, 



(A24) 



the term [dfp^i/da'] da' can be written as 



dfp,i ,„ dfp^i 



The integration over the variable v'^ = a in equation (A23), with variables v'^^ a fixed, can be carried out with 

d 

dp = d-y ^ 0, g^fpA (o^'' 7) da' = dfp^i (a', /3, 7) 



(A25) 
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as 



to give 



The two terms 



/ rf/p,i (a',/3,7) = -/3p_i n'p,i eajkk'i^o^k j da' 

J an J an 



/p,i (a, /3, 7) = "'p,i <^ajkk'i,jr'o,k ["0 - a] + /p,i (ao, 7) - 



(A26) 
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ao and /pj (q;o,/3,7) are independent of a. These terms can only assume values 



of f ■ = /? or f ■ = 7. By re-designating a independent terms, 



equation (A26) can be rewritten for bookkeeping purposes as 

/p,i (a, /?, 7) = 7) - Pp^i 



n'p,i 



^ajkk'ijr'Qi.a. 



(A27) 



(A28) 



Substituting the result into equation (A22), a differentiation with respect to the variable v'^ = /3 gives 



d_ 

dp 



/ip,i(/3,7)-/3p"l 



n'p,i 



^ajkkijrQi^a 



n'p,i 



or 



— Vi(/3,7) = -/3p,i 



The integration of both sides with respect to the variable I'l = P yields the result 

d 



/ ^Vi(/3',7)d/3' = -/ 



e«fc^ij''o,fc / dp' 



have been taken 



where the dummy variable P' is introduced for integration purpose and the terms e/jjfeA;^ fe, /9p,i and n'p,i 

out of the integrand because they are independent of /?'. Following the same procedure used in equations (A24) through (A25), 
the integrand [dhp^i/dp'] dP' on the left hand side of the integral is 

-^hp,i{P',^)dP' = dhp,^ (/3',7). 

Consequently, hp^i {P, 7) is given by 

hp,i {P,l) = Pp^i ' ^P3kk'^A,k [/^o -P] + hp,i (/?o,7) • (A29) 

^jSjkk'i^jv'Q f./ |/3p,i «'p,i II |j Po and /ip,i (/3o, 7) are independent of p. The p independent terms can be re- 



The two terms 
designated as 



£f3jkk'ijr'o^kPo + /ip,i (/3o,7) • 



9p,i (7) = Pp,l 
For bookkeeping purposes, equation (A29) is rewritten as 

-1 

hp,i (P, 7) = 5p,i (7) - Pp} n'p.i e/3jfcfci^^-ro_fe/3. 
Substitution of hp^i (/?, 7) into equation (A28) gives 



fp,i {a, P, 7) = 9p,i (7) - 0p,l 



n'p,i 



[^0jkk'ijr'o^k0 + eajkk'ijr'Q^k<^] . 



(A30) 



(A31) 



(A32) 
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Substituting fp, i {a, 0, 7) into equation (A22) once more, and performing the differentiation with respect to the variable f • = 7, 
where 7 7^ a 7^ /3, we obtain 



d_ 



^ijkk'i^jTQf. — 



or 



-^9p,i (7) = 



where the differentiation have been changed from dtod because gp^i is a function of single variable. The integration of both 
sides with respect to the variable 1^1 = j then gives 



J^5.,i(7')d7' 



70 



70 



have 



where the dummy variable 7' have been introduced for integration purpose and the terms ^-yjkk'ijrQ^i^, Pp,i and 

been taken out of the integrand because they are independent of 7'. Knowing [dgp^i/dj'] d-y' = dgp^i, the integration is carried 
out to yield 



9p,i (7) = Pp,l n'p^i ^jjkk[yo,k [70 - 7] + 9p,i (70) • 



(A33) 



The two terms 
as 



^7jfc^ij^o,fe/ '^'j''^!}] 7o and (70) are independent of 7. The 7 independent terms are re-designated 



'^ijkk'ijr'o^klo + gp,i ilo) ■ 



n'p^i 



^7ife^i,i''o,/c7- 



60 = Pp,l n'p,i 
For bookkeeping purposes, equation (A33) is rewritten as 

9p,i (7) =bo- P~l 
Substituting gp^i (7) in equation (A32), the result for fp^i {a, (3, 7) is found to be 

!1 ^ f'l = a = x', 
3^1/^=7 = 2'. 

The cross product expressed in terms of the Levi-Civita symbol is expanded to give 

^x'jkk'ijTQ,. = k[j^y,rQf.^^, — k[f.^^,r'Qj^y,, 



(A34) 
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(A36) 



(A37) 



Sy'jkk'ijrQi^ — k[j^^,rQi^^^, k'i,.^^,rQi.^ 



(A38) 



(A39) 



It is important to understand that the functional /p,i in equation (A36) is a scalar field description of an infinite family of parallel 
planes characterized by the normal given in equation (A17), 



i=l 



Because the normal n'p.i is a cross product of the two vectors fc'i and i?'o, the surface represented by the scalar field /p i is a 
plane spanned by all the scattered wave vectors. The graphical plot of the functional /p,i is illustrated in Figure 25. The three 
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Figure 25: Parallel planes characterized by a normal n'j,,i = — ki'p.i Z)i=i ^ijkk'ijrQi^ei 



coefficients, 



^ajkk'ijrQf,, £0jkk[jrQ,., ^jjkk'ijrQ,., 



of the independent variables a, /3 and 7 define the slopes along the respective bases a, $ and 7. The integration constant 60 
provides infinite set of parallel planes whose common normal is n'p 1. For what is concerned with here, only one of them 
containing the coordinate origin is required. It is convenient to choose the plane with 60 = 0. With the plane of 60 = chosen, 
the scalar field of equation (A36) is rewritten for the sake of bookkeeping purposes: 
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. 3- 
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(A40) 



where —00 < {a, /3, 7} < 00. 

As mentioned before, /p i of equation (A40) is a scalar field where k'l and R'o are the two initially known vectors which 
span locally the reflection surface. Other than k'l and R'o, any member vectors of the spanning set for the incidence plane can 
also be used to determine the orientation of a surface. Then it is always true that 



fp,i (a,/?, 7; k'i,R'oj = fp,2 (a,f3,j; k'2,R'ij =■■■ = fp^N (a,P,^;k'N,R'N^ 



(A41) 



where the integer index N of k'^ is used to enumerate the sequence of reflections; and the integer index N of R'n is used to 
enumerate the sequence of reflection points. 

What is still yet undetermined in equation (A40) is the proportionaUty factor /3p,i. From equation (A20), the (3p^i has an 
algebraic definition 



— /p,i(a,/3,7) 



-1/2 



' 1 - 




= a, 


< 2- 




= /?, 


, 3- 




= 7. 
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The partial derivatives dfp^i/dy[ are calculated with the solution of fp^\ in equation (A40), 

3 



d 



n'p,i 



-1 Q 



which reduces to 



A 



-1/2- 



(A42) 



The two possible solutions are either /3p_i = 0, or the term enclosed in the outermost square bracket must vanish. Because 



^ijkk'i jr'Qu 



1/2 



the 



/3p,i = is a useless trivial solution, the second solution has to be adopted. Since 

equation (A42) is already satisfied for any value of /3p,i. Therefore, /3p,i is an arbitrary quantity, and it is simply chosen to be 
unity which makes the gradient function V'/p,i automatically normalized. The scalar field solution of equation (A40) is 
now restated with (ip^i = 1, 



/p,i (a,/3,7) = - n'pA ^eijkKjr'o^k^'i, i=< 



i=l 



' 1 - 




= a 


= x' 


< 2- 




= /? 


= y' 


. 3- 




= 7 


= z' 



(A43) 



where — oo < {v'l = a = x' ,1/2 = P = y' ,1^3 = 'J = z'} < 00. 

The intercept between the plane of incidence, defined by equation (A43), and the hemisphere is found through the algebraic 
relation 



which can be written as 

















[r'f, 










' 1 - 
















r'f - E W'f = 0, 


i = < 


2- 




= y' 


i=l 














. 3- 




= z' 



(A44) 



where the r\ is the radius of sphere and the index i denotes the radius of the inner surface. The intercept of interest is shown in 
Figure 6. 

One may be tempted to incorporate the surface vibration into equation (A44) through a sUght modification 



where the vibration have been introduced through the time variations in radius. Some have employed such a model in describing 

the "Casimir radiation," as well as the phenomenon of sonoluminescense mainly due to its simplicity from the mathematical 
point of view [21, 36]. In general, if one wishes to incorporate the vibration of a surface, the description of such system could 
be represented in the form 



[ra0',<^',t)]'-EK'(^',<^')f =0. 



i=l 



Since the radius function varies with 6*', (/)' and t, its treatment has to be postponed until the surface function can be found in 
later sections. In the present discussion, the hemisphere is regarded as having no vibration. 

Returning from the above short digression, the surface function of the sphere is expressed as the null function from equation 
(A44), 



V\ = X 



^2 = y 



(A45) 
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The intersection between the two surfaces, the plane of incidence defined in equation (A43) and the hemisphere defined in 
equation (A45), is found through the relation 

(x'. y'. z') - (x', y', z') = 0, (A46) 

which is equivalent to setting the scalar function fp^i{x',y',z') = 0. Substituting expressions for fp^i {x' ,y' , z') and 
fhemi (x' , v' , z' ) given in equations (A43) and (A45), we arrive at 

-1 



i=l 



^ijkk'ijro^k'^'i > - [r-] =0, i = < 



I 3 ^ ^ = 7 = 0'. 



(A47) 



It is convenient to rewrite [r -J^ in the form 
Equation (A47) can then be written as 



i=l 



( I ^ r' — r' , 



(A48) 



E 









n'p,i 



-1 



0, i = I 



' 1- 




= a = 






= r'- 


< 2- 




= /3 = 


v'\ 




= r' 


3 - 

v 










= r'. , 



(A49) 



Since the first two terms are already known, we can set each braced term equal to zero, 

Wif - eyfefcij-r^, fei/,' - [r^_,]^ =0, i = 1, 2, 3. 



(A50) 



The above relation, equation (A50), is valid in determining the set of discrete reflection points. The solutions of this quadratic 
equation are 

. 1/2 



+ Ki 



i = 1,2,3, 



where the summation over the indices j and k is implicit. The restriction of t-- being real imposes the condition 

+ KiY>0, i= 1,2,3. 



"1 




2 


n'p^i 



(A51) 



(A52) 



In spherical coordinates, the three radial vector components r'^ ^, and 3 are 

r'ii = r'^sinO' cos(j)' , r-2 = J'iSin^'sin^', r-3=r^cos^', (A53) 

where r ■ ^ = "i^'i x' ' ''^'i 2 = ^j' y' ^^'^ '^i 3 = z'- ^^^^ the terms r-, 9' and (f)' are the usual radial length, the polar and the 
azimuthal angle. This guarantees that R' = X]i=i ''i i^i is on the sphere, and justifies the step taken in equation (A50) since we 
are only interested in the conditions of the discriminants expressed by equation (A52). With ^ redefined in terms of spherical 
coordinates, the reality condition of in equation (A52) becomes 



[r'if sin^ e' cos^ </>' > 0, 



(A54) 



^mk'iyo,k 



+ [r'ifsm^ e'sm^ </.' > 0, 



(A55) 
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n'pA 



+ [r'f cos^ e' > 0. 



(A56) 



Equations (A54) through (A56) provide allowed range of v'^ to ensure its value being real. The solution of equation (A56) is 

2 



cos^ e' > - 



' 1 




-1 




n'p,i 





which leads to two inequalities, 



cos^' > i-—7 
- 2r 



(A57) 



(A58) 



These two inequalities cannot be satisfied simultaneously by the two vectors ^'o and k' i . We have to look for sin 9' by combining 
equations (A54) and (A55) to give 



which yields 



sin'' e' > — [r: 
4 ^ ' 



The solutions are again two inequalities, 



sin 61' > — 

- 2r 



n'p,i\\ {[^ijkk'ijro^k]^ + [e2mnk'i^mro,nT] + [r'ifsm'^0' > 0, 

[bijkk'ijr'o^k]'^ + [e2mnk'i^rnr'o,nY] , 



sin 61' < - 



2r', 



(A59) 



(A60) 



(A61) 



which cannot be simultaneously satisfied by the vectors R'o and k'l . We have to combine equations (A57) and (A59) to give 

tan^ 6»' > [e39rfci,,ro J | [eijkk[jro^^,f + [e2mnk[^^r'f^J^^ , (A62) 
which leads to another two inequalities, 

1/2 



UnO' > [esqrk'i^gr'o^^] [[^Ijkk'ijr'o^k] + b2mnk[^rnro,n\ I 

tan6»' < - [e3qrk[yo,r]~^ [bijkk'ijr'o^k]'^ + [e2mnk'i^rnr'o,nT] 
In the specified range for 6' , < 6' < n, the tangent function has the limits 



1/2 



lim < 6»' < -TT - lei 



< tan 61' < oo, 



lim (-n+ lei < 6' < n] ^ -oo < tan 61' < 0, 

e^o \2 ) 



(A63) 
(A64) 

(A65) 
(A66) 



where e is infinitesimal quantity introduced for limiting purposes. Since there is no guarantee that e3grfci,g?'o,r > in equations 
(A63) and (A64), one has to consider both cases where caqrk'i grQ ,. > and iaqrk'i gVQ ,. < 0. Therefore, for the positive 
denominator case where e3grfci,q''o,r > 0, we have 



e3grfci,,r[,_, > 0, lime^O {O < 0' < ^TT - \e\) , 

( -ifr 1^ "^''''^ 

< I tan6»' > [e3qrM,9''0,r] I [eyfefci,j''0,fej + [e2mnM,m''0,n] ' 



(A67) 



< oo; 
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l/2^ 



(A68) 



For the negative denominator case where e3grfci,q''o,r < 0) we rewrite equations (A63) and (A64) in the form 



tan 61' > - \e3qrk[,qr'o,r\ ^ {[eiifcfcij^o.fe]^ + [e2mnA;i,m^o,n] ^} 



1/2 



1/2 



(A69) 



(A70) 



The tangent function in the domain < ^' < tt has a discontinuity at ^' = 7r/2 , the inequality (A69) has the hmit > tan ^' > 
— oo, and the inequahty (A70) has the Umit oo > tan^' > 0. Therefore, the limits for a negative denominator case where 



< 



^3qrk[yo,r < 0, lime^O (O < ^' < ^TT - |£|) , 

tan 61' < \e3qrk[yo^^\~'^ | ey^fci^^r^ ^ + [e2mnfci_„r^)_„] ^| 
<^3qrk[^gr'o^^ < 0, lime^o (Itt + \e\ <6' <-k), 



1/2N 



(A71) 



< oo; 



OQ<[iaXie' >-\esqrki^qr'o^^\ { eijkk'ijr'Qk + [e2mnk'i^rnro,n\ \ <0 



1/2^ 



(A72) 



Comparing equations (A67), (A68), (A71) and (A72), we see that two of them are identical when rewritten in terms of the 
later convention where e3qrk[ qr'Q^^ < is expressed as — |e3grfci,g''o,r| — 0- The two tangent function inequality limits are 
summarized below for bookkeeping purposes: 



lime^o (0 < e' < Itt - \e\) , 



0< tan6''<|e3q^fci^,r^)_^| I ^ukk'^r'^ k + [e2mnA;;_„r^, J \ < oo 



1/2^ 



(A73) 



lime^o (57^+ <e' <Tr), 



00 < ^tan6'' > - |e3g^A;;_gr^,_^| ' | [eyfefc^ .,.r(, ^.j + [e2mnk[^rnr'o,nY \ )<0 
The corresponding arguments for inequahties in (A73) and (A74) are 

lime^o (0<^'<57r-|£|), 



1/2^ 



(A74) 



6»' = arctan |e3grfci,gr^,_^| \ <^ijkk[yo,, + [e2mnfc;,„r^, „] 



T 2 



1/2^ 



(A75) 



lime^o (sTT + |£| <9' <7r), 
6' = arctan ^- \e3qrk[^/o^r\~^ | eukkiy^ ,^ + [e2mnk[^rnr'o,nT^ 



1/2^ 



(A76) 
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The spherical coordinate representation is incomplete without the azimuthal angle (j)' . We have to solve for the allowed range for 
the azimuthal angle (j)' by combining equations (A54) and (A55). From equation (A54), we have 





"1 




-1 ^ ^ " 


> - 


2 




eijfefcij''o,fe 



and from equation (A55), 



They are combined to give 



sin^ <j)' > 



"1 




-1 


2 







[r^sin^^'] 



[rl sin2 0'] 



J 2 ii ^ r I / I ^ r i / / n 



2 A' 



COS 



The two inequalities are derived from the last equation, 

tan^' > [eyfefcijro fc] ^ e2mnk'i,m.'^o,n' tan0' < - [eijfcfci^j-ro^fc] ^ e2mnk[,mro,n■ 
ln the range of (/)', < 0' < 27r, the tangent function has the limits 



lim (^0<(l)' <^n- \e\j [0 < tan^' < oo] , 

» [-0O < tan^' < 0] 



lim ( -TT 
£^0 V2 



lei < (A' < TT- 



lim I TT + \e\ < 6' < -IT — \e 

e^o \ ' ' - ^ - 2 ' 



lim -.+ H< 



^> [0 < tan^' < oo] , 
' <2n- |e|^ ^ [-0O < tancj)' < 0] 



(A77) 



(A78) 



(A79) 



(A80) 



(A81) 



(A82) 



where e is an infinitesimal number used in the limiting process. Because discontinuities occur at (p' = 7r/2 and (j)' = 37r/2, the 
inequalities (A78) has the hmits 



< (tan(f>' > [eijkk'ijr'o^k] ^ e2mnfci,„ro^„) < oo, 



and 



-00 < (tancj)' < - [eijkk'i^jr'Q^^] ^ e2mnk'i^rnro,n) < 0- 
The ranges for inequalities in (A79) through (A82) can now be expressed explicitly as 

lime^O (0 < < ItT - |£|) , lime^o (tT + |£| < 0' < §77 - |£|) , 



< tand)' > 



(A83) 



lime^o (^TT + \s\<(f)'<Tr-s), lim^^o (|7r + |£| < 0' < 27r - \s\) , 



—00 < I tand)' < 



^Ijkk'ijrQ i. ^2mnfci,TO''0,ra ) < 0" 



(A84) 
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The solutions for 6' are 



lime^o (0 < < Itt - \e\) , \im,^o (tt + |£| < 0' < §77 - |£|) , 



t>' = arctan 



(A85) 



lime^o (sTT + |e| <(/)'< TT - e) , lim^^o (Itt + \e\ <(j>' <2it- \e\) , 



arctan 



(A86) 



In order to have f • values being real, the allowed range of 9' is determined by equations (A75) and (A76) and the allowed range 
of 4>' is determined by equations (A85) and (A86). Having found valid ranges of 9' and (p' in which v'^ is real, the task is now 
shifted in locating reflection points on the inner hemisphere surface in spherical coordinates. To distinguish one reflection point 
from the other, the notation i/- is modified to i^l ^ i^'i i in equation (A51). The first index 1 of ^ denotes the first reflection 
point. In this notation, the second reflection point would be 1/2 i and the A'^th reflection point, v'j^ ^. Then equation (A53) is used 
to rewrite r- ^ in terms of spherical coordinates. The Cartesian coordinate variables x' , y' and z' in equation (A51) are expressed 
as 



'1,1 



= x-,= 



j^O,k 



± 



-ijk 



1/2 



k'l'sin- 



e[ cos^ 



(A87) 



\2 = yi = 



' P>1 



<^2jkK,jr'o^k ± 



£2jkk'i,j''^0,k 



1/2 



+ [r'fsm^9[sm'' 



(A88) 



''1,3 = Zl = 



n'pA 



e3jkK,jr'o,k ± 



n' 1 



1 2 



^3jkK,j''''o,k 



1/2 



(A89) 



Although the first reflection point on hemisphere is fully described by R'l in equation (A12), it is not convenient to use R'l in 
its current form. The most effective representation of -R'l is in spherical coordinates. We set 9' = 9[ and (j)' = that describe 
the same reflection point R'l on the hemisphere. The subscript on the angular variables 9[ and (p'l denotes first reflection point 
on the hemisphere surface. In terms of Cartesian variables x', y' and z' , the first reflection point on hemisphere is given by 



1 ^'1,1 = x'l, 
R'l {x[,y[, z[) = ^ u[iei, z = <| 2 ^ u'^^ = y'l, 



(A90) 



i=l 



3 K.3 



The same point on the hemisphere, defined by equation (A90), can be expressed in terms of a parametric representation of 
equation (A12), 



3 r 



R'l (a,j-;-R'o,fc'i) =^ 



k'l 



k'lA 



— ^ ] ^l,i^i; 



where 



"^i.i = ^o,i + k'l fci_j, i = 1,2,3. 



(A91) 



(A92) 



Both representations, ^'1 {x'i,y'i, z'^) and -R'l ^^i,p; R'o, k'l^ , describe the same point on the hemisphere. Therefore, we have 

3 

^'1 {x'i,y'i,z'i) = R'l ^'0, k'l) ^ ^ Wi,i - Ti,^] = 0. 
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The components of the last equation are 

<i-TfM = 0> i = 1,2,3. 
Substituting expression of ^ from equations (A87), (A88) and (A89) into the above equation, we obtain 



"1 




2 


n'p,i 


"1 




2 





1/2 



+ [r'f sin2 9[ cos2 ct>[ } - Ti,i = 0, 



1/2 



Ti,2 = 0, 



(A93) 



(A94) 



(A95) 



p,i 



^3ifefcl,j^0,fe 



2 ^ 1/2 

,/l2 „2 /)/ 



[r^lcos^e'i^ -Ti,3 = 0, 



(A96) 



where Ti^j is defined in equation (A92). To solve for 6[, equation (A96) is first rearranged, 

^ 1/2 



± 



+ [r'if cos^ 0[ \ = Ti,3 - § n'p,i eajkk'ijr'o^k 



Square both sides and solve for cos-^ 9'i, the result is 



cos' e[ = [rr' 



n'p,i 



^3jkk'i,i''''o,k 



(A97) 



For reasons discussed earlier, 6[ information from the sine function is also needed. Following the earlier procedures, equations 
(A94) and (A95) are combined to yield the relation. 



1 

+ 2 



n'p,i 



£ljkKj''^0,k ^ 



e2mrafcl,TO^0,n ^ 



The equation is not easy to solve for sin^ . Fortunately, there is another way to extract the sine function which requires the 
knowledge of . The solution of 6[ is postponed until a solution of ^'^ is found. To solve for , it is desirable to solve for 
cos^ and sin^ from equations (A94) and (A95) first. Rearranging equation (A94), 



"1 




2 


n'p,i 



2 ^ 1/2 

+ [r'f Sin' 0[ cos' cl^'A - Ti,i 



n'p,i 



£2mrafcl,TO^0,r; 



1/2 



+ [r^]'sin2 0[ sin' <^'i 



Ti,2 = 0, 



± 



1 1 

+ [r'f sin' e[ cos' cl^'A =^hi-^ 



n'p,i ^ijkkijVQi^ 
and followed by squaring both sides, then cos^ <p[ can be found to be 

-2 r ^ ~i 

cos' (j)[ = [ri sin 0[] ^ - Ti,i n'p,i eijkk'ijVQ ,. 

Similarly, rearranging equation (A95), 



n'p,i ^ijkk'ijrQ,., 



(A98) 



± ■ 



n'p,i <^2jkk[jr'o^k 



1/2 



+ [r^] sin^ e[sin'(t)[ 



= Ti,2-^ 



e2jkk[jro^k 



and squaring both sides, then sin'^ can be found to be 



sin 



2 J.I 



r,- sinf 



T?,2 - Tl,2 



n'p,i 



e2jkk[jr'o^k 



(A99) 
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The function tan^ (j)'-^ can be obtained by combining equations (A99) and (A98), 

-1 



tan 



Tf 1 - Ti, 



j' 0,k 



-1 






-1 




1,2 ~ J- 1,2 


n'p,l 


e2mn^l,m''o,n 



(AlOO) 



Finally, the azimuthal angle 4>'i is found to be 



arctan 



± 



V 



T?,2-Ti,: 



l/2^ 



(AlOl) 



The restriction of ^'^ being real imposes the condition 



- Ti,i 



n'p,i 



^ijkk'ijrQf. 



or 



T2 'Y' 
1,2 ~ ^ 1,2 



eij/cfcl,j''o,fe; 



where <;>0. Following equations (A77) through (A86), the following results are obtained: 

lime^o (0 < (Ai < - |£|) , lime^o (tt + |e| < < Itt - |e|) 

l/2\ 



= arctan 



(A 102) 



lime^o (^TT + kl < 0'i < TT — e) , limg^o (f''' + kl < <f>i < 27r — 

ll/2\ 



: arctan 



"'~i_2-T"i,2||»'p,l|| e2m,.fc'i,^ro,„ 
Tff,l-'^i.i||"'p-i|r^^ij'=*='i,3''o,fc 



(A103) 



Having found the solution for cp'^ , we can proceed to finalize the task of solving for the polar angle 9'i . Combining the results for 
cos^ (^i and sin^ found in equations (A98) and (A99), sin^ 9[ can be found to be 



2 a' 



n'pA 



{'^lA^ljkKjr'o^k + '^l,2£2mnk'i,mr'o,n} 



(A104) 



The function tan"* 6[ is constructed with equations (A97) and (A104), 



Then 6[ is given by 

6'i = arctan 



"^"1,1 + "^1,2 - n'p^i ['^hieijkk[yo^k + ^l-2<^2mnk'i^mr'o,n] 



'^1,3 -'^1,3 n'p,i esgrk^^o^^ 



(A105) 



± 



V 



T?,l + T?,2- n'pA [^lAeijkkljr'Q^k + '^h2e2mnk'i^rnr'o,n] 



'I'l,3-'^1,3 «'p,l ^3qrk[^qro, 



l/2\ 



(A106) 



Following equations (A62) through (A76), we arrive at 

lim^^o (0 < e[ < Itt - |£|) 



9[ = arctan 



^l,l+'^l,2-|K'p.i|| ^{^i.f^ijkk'i Jr'g^+rl_2e2„^,^k[ ,^r'o^^} 



1/2N 



(A107) 
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O'l = arctan 



lime^o (571- + kl < 0[ < w) , 



l/2^ 



(A108) 



The allowed angular values are all defined now: by equations (A102) and (A103); and 9[ by equations (A107) and (A108). 
The initial reflection point on the inner hemisphere surface can be calculated by the equation: 



i=l 



1 — > i/^ J = r ■ sin^i cos t/i^, 

2 ^ ^[2 = r'^ sin 9[ sin ^[ , 

3 ^ K,3 = r'iCOS0[. 



(A109) 



We still have to determine the maximum wavelength that can fit the hemispherical cavity. It is determined from the distance 
between two immediate reflection points once they are found. We have to find expression describing the second reflection point 
R'2- In Figure 6, the angle ^1,2 satisfies the relation 



V'1,2 + 6'2 + 6»r- = TT. 

Angles 62 and 6r are equal due to the law of reflection, consequently 

02 = 0^ = 0^, 7/.i,2 =7r-2^i. 



(Alio) 



(Alll) 



It is important not to confuse the angle 9i above with that of spherical polar angle 6i which was previously denoted with an index 
i to indicate particular reflection point R'i. The 6i in equation (Alll) is an angle of incidence, not a polar angle. In order to 
avoid any further confusion in notation, equation (Alll) is restated with modifications appUed to the indexing convention for 
angle of incidence, 



i>i,i+l = TT - 29inc- 



(A112) 



The relation that cormects angle of incidence to known quantities k'l and R'l is 



where 



R'l 



k'l ■ R'l ^Y^k[ii'[i = r'i k'l cos9inc, 

i=l 

r'i, and the index i is not summed over. The incident angle 9inc is given by 

k'l 



9v 



arccos 



i=l / 



Substituting the expUcit expression of z/J ^ : 

u'i i= x'i= r'iSm.9'iCos(l)'i, u'i 2 = y'l = T'iSm.9'iSva.(j)'i, i/( 3 = = r- cos^i, 
into equation (Al 13), the incident angle is evaluated as 



9inc = arccos 



smf 



k', cos 6'-, + A;',, sin 6'-. + k', cos 9'^ 



V 



k', 

Vl 



k', 



(A113) 



(A114) 



(A115) 



where k'l i = k'^,^ , k'12 = fc^/ and fc^ 3 = A;^, . The second reflection point R'2 has the form 

3 

^'2 (^-2,1' ^^2,2, ^'2,3) = ^2,z 1^1,2, 1^1,3) ^i, i = < 



1 > 1/2/1 — -^2' ^'1,1 ~ -^l' 



^2,2 ~ 2/2' ^1,2 ~ Vl' 



(A116) 
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The relation that connects two vectors R'l and R'2 is 

R\-R'2 = [r'fcos^i,2, 



(Am) 



where 



R' 



(A112),as 



r ■ for a rigid hemisphere. Equivalently, this expression can be evaluated using ^pi^2, given in equation 

3 

[r'f cos (tt - 2einc) - J2 = 0- (A118) 



Equation (Al 18) serves as one of the two needed relations. The other relation can be found from the cross product of R'l and 

R'2, 



R'l X R'2 = '^^ijk'^'l,y2,k^i- 



(A119) 



1=1 



Since R'l and R'2 span the plane of incidence whose unit normal is given by equation (A17), 

i=l 

the cross product of R'l and R'2 can be equivalently expressed as 

11-1 ^ 

R'l X R'2 = -ri,2 n'p,i| ^eijkk'ijr'o^k^i, 

i=l 

where ri^2 is a proportionahty factor. The factor ri^2 can be found simply by noticing 



(A120) 



R\ X R'2 = ri,2n'p,i 



R'l X R'2 



:ri,2 



= ri,2, 



which leads to 



1,2 



R'l X i?'2 



[r-]^sin (tt - 20inc) ■ 



Equations (A119) and (A120) are combined as 



E 

i=l 



^ijk^l,j^'2,k + ^1,2 





-1 


n'p^i 


^ijkkijrQj. 



ei = 0. 



The individual component equation is given by 

^ijkK,j'^2,k + ^1,2 



n'pA 



^ijkki^jTQf. — 0, i — 1,2,3. 



(A121) 



(A122) 



(A123) 



Equations (A 123) and (Al 18) together provide the needed relations to specify the second reflection point R'2 in terms of the 
known quantities, -R'o, k'l and -R'l. It is convenient to expand equations (A118) and (A123) as 



-[dri,2/deinc]/2 
- — ^ 

[r^] cos (tt - 26'i„c) "K, 1*^2,1 - ^1,2*^2,2 - ^,3^^2,3 = 0. 

-1 

^1,2^^2,3 - ^^1,3^^2,2 + ^1,2 "'p,l eijfc^'l.j^'o.fe = 0' 

-1 

'^i,3'^2,i-'^M''2,3 + ri,2 n'p,i ^2jkk'ijr'ok = 0, 



(A124) 
(A125) 

(A126) 



n'p,i 



Equations (A124) and (A125) are added to yield 

^1,1^2,1 + Wl,2 + ^2,2 + [2^1,3 - ^1,2] ^"2,3 = ri,2 

Equations (A 124) and (A 126) are added lo give 

[^^1,1 - 2^1,3] ^2,1 + ^'1,2^2,2 + [2^1,3 + 2^1,1] 2^2,3 = ri,2 

Similarly, equations (A124) and (A127) are combined to give 

+ ^1,2] ^2,1 + Wl,2 - ^'2,2 + 2^1,3^2,3 = ri,2 

Define the quantities 

Q^l = K,2 + K,3> "2 = K,3 - 2^1,2. Cl = ri,2 
"3 = - a4 = I^'i,3 + '^1,1' C2 = ri,2 



n'p,i 



n'p,i 



1 dr.: 



2 d6in 



1 dr.: 



2d6i, 



1 dr.: 







n'p 


1 






n'p 


1 






n'p 


1 



2dei, 



^2jkk'ijr'Qi. 



1 rfFi 



2 dOinc ' 

1 <iri,2 



where Fi 2 is defined in equation (A121). Equations (A128), (A129) and (A130) form a reduced set 

^,1*^2,1 + «li^2,2 + 0:21^2,3 = Cl, a3t'2,l + l'l,2^'2,2 + ^41^2,3 = C2, a5'^2,l + OCqv'^^2 + 2^1,3^^2,3 

In matrix form it reads 

a5 Q'6 iy'i3 





v'2.1 




"Cl" 








C2 




v'2,3 




C3 



Mo 



and its determinant is expressed as 

det (Mo) = {v^^ + v{^2 + y'rA { [2^1,2]' + [K,3]'} = [r-^]' [^^m + ^^1,2 + ^^m] • 
Three new matrices are then defined here as 



Ml = 



Cl Oil 0L2 
C2 K,2 "4 
C3 "6 K,3 



Mo 



i^M Cl "2 
"3 C2 "4 
"5 C3 K,3 



M, 



2^1,1 "1 Cl 
"3 K,2 C2 

"5 ae C3 



^2,3 



det {m^ I det (Mo) 



The variables 1/21, ^2,2 2^2,3 solved with the Cramer's Rule as 

v'^ -^ = det (Mi) / det {m^ , v'^^^. = det (M2) / det (Mq) 
Explicitly, they are given by 

^2,1 = "2,1 = H,1 - 2^1,2 + "'1,3] Cl + H,2 - "'1,3] - ["1,2^- ["'1,3]^} C2 

+ {2^1,1 ^i,2 + 2^1,3] + ^2]'+ [2^1,3] '}C3) Wl,l + l^'l,2 + ^Ur^ [r'i\~^ , 
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^'2,2 = ^2,2 = {{K,2 + ^,3] + H,lT + Whs]^} Cl + {<,2 H,l ' ^,3] + ^,2]^} C2 



(AOS) 



where 



^ ^2,3 = ({ - ^'1,2] <,s - Ki]'- [^,2]'} Cl + { H,i + <,2] ^1,3 + + K2T} C2 

+ { Wl,2 - i^l,3 + [^1,3] ^} Cs) + Z^l,2 + I'ls] , 



K,i + K,2 + K,3 + 0, 



^'\,2 = y'l {i'i,Oi,(t)i) 

2^1,3 = z'l ir'i,Oi,(f)i) 



^2,1 (K,1'K.2=K.3) = 4 
1^2,2 (K,1,K,2>K,3) =y'2 
^'2,3 ^1,2> ^1,3) = ^^2 



(A136) 



(A137) 



In the above set of equations, £'2^^ has been used to indicate that 1^2,1 is now expressed expUcitly in terms of the Cartesian 

coordinates (i^J ii '^i 2j ^1 3) instead of spherical coordinates corresponding to the second reflection point, (r^, O2, 4>2) ■ The 
second reflection point inside the hemisphere is then from equation (Al 16), 

3 

-R'2 (i>2,l>^'2,2>^'2,3) = ^1^2,^, 

i=l 

where i>2 j, i = 1, 2, 3 are given in equations (A134) through (A136) with restriction given in equation (A137). In general, aU 
subsequent reflection points R'n can be expressed in generic form 

3 

R'n {i^'N,l^^'N,2^^'N,3) = ^N.i^i 



through iterative applications of the result i>2 , , i = 1,2.3. This however proves to be very inefficient technique. A better way 
is to express R'n in terms of spherical coordinates. Because R'2 belongs to a spanning set for the plane of incidence whose 
unit normal is n'p^i defined in equation (A17), the component relations i>2 j of equations (A134), (AOS) and (A136) satisfy the 
intercept relation given in equation (ASl), 



^2,i 



^ijkkijrQ j. ± 



1/2 



i = 1,2,3, 



where ^ = r- sin 62 cos ?!>2, ''i 2 = ''i sin 62 sin ^2 and r - 3 = cos 62- Here the subscript 2 of angular variables d'2 and ^2 
denote the second reflection point. In terms of the angular variables, using the above expression for i>2 i, the i/l^i, v'22 and £'2 3 
are expressed as 



T ■ 



2 \ 1/2 

+ H]'sin2^^cos2</.^l 



^likk't^o^k - ^2,1. (A138) 



^2jkk[yo,k 



1/2 



+ K]'sin2^^sin2.^'2 



n'pA 



^2jkk'ijr'o^k - ^2,2. (A139) 



n'p,i <^3jkk[yo,k 



1/2 



+ [r'l cos' 



2 a' 



n'p,! ^3jkkijrQi. Z>23. 



(A140) 



Square both sides of equation (AOS), cos (p2 can be solved as 



cos^ (b'o = Ir'i sin ( 



n'p^i 



^ijkk'ijr^^^ ^2,1 



J^'p.i ^ijkk'ijr'fju 



(A141) 
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Similarly, square both sides of equation (A139), sin^ 02 be solved as 



sm 



2 A' 



r,- smf 



1 ^ 


1 






2 





^2jkk'i,j''''o,k 



(A142) 



The last two equations are divided to give 



tan^ (b'o = 







-1 


2 






-1 


2 


1 

2 


n'p,i 


^2mnki^rn'''0,n ~ ^2,2 




1 
2 




£2mn^l,TO''0,n 





















^Ijkk'ijVQi. 1^2,1 



The azimuthal angle 02 is obtained as 



/ 



i)2 = arctan 



± < 



V 



1 

2 


n'p,i 


-1 

e2mnfcl,TO^0,n ~ ^2,2 


2 


1 
2 

















l/2\ 





1 






2 


n'p,i 









^Ijfc^'lJ^O.fe ^^2,1 



^ijkk'ijrQi^ 



Following the procedures used in equations (A77) through (A86), the following results are arrived at 

lim^^o (0 < 0^ < - |£|) , lim^^o (tt + |e| < (^2 < ^tt ~ \e\) , 



02 = arctan 



■||n'p,l|l e2™,ifel,„i-^, „-i>^ 2J -[||ln'p,l|l e2™,ifel,„r-^, „J 
[5|K'p,i|r^eijfcfel,j'-^,fc-'>2,i] -[sIK'p.iir^eijfcft'i ,,.r^_^] 



1/2N 



(A143) 



lime^o (^TT + \s\ < (f)'2 < IT - e) , linie^o (Itt + |e| < 0^ < 27r - 



= arctan 



n'p.l £2^ 



2||"'p,l|| ^l3'«^l,j''o,fc ^2,1] 



l/2> 



2||"''p,l|| ^lj'='^l,j'"o,fc] 



(A144) 



The solution for 02 forms a generic structure for any subsequent reflection points on the inner hemisphere surface. The Nth 
azimuthal angle 0^ is found following a prescribed sequential steps 



K ^ 02 ^ 03 



(A145) 



By reversing the direction of sequence, 0^ can be expressed in terms of the initial azimuthal angle 0'i, 0^ = 0^ (0'i) . The 
polar angle 62 of the second reflection point can be found by squaring equation (A140), which yields 



n'p,i e3jkk[ jr'o k - i^'2,3 



1 ^ 


1 






2 


n'p,i 



(A146) 



Add together equations (A141) and (A142), sin 62 can be solved as 





2 


1 




'>2,1 




2 


n'p,i 





'1 




-1 


+ 


2 


n'p,i 





e2mra'C]^ ,„rg „ — l'2^2 



1 ^ 


1 






2 


n'p,i 



n 2' 



(A147) 
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By dividing equations (A146) and (A147), we get 

-1 



h- 1 




2 







'1 




-1 




2 









2 


"1 








2 






2 


1 








2 




-^2,3 


2 


"1 






2 





n 2 



The polar angle 02 is given by 





'1 




-1 


arctan ^± 


2 







1 ^ 


'1 






2 


n'p,i 



^2mnki^jnrQ^ri 



+ 



f^'p,l ^2TOnfcl,TO''0,n ^2,2 





'I 






2 


fl'p,! 



1 ^ 


"1 








2 


n'p 


1 


1 ^ 


1 








2 


n'p 


1 



n 2 



1/2 



-1/2N 



eSgr^l.q'^O.r ^^2,3 

Following the procedures given in equations (A62) through (A76), the following results are obtained 



6*2 = arctan 



n'p,i 



0'2< 




-N) 






2 r 


0,k ^2,1 












1 






2 


n'p,i 



1 ^ 








1 






2 


n'p,l 



^9 = arctan 



+ 



^2mnki m''Q,n ~ ^2, 
'^'p,l ^3grfcl,q''o,r ~ ^2,3 

lime^o (^TT + |e| < 61^ < tt) 
^ -1 

'^'p,i ^ijkkijrQ^. — 1^2,1 

^2mnki^rn''"0,n ~ ^2,2 
n'p,! ^3qrk'iqr'Qj. — i>2^^ 



>> 1/2' 



^2mnki raVQ . 



-1/2N 



1 ^ 








1 

2 


n'p,i 



n'p,i 



-1 2 






-1 




1 
2 




£2mn^l,TO^0,n 



1/2 



1 ^ 






-1 




1 

2 


n'p,i 


^SqrkigTQ,^ 



-l/2\ 



(A148) 



(A149) 



The above result of 62 forms a generic structure for any subsequent reflection points on the inner hemisphere surface. The N\h 
polar angle 6'^ can be obtained by foUowing the sequential steps 



9[ -> 6*2 



(A 150) 



Equivalently, reversing the direction of sequence, O'j^, can be expressed as a function of the initial polar angle 6[, = 9'^^ (9i) . 
With angular variable ^3 defined in equations (A143) and (A144), and 6*3 defined in equations (A148) and (A149), the second 
reflection point on the inner hemisphere surface is given by 



(A151) 



3 




' 1 - 


^ ^2,1 


= r[ sin 61^ cos(?i^, 


2{r'i,e'2A'2)=Y.'''^,i 
i=l 


2) 6i, i — < 


2- 


^ ^'2,2 


= r'i sin 9'2 sm(j)'2, 






, 3- 


"2,3 


= r[ cos 9'2. 
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¥i+l,i+2 



Figure 26: The two immediate neighboring reflection points R'l and R'2 are connected through the angle V'l.a- Similarly, the two distant 
neighbor reflection points R'i and R'i+2 are connected through the angle ,V'i+i,i+2 • 



As shown in Figure 26, two reflection points R'l and R'2 are related through ^1,2, which is the angle measured between the 
two. Since R' j. where the index j = 1, 2, • • • , N„iax and Nmax is the last count of reflection before a repeat in cycle, belongs 
to a spanning set for a plane of incidence whose unit normal is n'p^\ given in equation (A17), all reflections occur on the same 
plane of incidence. The task of determining the A'^th subsequent reflection point R' ^ is therefore particularly simple. The needed 
connection formulae between the initial reflection point R'\ and the JVth subsequent reflection point i?V is found through both 
scalar and vector cross product relations similar to those given in equations (A117) and (A119). In order to generalize the 
previous result for R'2 to R' n-, recall the set from equations (Al 17) and (Al 19), 




<^ijkVl,3^'2,k + Tl,: 



1,2 - 




<,A,i = 0, 






-1 


ri,2 







(A152) 



ii = 0, 



where Fi 2 = [r ■] sin ipi^2 and V'1,2 = tt — 20inc- Because R'n, R'i and R'2 belong to a same spanning set forming the plane, 
it is true that 

^'1 X R'n oc R'i x R'2. 

Therefore, we can write 

^'1 X R'n = T[^nR'i X R'2 = r;_jvri,2n'p,i = Ti.jvn'p,!, 
where Fi^at = F'^ jyri^2 is a proportionality factor. Comparing the results, 

{^'1 X R'n = ri,jvnVi} ^ {^'1 x R'2 = ri,2n'p,i} > 
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one obtains a set of relations similar to equation (A152) for R'n, 

[r'f cos Vi,jv - Ei=i K,i'''N,i = 0' 

-1 



R'l ■ R' 



N 



(A153) 



Ci = 0. 



In the above expression the index N onu'j^ ^ and v'j^ j. denotes components corresponding to R'n and tpi,N is the angle measured 
between R'l and -R'jv- The proportionahty factor Fi^jv is found to be 



R'l X R' 



N 



AT 



N 



[r-]^sin^i,jv. 



R'l X iJ'jv = ri,jvn'p,i 

which yields 

Tl^AT — 

The angle V'l.w is contained in f^Vi,2,i/'N-i,iv ^s shown in Figure 26, 

^K'i.H'k = ^'/'i,2,'/']v-i,n = "01,2 + V'2,3 H h V'iV-2,JV-l + '0JV-l,Ar- 

For each tpi,i+i > the sum of inner angles of a triangle gives 

1pl,2 + O2 + Or = TT, 
tlj2,3 + 0s + Or+1 = TT, 



V'JV-2,JV-1 + + Or+N-2 = TT, 

'4>N-1,N + On + Or+N-1 = TI". 

The law of reflection gives 

O2 = O3 = ■ ■ ■ = On-1 = On = Or = Or+l = ■ ■ ■ = 6r+N-2 = ^r+JV-1 = ^inc- 

Hence, the angles are found to be 

V'1,2 = "4)2,3 = ■■■ = tl)N-2,N-l = 1pN-l,N = TT - 2^i„c- 

The angle fivi,2,V'N-i,N expressed as 

^^j/'i,2,Vn-i,n = i^\,2 + '4>2,3 H h tlJN-l,N = [tT - 2^j„c] + [tT - 2^j„c] H h [tT - 261,, 

or 

V'l.iV = fij/'i,2,V'N-i,N = [-/V - 1] [tT - 20inc] ■ 

Hence for Fi^jv, we have the result: 

Ti,N = [r'ifs\n{[N-l][K-20inc]). 
The angular variables 6'j^ and 0^ corresponding to Nih reflection point R' n are given as 

lim,_o {^<0'n<\tt- \e\) , N>2, 



^'n>2 — arctan 



+ 



n'p,i 



n'pA 







-1 








2 


1 

2 






0,k 












2 








e2mnfci,TO^0,n ^iV,2 




1 
2 





-,2^ 1/2 



1 ^ 








1 






2 


n'p,i 



^2mnki^rnTo,n 



^Sqrk'iqr'Qr 



-1/2N 



(A154) 



(A155) 



(A156) 
(A157) 



(A158) 
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lime^o (sTT + |e| < 6*^ < tt) , > 2, 



^'n>2 ~ arctan 





1 




-1 


(-{ 


2 







1 ^ 








1 






2 





+ 



1 ^ 








1 






2 





n 2 



Af,3 



1 ^ 








1 






2 





1/2 



-1/2N 



(A159) 



lim^^o (0 < <t)'N < - |e|) , lim^^o (tt + |£| < (^'n < - |e|) , N >2, 

1/2N 



''N>2 



arctan 



[|||"'p,i|| ^e2mnfe'i,„ro^„-£';v,2] -[l||"'p,i|| ^e2m,.fei_„r^_„ 



(A160) 



lim^^Q {Itt + \e\<(j)'j^ < IT - e) , linie^o (Itt + |£| < (/)^ < 27r - |e|) , iV > 2, 

12 >! l/2> 



^jv>2 = arctan 



r2~ 



(A161) 



where {I'^vi • ^ = 1)2,3} are given in equations (A134), (A135) and (A136) together with the modification j>2 i 
Ci (ri,2) ^ Ci (ri,iv) , C2 (ri,2) C2 (ri.iv) and Ca (ri,2) ^ C3 (ri,jv) , where ri,jv is given in equation (AIS?). With 
angular variable 6'^>2 defined in equations (A158) and (A159) and 4''n>2 defined in equations (A160) and (A161), the Nth 
reflection point on the inner hemisphere surface is given by 



R'n {r[, 0'^, ct>'j,) = ^'N,i ir'i, e'^, cj>'j,) 



I = < 



i=l 



' 1- 




— r[ sin O'pf cos 


i>'N, 


< 2- 




= r ■ sin sin ( 




3- 

v 




= r^cos6'^, 





(A162) 



where the initial reflection point R'l is given in equation (A 109). 

For a sphere, the maximum number of reflections are given by the equation 

Ns,max''pN-2,N-l = 277, 

where ?/)]v-2,]v-i is the angle between two neighboring reflection points R'n-i and R'n-2; the subscript Ng^^nax denotes the 
maximum number of reflection points for a sphere. The above result is a statement that the sum of aU angles is equal to 2n. 
Application of the rule shown in equation (A155) for ipN-2,N-i gives 



^max ^^inc\ — 27r > ^s^max ' 

where 6inc is given in equation (A115). In expUcitform Ns^max is given by 



2-K 



(A163) 



iV., 



2tt 



TT — 2 arccos 



sin 6\ k'l cos S'-i + k',, sin 6\ + k', cos 6\ 



V 



ki, 



k', 



(A 164) 



Appendix B: MAPPING BETWEEN SETS (r, 9, 4>) AND (r', 9', 4>') 



In this appendix, the original derivations and developments pertaining to the mapping between the sets (r, 6, (f)) and (r', 9' , (p') 
used in this paper are described in detail. 
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For a sphere, the natural choice for origin is the sphere center from which the spherical coordinates (r^, 0', cj)') are prescribed. 
For more complicated configurations, as shown in Figure 7, the preferred choice for the origin depends upon the problem in 
hand. For this reason, this section is devoted in deriving a set of transformations between (r^, 0' , </>') and (r^, 9, 4>) , where the 
primed set is defined relative to the sphere center, and the unprimed set is defined relative to the global configuration origin. In 
Cartesian coordinates, the two vectors R and R' describing an identical point on the hemisphere surface are expressed as 



(Bl) 



where 



Ul = X 
11^3 = Z 



^2 = y 



x'- 




" e 1 








y' 
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Here R and R' are the position vectors of the same location relative to the system origin and the hemisphere center, respectively. 
The two vectors are related through a translation, 



R (i/i, V2, V^) = Rt (z^T,1, 2^T,2, «^T,3) + R' {'^[,'^2, ^'s) = I'^^'' + '^'i^ 



(B2) 



where Rt = J2^=i ^T,i&i is the position of hemisphere center relative to the system origin. Equation (B2) can be written as 



Eh- 



i=l 



and the component equations are 



!^i = 0, 



0. 



1,2,3. 



(B3) 



(B4) 



It is to be emphasized that in the configuration shown in Figure 7, the hemisphere center is only shifted along y by the distance 
i/r.2 = a, therefore VT,i^2 = 0. Nevertheless, the derivation is done for the case where VT,i ^ 0, i = 1, 2, 3 for general purpose. 
In expUcit forms, they are written as 



J^i - vt,i —v[=Q, Vi — Vt,2 -1^2 = 0, vz — vt,z - ^3 = 0. 

In spherical coordinates, 

y\=ri sin 6 cos (f) = r, Ai {6, (j)) 
f2 = r, sin sin (/) = A2 {6, 4>) 

1^3 = n cos 9 = Ti A3 (9) 

equation (B5) is written as 



(B5) 



; =r^sin6i'cos.?i' = r^A'i(6l' -^'^ 
i =r'ism9'sm(l)' = r'iA'2 [9' 



(B6) 



Tj sin 9 cos 4> — vt,\ — r[ sin 9' cos (j)' 



ri sin 9sm<p— 1/^,2 — sin 9' sin </>' = 0, 



Vi cos 6 — UT,i — r'i cos 6' = 0, 



(B7) 
(B8) 
(B9) 



where the Cartesian variables {ui, v'^ : i = 1,2, 3} were expressed in terms of the spherical coordinates. The cos^ and sine 
functions are obtained from equations (B7) and (B8), 



COS0 = 



r; sm f cos < 



n sm ( 



sm0 = 



i'T,2 + sm a sm < 
n sin 9 
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The azimuthal angle is given by 

<p = <p 4'\ ^T,i, vt,2) = arctan 



i^T,2 + r ■ sin 9' sin 4>' 
t^T,i + t[ sin 6' cos ^' 



(BIO) 



where the notation (j) indicates that (j) is expUcitly expressed in terms of primed variables. Combining equations (B7) and (B8), 
we have 



which leads to 



Ti sin 6 [cos (j) + sin <f)] — vt,i — i'T,2 — sin 9' [cos </>' + sin (f)'] = 
t^T,i + t^T,2 + sin 9' [cos <p' + sin <p'] 



smu = 



Ti [cos (j) + sin (f)] 
From equation (B9), we have 

cos 9 = [z/T,3 + r'i cos . 
Combining the above results for sin 9 and cos 6; and, solving for the argument 6, 

' vt 1 + i^T 2 + r'j sin 9' [cos (j)' + sin ( 



9 = arctan 



cos (f) + sin 4>\ [vT,3 + Tj cos ( 



(Bll) 



where is to be substituted in from equation (BIO). For convenience, the above result for 9 is rewritten explicitly in terms of 
primed variables, 



9 (ri,9',(t)',RT^ = arctan 



{i^T.i + i^T,2 + sin 9' [cos (j>' + sin (f>'] } [vt.s + t[ cos 9'] ^ 
cos 1 arctan ( ''^'^"J"'',' '^'"f/'" ) ) + sin ( arctan ( "T.^+r^smO sm (j> \ \ 

V \ \ I'T.i+rJ sin 9' COS <p' y y ir.i+r^ sinS' cos <p' y y , 



(B12) 



Here the notation 9 indicates that 9 is explicitly expressed in terms of primed variables. The magnitude of a vector describing 
hemisphere relative to system origin is found from equation (B2), 



n 



1/2 



R 



i=l 



^ A'i(6»',^') =sin6''cos^!)', 
{9',(p')=sm9'sm(l)', 
L A'3{9') =cos9'. 



In terms of spherical coordinates, the position vector is expressed as 



1/2 



R 



(rl, A, A', Rt) = \Y, [yT,i + AK'^f \ ^^e- 



. j=i 



i=l 



h {9,4>) 



A, 



= sin 9 cos 4>, 



sm 9 sm (p, 



cost/. 



(B13) 



Appendix C: SELECTED CONFIGURATIONS 



In this appendix, the original derivations and developments in this paper pertaining to the selected configurations: (1) the 
hoUow spherical shell, (2) the hemisphere-hemisphere and (3) the plate-hemisphere, are described in detail. 
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1. Hollow Spherical Shell 



For the reflection problem in a sphere as shown in Figure 8, the natural choice for a system origin is that of the sphere center, 
R' = 0. The A'^th reflection point inside the sphere is given by equation (A162) as 



R 



i=l 



K.N. 2 = risine' ,^ sin 



where the label ,s have been attached to denote the sphere. Keeping in mind the obvious index changes, the angular variable 9'^ jy 
is defined in equations (A158) and (A159), and 0^ ^, in equations (A160) and (A161). Staying with the notation of equation 
(B13), R's N is rewritten as 



Rs,N (r-^, A;,;v) =riYl K,N,A, < 

i=l 



sin 9'^ jvcos (p'^ j^, 



K,N,3 K.iv) =COS^^_;v. 



(CI) 



where the relations yT,s,i = 0, and J2i=i [-^'s at i] ^ = 1 are used. 

The maximum number of internal reflections for a spherical cavity before a repeat in cycle is given by equation (A163), 



2n 



where 9inc is given in equation (Al 15). 



The distance 



between two immediate neighboring reflection points on a sphere is 



It should be noted that 



R 



s,2 



L 



Rs,2 [r'„A',^2)-Rs,i (r',,A's,i) 



R 



s,3 



R. 



J =3,--- ,7V«, 



(C2) 



(C3) 



The only reason that Rg.i and Rs.2 are used is for the purpose of convenience. 

To compute the resultant wave vector, k' inner , acting at the point , the incident wave is first decomposed into components 
paraUel and perpendicular to the local normal vector, — i?/^ i , of the inner surface 



k'i ' R' s.l 



R'sA + 



X R 



s,l, 



where the subscript (+) of denotes the particular contribution where the incident wave k'i is approaching R' s,i from R' sfi- 
From equation (A14) of Appendix A, the corresponding reflected wave vector can be expressed in terms of the incident wave as 



J- -|- — f 



R's,i X { [k'i ■ R' 
R's,i ■ { [k'i ■ R's, I 



s,l 



R'sA + 



X R 



R's,i + 



Ft' s,l ^ ^'i 



R's,l X 
R'sA 



[k'^ ■ R's 



R'sA+R's.i X 



R's.l X fc',| X R's,l]\ X R's,l 



{ [r's,i X fc*',] 



k'i ■ R's,l + R'sA ■ { [R'sA X k'^^ X R's,l}] R's,U 



where ar,± = a^^w = 1 and n' - 
expression is simpUfied to 

P , 



R's 1 • Because R's i -L I [R's,1 X fc'ij X R's, i| and R's,i \\ [k'i ■ -R's.ij -R's.i, the above 



R'sA X 



{[R's 



,1 X k'i 



X R'sA } X R'sA 



k'i ' R^s.i 



R'sA. 
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The changes in resultant wave vector k' inner at the point R's,i due to k'i^+ at location R's,o is given by 

^k'inner,+ -R's,l ) -R's,0^ — ^'r,+ ~ 



{ 



R'sA X k\ 



R's.i X k'. 



X R 



'-A 



X R'sA 



X i?'. 1 - 2 



k\ ■ R'sA 



R'sA 



k'i ' R' sA 



R's,l, 



where 



R's,i X { [r's,i X k'i] X R's,i}] X ^',,1 = 



For the incident wave traveling in the opposite direction, i.e., approaching R's,i from R's,2, one has 



k'i ■ 



-k'. 



-k 



k'^ 



-h' ■ — —h'- , 



where the subscript (— ) onk'i- denotes the particular contribution where the incident wave k'i is approaching irovtvR' s,2- 
In this case, the changes in the resultant wave vector k' inner at the point R' s,i due to k'i - at the location i?'s,o is given by 

^k inner,— ^5 R 5,1) R 5,0^ — k — k i^— — k j^-j- ~1~ k r,-\- — ^^k inner, -\- R 5,1) R 5,0^ * 

The resultant wave vector k' inner acting at the point R' s,i due to incident wave approaching R' s,i from R' s,o and the other 
incident wave approaching R' s,i from R' s,2 is therefore 

^k inner ^) R s,l ) R 5,0^ — ^k inner, -\- ^1 R s,l^ R s,0^ ~t~ ^k inner,— (^1 R s,li R s,0^ 



= -4 



• R'sA 



R's,i, 



where the subscript b of fc'^ 5 denotes the wave vector for ambient fields inside cavity. 



The wave number 



i,b 



that can be fit in the bounded space of a resonator is restricted by the boundary condition 

1 









-1 


k'i^b 


= mr 


L 


= mr 



The scalar product of k'i^b and R's,i is 



k'^,b ■ R'sA 



COS 9r, 



where the angle between the two vectors k'i^b and Rs,i is equal to the angle of incidence 6inc, as shown in equation (Al 15). The 
momentum transfer is proportional to 



^k' inner R' s,!-) R' s,Q^ 



Ann cos Oj. 



Rs,2 (r^A;J-i?,,i (r'„A',^,) 



R's,i, 




(C4) 



Similarly, the resultant wave vector k' outer acting at point R's,i + aR's,i on the outer spherical surface, where a is the sphere 
thickness parameter, is given by 



'outer (; -R's,l + a-R's,l) = -4 k ij ■ R! s,l R's,l, 



where the subscript / of k'ij denotes the wave vector of the ambient fields in free space, and the factor 4 is due to the fact there 



are two incidence wave vectors from opposite directions. The free space wave number 
due to the boundary. And, the scalar product of kij and R's,i is 



k'i. 



has no quantization restriction 



ki,f ■ R's,i = 
The momentum transfer is proportional to 



k 



cos (tT - 6inc) = - 



(lR's,l 



aR' 



5,1 



k'i. 



COS 0incR's,l, 



COS 6incR's,l- 

< Oinc < 7r/2, 
n=l,2,--- . 



(C5) 
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2. Hemisphere-Hemisphere 

For the hemisphere-hemisphere configuration, the preferred choice for a system origin is that of ^ = 0. The A'^th internal 
reflection point is given by equation (B13), 



1/2 



(C6) 



1=1 



where the label h denotes hemisphere; and 



A?i,JV,l (Oh,N,(i>h,N^ = sin 9 COS (j)h,N, 
^h,N,3 (^h,N^ = C0s6h,N- 

The definitions for A'^ j, i = 1, 2, 3 are identical in form. The angular variables (^dh,N, 4'h,N^ 
and (B12), where the obvious notational changes are understood. The implicit angular variable 6'^ ^ is defined in equations 
(A158) and (A159); and 0'^ jy, defined in equations (A160) and (A161). 

We have to determine the maximum number of internal reflections of the wave in the hemisphere cavity before its escape. 
Three vectors, R'o, ^ik'i and R'h,i = R'i, shown in Figure 26 of appendix A satisfy the relation 

R'h,i=l — R'o = Ci=lfc'i=lj 

where the notation h of R'h,i=i denotes the hemisphere. The path length squared is given by 



are given in equations (BIO) 







2 




2 




R'h,i=l — i 






R'h,i 


+ 


R'o 



2R'h,i ■ R'o = [r'il + 



R'i 



2r' 



R', 



cosV'0,1- 



Since 



R'h,i=i — R'o 



V'0,1 = arccos ( 2 ^ 



= , the angle V'0,1 is found from the last equation to be 

1 



R'o 



+ [r. 



R'i 



R', 



1 r ' 



where ^1 = ^i,p is given in equation (Al 1). The angle V'1,2 measured between the two vectors R'h,i and R'h,2 is 



V'1,2 = arccos (^[r-] ^ R'h,i ■ R'h,2^ = arccos (^Yl^'h,hi^'h,2,ij 



where R'h,i and R'h,2 have been exphcitly written for iV = 1, 2 in equation (CI), or equivalently, 

V'1,2 = TT - 29inc 

from equation (Al 12). For a hemisphere, it is convenient to define a quantity 



'^h,max = -, [tT - V'O,!] 

Vl,2 



or explicitly. 



1 



TT 2^^2nc 



TT — arccos 



R'o 



+ [r. 



R'o 



R'o 



>2 



(C7) 



where ^i_p is given in equation (All) and 0inc is given in equation (Al 15). The maximum number of internal reflections is then 
simply 



.max — \^n,max\ri ? 



(C8) 
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where the notation 



The distance 



^h.rnax]Q IS the greatest integer of Ih^max and it is defined in equation (C7). 
between the two immediate neighboring reflection points of the hemisphere is 



It should be noted that 



Rh,2 — Rh,l = Rh,j — Rh,j-1 , j = 5,- ■ ■ ■,^h,r, 



(C9) 



(CIO) 



and the only reason that Rh,i and Rh,2 are used is for the purpose of convenience. 

The change in wave vector direction upon reflection at the point R'h,i inside the resonator, or at the location R'h,i + o,R'h,i 
outside of the hemisphere, is given by results found for the sphere case, equations (C4) and (C5), with obvious subscript changes. 



^k'inner R' h,l, R' h,0^ = 



4mr COS ( 



R'k 



< Oinc < 7r/2, 



(CU) 



n = l,2,- 



and 



outer 



; R'h,i + aR 



COS O^ncR'h,!, 



< 6»i„c < 7r/2, 
n = l,2,-- - . 



(C12) 



The above results on Ak' inner {] R'h,i, R'h,o^ and Ak' outer R'h,i + (iR'h,ij have been derived based upon the fact that 
there are multiple internal reflections. For a sphere, the multiple internal reflections are inherent. However, for a hemisphere, 
it is not necessarily true that all incoming waves would result in multiple internal reflections. The criteria for multiple internal 
reflections are to be established. For a given initial incoming wave vector k'l, there can be multiple or single internal reflections 
depending upon the location of point of entry into the cavity, R'o- Shown in Figure 9 of section (3.B.2) are two such reflections 
where the dashed vectors represent the single reflection case and the non-dashed vectors represent the multiple reflection case. 
Because all the processes occur in the same plane of incidence, the relationship R'f = —XqR'q with Aq > has to be true. 
Therefore, we will have 



= R'o + ipk'i, R'f = -XoR'o = R'l + k'2. 
After eliminating R'l from the last two equations, we find 

^'o = -[l + Ao]"' kpk'i + k'2 



(C13) 



(C14) 



The direction of the reflected wave vector k'2 cannot be arbitrary because it has to obey the reflection law. The relationship 
between the directions of an incident and the associated reflection wave is shown in equation (A 14). Designating n' = — i?'i/r 
k'r — > k'2 and k'i — > fc'i , the reflected wave vector k'2 can be written in the form 



k'2 oc ari [r- 



R'l X k'l 



X R'l - ani [r-] R'l ■ k'lR'i. 



By introducing a proportionaUty factor A2, it becomes 

-2 



k'2 = X20ir,± r-j 



R'l X k'l 



2 ^, 



xii'i-A2Q,i| [r,;] '■R'l-k'iR' 



The goal is to relate i?'/, or Aq, in terms of R'q. Substituting the expression for R'l from equation (C13), we arrive at 



k'2 = -^pA2Q!^,|| [r-] ^ k'l + CpA2 [r-] ^ ja^,. 



+ A2 \r: 



R'o X k'l 



X R'o — Ur- II R 



R'o X k'l 
'o-k'iR'o} 



X k'l — ttr I 



R'o ■ k'lk'i + 



k'l 



R 



(C15) 
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Finally, equations (C14) and (C15) are combined to yield 



R'o X fc'] 



X fc'i — a,.. 



+ar\\R'(i ■ k'lR'o - art 



R'o X k'l IxR'o- A2 ' [r'f [1 + Ao] R'o = 0. 



Utilizing the formula 



Ax B 



c = EU{ 



A-C 



Bi- 



B-C 



'j A( I e;, the cross products are evaluated as 



R'o X k'l 



X k'l 



k'l 



R'o X k'l 



R'o X k'l 



i?'o • fc'i 



^'0 = Ell 



R'o 



k'l ■ R'o 



^0,1 r ei, 



k'l 



k'l = Ef=i k'l.ei, R'o = E?=i r'oA. 



3 „/ 



We can rewrite equation (CI 6) as 



1=1 



R'o ■ k'l 

k'l 



k'l 



k'i,i + 



k'l 



'0,1 



R'o ■ k'l 



k'l 



ki,i 



'0,1 



R', 



X^'ir'f k'l 
' k[, - k'l ■ R'or'oj 



k'i,i I + Q-r,\\R'o ■ k'lr'o^i - ar,± 



^2'[r'if[i + Mr'oAei=0, 



which leads to the component equations 





















("'■,11 


R'o ■ k'l 


k'l 


'k'i^i + 


k'l 


r'o,i 




R'o ■ k'l 


k'l 



k'lj 



k'l 



'0,1 



^2 ^ [r'if k'l k'l^i ) + ar^\\R'o ■ k'lr'o^i - ar,± 



R'i 



' k'l^i - k'l ■ R'or'o^i 



\2^Kf[l + \o]r'o,=Q, 



where i = 1, 2, 3. For an isotropic system, ar^i. = ar,|| = ^r, the last equation reduces to 

-1 



^l^rk'i^l + 



2ar 



k'l 



k'l 



k'li 



+2a,i?'o • k'lr'o^i - a. R'o k',j - X^^ [r'f [1 + Aq] r'o, = 0, 



(C16) 



where Z = 1, 2, 3. Because there are three such relations of the above, all three component equations are added to yield 

3 r 

'2\\k'i\\r'o,i-^2'[r'f 

=1 L 



=1 



E 



k'l 



k'l, 



.1=1 



^ |2i?'o • k'lr'o^i - \\R'of k'l^i - X^' [r^]' [1 + Aq] r'oA = 0, 
1=1 ^ 



(C17) 



where the both sides of the above equation have been multiplied by 



Ei=i k'l I 



Since is just a positive root of ^1 of equation (A9), it satisfies the equation 

1 



and ar = I have been chosen for simplicity. 



1=1 



1=1 



1=1 
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where the index i in equation (A9) have been changed to I. The reflection coefficient have been set to a unity in equation 
(C17) for the very reason that = 1 had already been imposed in equation (A9). Because = ^p, and the fact that coefficient 
of ?p = Ci — 1 in equations (A9) and (C17), the two polynomials must be identical. Therefore, subtracting equation (A9) from 
equation (C17), we obtain 



L;=i 



E 

1=1 



2 fc'i 



'o,( 



-1 r„'i2 



k'l 



^ |2#o • k'Wo^i - ||i?'of fci,; - A^-i [r'f [1 + Ao] r'^^i - [r'^^if] = 0. 
1=1 ^ J 



Because is a particular value for the root of ^i, for the case where ^ 0, the above equation is satisfied only when the 
coefficients of the different powers of vanish independently. This is another way of stating that each coefficients of the 
different powers of in equations (A9) and (CI 7) must be proportional to each other. For the situation here, they must be 
identical due to the fact that coefficients of = = 1. Hence, we have the conditions: 



E3 y 
1=1 H,i 1^1=1 



k'l 



S;=i '''0,1^1,1 ~ 0' 



(C18) 



From the first expression of equation (CI 8), we find 



k'i,i-K'[r'if[^ + Mr'^,i 



'0,1 



0. 



a: 



2 k'' 



H -1 



=1 



E 



\k'i 



^1,1 E 



'0,1- 



(C19) 



Solving for Aq from the second expression of equation (CIS), we have 

2 



Ao = 



U=i\k',^i+2R'o-k'ir',^i 



R'o 



m=l 



- A" V 



0,( 



•^2 J2l=l ''"'q.I 

or, substituting the expression of A2 ^ given in equation (C19), we have the result: 



Ao 



E ^1." 



{3 3 
EE 
j=i 1=1 



k'l 



^1,1 ^l,r. 



1=1 



I I 

'^0,l'^0,j 







2 




-1 3 




+2i?'o • k'lr'^^i - 


R'o 


ki,i — 2ro_( 


E^M 

.1=1 


E 

i=l 


\k'i 



1=1 

3 

~ ^hi ^'l,r. 



(C20) 



Referring back to Figure 9, the term Aq is connected to R' f through the relation R' f = —XqR'q. Therefore, the criterion for 
waves to have multiple or single internal reflection is contained in the controlled quantity Aq. The vector R'q is a quantity that 
must be specified initiaUy. Because Aq is a positive definite scalar, we can rewrite it as 

( 3 ^-l/2 



Ao 



R' 



R'o 



^^'/ EK.] 



Substituting the above definition of Ao into equation (C20), the quantity U?'/ can be solved as 



R' 



R', 



.n=l 



{3 3 
EE 
3=1 1=1 



k'l 



^'1,1 ^'ly 



m=l 



0,1' 0,j 







2 




-1 3 




+2R'o ■ k'lr'o^i - 


R'o 


ki,i — 2rQ ; 


.1=1 


E 

i=l 


\\k'i 



1=1 

3 

~ ^l,r' 



m=l 



'0,i 



(C21) 
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Because the hemisphere opening has a radius r[, the following criteria are concluded: 

< r-, single internal reflection, 




(C22) 



where 



R' 



> r[, multiple internal reflections, 
is given in equation (C21) and r • is the radius of a hemisphere. 



3. Plate-Hemisphere 

The description of a surface is a study of the orientation of its local normal n'p, which is shown in Figure 10. In spherical 
coordinates, the unit vectors n'p, O'p and (j)'p are expressed as 



p — Z^i=i ^^p,i^%i Op — z^i^i Q0, ei, (p p — Z^i=i sing/ g^i ei, 



(C23) 



where 



(Op, ^p) = sin e'^ cos {O'p, </>;) = sin 6'^ sin A^g {9'^) = cos (C24) 

It is easy to show that the set of unit vectors (n'p, 6'p, $'p^ forms an orthonormal coordinates. Therefore, the points on plane 
can be described by a 2D coordinate system made of 6'p and ^'p, 

3 



'pA 



<6L dA: 



" d9'p sin O'p d(j)fp 



(C25) 



If the plane's orientation constantly changes in time about its origin, the points on the plane experience the velocity dR'p/dt, 



3 

E 

i=l 



(92 A' - 



sine;, \de'pd<t>'p 




sin6»: 



a^A' 



pd<t)'pde'p sine;, a [^/^ 



(C26) 



where 0'^, (j)'p are the angular frequencies and i^'p,0'^, f'p,<t>'^ are the lattice vibrations along the directions 9'p and (f)'p, respec- 
tively. Here, it is understood that A^ 3 is independent of cjj'p . Therefore, the differentiation of A^ 3 with respect to the cjj'p vanishes. 
For the static plate in which there are no lattice vibrations, y'pfi'^, and y'p,<j,i^ vanishes. 

For the case of plate-hemisphere configuration shown in Figure 11, the points on the plate are represented by the vector Rp 
relative to the system origin. Making the correspondence in equation (B2), R ^ Rn, Rt — > Rt,p and R' — * n'p, the two 
angular variable sets {O'p, (j)'p) and {6, <j)) are connected through the relations given in equations (BIO) and (B12) with — > 1. 

= 1. Therefore, we obtain 



Here r,' — > 1 because 



{O'p, 4>'p, yT,p,i,yT,p,2) = arctan 



yT,p,2 + sin6'pSm(f>'p 



9p (Op, (j)'p, Rt,p^ = arctan 



, yT,p,i + sm cos (pp ^ 



cos (arctan ( ^™ ) ) -|- sin (arctan ( '^^•''•^+'^'"^p 



sin (p'j, 
cos 0; 



(C27) 



(C28) 



where the subscript p of <f)p and 9p indicates that these are the spherical variables for the points on the plate shown in Figure 1 1, 
and they are not that of the hemisphere. The vector Rp becomes 

3 

Rp = Rt,p {yT,p,l, yT,p,2, yT,p,3) + R'p = E 



yT,p,i + Vpfii —^TTT. r 



p dO'p smO'p d(i)'p 



(C29) 
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The magnitude Rp is given by 



3 



+ 1^: 



dA' 



P'^p d0L sin 6*: 



1/2 



In terms of the spherical coordinates, Rp is expressed by 

^ ^ _ r 3 



i=l 



dA' 



sniO'p d4>'p 



1/2 



where 



Ap,i (^9p,(j>p^ = sin 0p cos (j)p, Ap,2 (^p,</'p) = sin sin (^p, Ap,3 (^^p^ = cos^p. 



(C30) 



(C31) 



Here, the subscript p in i?p indicates that the vector Rp describe the points on the plate. If the plane's orientation constantly 
changes in time about its origin, then the same orientation change observed relative to the system origin is given by the velocity 

dRp/dt, 



dA' 



-\ 2 



-1/2 



OA: 



p,k 



sin 9: 



p "^p 



3 3 

EE 

j=i k=i 

d'A'^,k 



" d0'„ 



sind'p \de'pd(j)'p 



- cot e. 




vfi' 



d'K.k 



d^A 



p.k 



p'pd0'p sinO'p d [(j)'^] 



f>'p + v 



dA'^.k ^ 5a;; 



dO'p 



sin 9'p dcp'p 



A. 



■P,3 



+E 



dA^ 

dO'p sin^;, d<p'p 



9Ap,j dOp ., 



dA 



p j W^Jp 



aft 



p ^-rp 



Pp d<P'p '. 



(C32) 



where it is understood that A^ 3 and Ap^ are independent of (p'p and ^p, respectively; and as a consequence, their differentiation 
with respect to (j>'p and (pp vanishes. Here 9'p, (f)'p are angular frequencies and i'T,p,i is the translation speed of plate relative to 
system origin. Also, i^'p.o'^, t^'p,4>'^ are lattice vibrations along directions 9'p and ^'p, respectively. For a static plate in which 
there are no lattice vibrations, i^'p.e'^ and i^'p.4,'^ vanishes. 

A cross-sectional view of the plate-hemisphere system is shown in Figure 12. The initial wave vector k'i traveling toward the 
hemisphere would go through reflections according to the law of reflection and finally exit. It then continues toward the plate 
and reflects from it. Depending on the orientation of plate at the time of impact, the wave would either escape to infinity or 
re-enter the hemisphere to repeat the process all over again. 

The equation (C25) defines points on a plate, as shown in Figure 1 1 of section (3.B.3), relative to the plate origin. If Sp is a 
set of points on a plate whose members are defined by R'p of equation (C25), the wave reflection dynamics off the plate involve 
only those points of Sp in the intersection between the plate and the plane of incidence whose unit normal is n'p^i given in 
equation (A17). In order to determine the intersection between the plate and the incidence plane, the plate is first represented by 
a scalar field. From equation (C23), the unit plate normal is 

3 

Tl'p — ^ ^ Ap^jCi- 



The scalar field corresponding to the unit normal n'p satisfies the relation 

3 ^ 3 3 



i=l 



V7pK,-^,-^)-E^^f7 = E^;,^^^ ^ E 



dv', f'* 



0. 
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The individual component of the equation is given by 

^ - A' - n ?• - 1 2 s 

where A^^ is independent of t*!. An integration with respect to u[ yields the result 



(C33) 



where the integration constant is set to zero because the plate contains its local origin. The intersection between the plane of 
incidence and the plate, shown in Figure 10 of section (3.B.3), satisfies the relation 



1=1 



'^'i = 0, 



(C34) 



where {i/'^,!/!^, v'^) is given in equation (A43), and u'^ is a scalar corresponding to the basis e,, of course. We have, from 
equation (C25), 



^' = ^-'"--^ + ;^^' ^ = 1.2,3. 



Substituting v'- into equation (C34), v'^ g, is solved as 



(C35) 



^3 9K^ 
V ,1 Z^i=l dd>' 









n'p,i 



sin 6*' 3 aA;, 
2^1=1 



OS' 









n'p,i 



(C36) 



where the summation over indices j, k, m and n is implicit, also the quotient ;/ ^, / sin 9p has been moved out of the summation. 
The R'p given in equation (C25) is then rewritten as 



" sin ft' ^ 



p i=i 



dA' 







n'p.i 


-1 




2^1=1 86'^ 




n'p,i 


-1 


de'p 



ei. 



{Cil) 



Similarly, Rp given in equation (C30) is rewritten as 



^T,p,i + 










n'p,i 


-1 

^ijkkijTQi^ 






d(j)fp 3 9a;, 

\^ Z^Z=l 80'^ 






-1 

eimn^l,rn.^O,n 





2 ^ 



1/2 



(C38) 



where Ap^j is given in equation (C31). If Nh,max is the maximum count of reflections within the hemisphere before the wave 
escapes, the direction of the escaping wave, measured with respect to the system origin ^ = 0, is 



(C39) 



where /civ^^^^+i = ^N^^^a^+ik' Nn,^aa:+i- Similarly, by the correspondence i?/,,]Vh,„„,+i ^ Rh,i+i and i?/,,]v^^„^ -R?i,i+2 
in Figure 12, the direction of the escaping wave vector k is equivalently described by the relation 



C^JVfe ™„^+i = Rp- Rh,Nh,, 



(C40) 
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where C is an appropriate positive scale factor. Combining equations (C39) and (C40), Rp is solved as 

Rp = C-R/t.JVh.max + l + [1 - C] Rh,Nh,max- 



(C41) 



Because both i?,, and k 



-1 belong to a spanning set for the plane of incidence whose unit normal is n'p,i given in equation 



(A17), we observe that the following relationship 



r,+i = in'p^i 



(C42) 



hold, where 7 is a proportional constant. Substituting the explicit form for n'p^i from equation (A17) into equation (C42), it 
simphfies into the following equation 

3 r ^ -i 

C^ijkRh,Nh,max + '^,jkNh,^^^ + l,k + [1 - C] ^i0kRh,Nh,max,3^Nh,max + '^,k + 7 '^'p,l ejjfc^lj^O.fe = 0, 

and its component equations are given by 

C^ijkRh,Nh,max+l,j^Nh,max+l,k + [1 ~ C] ^ijkRh,Nh,max ,j ^Nh,max+1-,k +7 

where kNh,rn.ax+i,k = Rh,Nh,max+i,k — Rh,Nh,m<ix,k as described in equation (C39). Finally, the scale factor ( is solved as 

-1 





-1 


n'p,i 


^ijkkijTQ^. 




-1 

^ijkkijTQ^. = 



eijkRh,Nh,max,jRh,Nh,max + '>;k " ^ijkRh,N h,max ,3 Rh,Nh,max,k +7 



^ijkkijrQ^. 



X [^ijkRh,Nh,max,jRh,Nh,max + i,k - ^ijkRh,Nh„,nax,jRh,Nh,,nax,k - ^ijkRh,Nh,max+l,jRh,Nh,max + '^,k 

+ ^ijkRh,Nh,max + i,jRh,Nh,^^^,k] ^ , (C43) 

where i, j, k = 1, 2, 3. Here, the notation d have been adopted in place of It should be understood that for irrotational 3D 
vectors, Ci = C2 = Cs = C- For vectors in 2D and ID space, it is understood then ^3, ^2 are absent, respectively. In current form, 
equation (C43) is incomplete because 7 is still arbitrary. This happens because f^^, and p'p^^, of -Rp, equation (C30), still needs 
to be related to the scale parameter Substituting Q for ( in equation (C41), it is rewritten as 

Rp = CiRh,Nh,m.ax + '^ + [1 - Ci] Flh,Nh,m.ax 

or using equation (C6) to exphcitly substitute for Rh,Nh,max+'^ ^h,Nh,niax ^ — ^h,max + 1, N = Nh,max, respectively; 
and, regrouping the terms 



R 



1/2 



1/2 



A 



^h,Nh,max + hi + [1 - Ci] 
\ 

i 

) 



(C44) 



where Ci = C2 = Cs = C- The subscript i of r ■ is not a summation index. Equating the above result for .Rp with that of equation 
(C38), we arrive at 

/ r r- , _ . r .... 1 -1 . . 2^ 1/2 



E 




VT,p,j + -7- 





/ ^a;^ 




n'p,i 


-1 

^xyzk-Y^y^Q^z 










n'p,i 


-1 


dO'p 

J 



1/2 



l2 



Nh,max + i,j 



. j = l 



^h,Nh,max + '^,i + [Ci - 1] Ah,Nh,., 

1/2N 



79 



and its component equations are 



3 

E 







A' + 




-1 

p V r' 




sin 6'p 








-1 





1/2 



1/2 



•'^'l.JVfe.rnax + l.i + [Ci " 1] ^h,Nh,max,i 

1/2 

2 I 



0, 



where i = 1,2, 3. Introducing the following definitions for convenience, 



1/2 



= E =1 



1/2 



Cc 



v^3 9a;_^ 

Z^a:=l d<bL 



A' + 



X / A' + 

X ^Z^(=i I^^S.i + 

the relation shown in equation (C45) is rewritten as 

3 



n'p,i 



n'p,i 
-1 



f y r' 



E 



/9A' - M' ■ 



sm.6'p \ d(j)'p 



de' 



-\ 2 



1/2 



Ap,i — Ci^cAft,JV;,,„„^ + l,i 

+ - 1] Bf Ah,JV;.,_,Z = 0, 



(C45) 



(C46) 



where i = 1,2,3. There are three such relations, one for each value of i. It is convenient to combine additively all three relations 
to form 



'J Pi3 



1/2 



E ^p>» ~ ^^^i E Ah,iVh,„„,+i,i 

i=l i=l 

3 

+ [Ci-i]BcE^^.^''— * = o- 



i=l 



After regrouping the terms and squaring both sides, it becomes 



E 

j=l L 



OA Z^i A/,..v,..,„„,,+i,, - [0 - 1] A,,..v,. 
Ef=i Ap,i 



The summation labeled is rewritten as 



. , sin^ 0' 
j=i p 



dA 



The above equation is simpUfied into a quadratic equation of v' ^, , 



sin 6*^ \ d(j>'^ 



'^T,p,j 



. , sin^ 6' 
j=i p 



aA', aA',^' 

P-,3 Pi3 

9<P'p ^ dO'p 



+ E 

i=i 



sinf 



aA' 



p,j \ , . ^ 



E;=i Ap,; 



= 0. 
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The two roots v'^^^, are given by 
/ 

-E 



3 r 



V 



smf 



OA' , dA' , 



89: 



p 




3 r 



sin 6*' 



5A' , aA' , 



. , sin^ (9' 



dA' , 9A' 



-I 2 



+ M 



1/2' 



3 I 

^ sin^ 61; 



9A 



pj 



d0'p 



(C47) 



where A.^, i?;^ and are defined in equation (C46). It is understood that one does not mix summation indices of Aq, 
and with those aheady present above. The result for i/ ^, is still incomplete because the factor 7 in Qi needs to be fixed 
by normaUzation. Unfortunately, the translation property of the plate, VT,p,j , niakes it difficult to extract Ci out of the radical. 
Besides the stated difficulty regarding Ci, v'^ ^1 is still ambiguous in deciding which of the two roots correspond to the actual 
reflection point on the plate. Fortunately, for the plate-hemisphere system of Figure 1 1, the choice of system origin is arbitrary. 
One can always choose the plate origin to be the system origin and the translation of the plate can be equivalently simulated by 
a translation of the hemisphere origin in the opposite direction. Then, in the rest frame of the plate, the translational motion of 
the plate is zero, i.e., i'T,p,j = 0. In this frame, i/^^, takes on much simphfied form 



^P,^; =±sin 



V O^C E?=i A/.,Ar,,_+i,i - id - 1] Be Ell H 



Ei=i 



■9A1. 



OS', , 



1/2 



Ei=i ^p,i 



For the sign ambiguity in f' j,, , it can be quickly fixed by noting that for i'T,p,j = 0, equation (C45) yields 



Ci Ej=l 



ArA 



Br A 



[sin^^]-ME,Lri^-^C^l ELA^ 



1/2 



P,3 



(C48) 



where A^, Bi; and are defined in equation (C46) with i'T,p,j = 0. It is to be noticed that for a situation where i'T,p,j = 0, A 
becomes identical to A' in form. One can obtain A simply by replacing the primed variables with the unprimed ones in A'. For 
convenience, Q of equation (C43) is rewritten as 



Ci = C-U^+jC-^B^, 



where 



Furthermore introducing the definitions. 



Af; 



y3 



the u' J, of equation (C48) is rewritten as 



■9A' 



1/2 



Ell A: 



P,U 



(C49) 



(C50) 



Cp^C-^A^Ap + jCg'C-'B^Ap + Cg'B^B, 



/3 "-^7 -°7^/3- 



sin^p, 



'^T,p,j 



0. 
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Substituting the last expression for f^^, into equation (C38), we arrive at 



Rp — 



dk' , OA' , 



1/2 



^Ap,iei, (C51) 



where i'T,p,i = 0. The vector cross product Rp x k^^^^^^+i is given by 



Rp X fcjv^ „„^+i — ^ijkRp,jkN^^^^^+i,kei 



or 



3 

^ ^ ^ijkA-pjkNh „^^^+l,k^i- 



9A' , 9a; , 



-I 2 



1/2 



Finally, substituting above vector cross product Rp x ^at;. „„^+i into equation (C42) and regrouping the terms, it becomes 
^ I [c^'C-'A^Ap + jC^'C-'B^Ap + C^'B^Bp] <j ^ 



i=l 



3 rfi\' . OA' 

p,i _ (J p,] 

b'p « dO'p . 



><<^ijkApjkNh,ma^+i,k - 7 '^ijkkijrQ^k ) = 



where n'p^i = n'p^i X^i=i ^ijkk'ijrQ ^.Ci have been used. And for the component equations 



C^'C:;'A,Ap + ^C^^C-^B,Ap + C^'BcBp] I ^ 



^ 'dA'^ OA' 



1/2 



^^ijkK,o^Nh,^a.+hk - 7 '^'p.i eijkkijr'o^k = 0, 
where i = 1, 2, 3. There are three such relations and they are additively combined to yield 

1/2 



Cp'C-'A,Ap + ^fC^'C-'B^Ap + Cp'B^Bp] I ^ 



dA' , OA' , 

d<i>'p dO'p 

.. 1 3 



-I 2 



X ^ ^ijkApjkNh,^a^+i,k - 7 '^'p.i X] ^»i'=^i,i^o,fe = 



i=l 



i=l 



Finally, 7 is solved to give the result 



7 = 7o 



I E '^ij'^KAk - Cp'C-'B^Ap I J2 
j=i [j=i 

E *'y'=-'^p.j''^-'V'.,'"ox+i.fe^ ( Cp^C'^A^Ap + Cp^Bf^Bfj 

3 

^ijkApjkNh,,na^ + i,k 



dA' , 5A' 



1/2 



p,j P'l 



x<E 



, 1/2 
21' 3 



(C52) 
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The parameter z/^^, is now completely defined, 



(C53) 



where {A^,Bt^,C^) , {A^, B^, C\) , (A^, B^, C^) and 7^ are given by equations (C46), (C49), (C50) and (C52), respectively. 
With fp^, defined in equation (C53), the reflection point on the plate is obtained from equation (C38), 



4= E 



Y-3 9a;,. 






-1 










-1 





1/2 



(C54) 



where i^T,p,j = 0. It should be noticed that for a situation where i'T,p,j = 0, A becomes identical to A' in form, and A can be 
obtained simply by replacing the primed variables with the unprimed ones. 

To see if the wave reflected from the plate at location Rp re-enters the hemisphere cavity or escape to inflnity, we consider the 
reflected wave fcr,Arh,max+i! 



(C55) 



where the relation found in equation (A14) have been used. As always, it is convenient to express vectors in component forms. 
Making the changes in variables {I, m, n) k) , (m, q, r) {j, I, m) , n' ^ n'p and k'i — > kN^ ,„^^+i, the component 

result of equation (A 16) is used to get 



3 3 



kr,N^,ma. + l = EE {"^'-L ['^Nh,ma. + hiK,kK,k " K,i'^^h,ma. + hkn'p,k] 

1=1 k=l 

-<^r,\\n'p^kkNH,,nax+hkn'p,i} ei, 



(C56) 



where n^^ and are coefficients of the normalized n'p. All wave vectors entering the hemisphere cavity satisfy the relation 



ei -^0 = 0, ^0 = E ['^'^'^,1 + r'o,i] ei, 



(C57) 



i=l 



i=l 



where ^„ is a real- valued positive scale vector and Ro is the points on the opening face of hemisphere. The scale vector ^„ has 
the form 



i=l 



With the scale vector ^„ defined above; and, Rp and kr,Nh,max+'^ defined in equations (C54) and (C56), respectively, equation 
(C57) is rewritten in component form 



E 

i=l 



E 



p.s 



Y-3 fw;.. 




n'p,i 


-1 

^ijkkijrQ). 






K,i + 


n'p,i 


-1 

^Imnki^^rQ,^ 


d0'p 



2 ^ 



1/2 



^C^ ^AjA/s + joC^ ^C^ ^B^Ap + ^B^B/s Ap^i + ^ {ar,± [kNh,max+^,i 



k=l 



X'np,k'^p,k - n'p,^kNn,r^^^+i,knp^k] - ar.lp^.fc^Afh.^a.+i.fe^p.i} - '^T,h,i - ro,t) Ci = 0, 
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which yields the component equations, 

3 



V3 






-1 






K,i + 


n'p,i 


-1 


dO'p 



1/2 



fc=i 



-n'p,ikNu,„,ax+hkn'p,k] - <^r,\\n'p,kkNh,„^ax+hknp^i} - VT,h,i - r'o^i = 0, 
where i = 1,2, 3. Finally, ^^^j is solved as 

/ r . _ , r 1 n . 2 ^ 1/2 



E 



V 



dA' ^'=1 




n'p^i 


-1 

^ijkkijTQ,^ 




9<^; 3 ^a;, 
2^1=1 ae'^ 


K,i + 


n'p,i 


-1 


de'p 



(C58) 



where i = 1, 2, 3. 

The above result can be applied in setting the re-entry criteria. Notice that < I'i^ which implies Tq j < r^, where is 
the radius of the hemisphere. It can be concluded then that all waves re-entering hemisphere cavity would satisfy the condition 
= ^K,2 = ^k;,3- On the other hand, those waves that escapes to infinity cannot have all three ^,^,1 equaling to a same 
constant. The re-entry condition j = 2 = Ck,3 is just another way of stating the existence of parametric Une along the 
vector kr.Nh ^ +1 that happens to pierce through a hemisphere opening. When such a line does not exist, the initial wave vector 
direction has to be rotated accordingly to a new direction, such that in its rotated direction there is a parametric line that pierces 
through the hemisphere opening; it leads to the condition that aU three £^^,1 cannot equal to a same constant. The re-entry criteria 
are now summarized for bookkeeping purpose. 



(C59) 



= ^K,2 = ^k;,3 : wave reenters hemisphere, 

ELSE : wave escapes to infinity, 
where ELSE is the case where ^^^i = ^^,2 = Ck,3 cannot be satisfied. 

Appendix D: DYNAMICAL CASIMIR FORCE 



The original derivations and developments pertaining to the dynamical Casimir force are included in this appendix. It is 
referenced by the text of this paper to supply all the fine details. 



1. Formalism of Zero-Point Energy and its Force 

For massless fields, the energy-momentum relation is given by 

"^ns = ^Total = PC, 



(Dl) 



where p is the momentum and c the speed of light. The field propagating in an arbitrary direction has a momentum p' 



J2^=i P'i^i- '^^^ associated field energy-momentum relation is hence 



1/2 



0. 
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The differentiation of the above equation gives 



1/2' 



1/2 



. i=l 



The total differential energy dH'^^ is 



z— 1 * ^ z — 1 



where the relation p'^ = [n^ + ^] fiA;^ has been used. The total differential momentum is 



1/2 



3 ^-1/2 3 

2 



The combined result is 



3 

E 

j=i 



i=l 



1 



i=l 



dW. 



dk' 



-1/2 



dp', = 0. 



Because all the momentum differentials are linearly independent, their coefficients are zero, 

r 3 ^-1/2 



1 



dW 



dk 



, i=l 



There are three such equations. Then, additively combining the three relations, and rearranging the terms, we have 



E 

.1=1 



1/2 3 

> E 



Squaring both sides to get rid of the radical leads to 

E' ' 

.1=1 



dk' 



dk' 



E^^ 



i=l 



i=l 



E^^ 



^2 2 ^3 "1 ^ 

The summations X^^^i \p'i\ and are rewritten as 

E b^]' = ip'af + E (1 - b:]' = ip'af + E (1 - 

i=l --I 





1— ' 




( 


n« + - 













2 |- 


j=i . 





?'a + E(l-'^-)?'^ 
1=1 

Kf +2^(l-5,„) 



3 



1 



1=1 



1 



E(l-^ia)p^ 

.j=l 

■ 3 

E(i-'^-)^^ 



(D2) 



(D3) 



(D4) 



(D5) 



where p'^ has been replaced by [n^ + /lA:-. Substituting the result into equation (D5) and rearranging the terms in powers of 
p'^, we have 



E(i-M 



/icl A;' 



i=l 
3 

E (1 - 



he] k'iP 



'iPa 



1 



h] [k'f=0. 
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Defining the following quantities, 



(D6) 



the above quadratic equation is rewritten as 



1 



he 



Sa,2P;-C2 2 + C^_,C„,3 = 0. 



Finally, the root p'^ is found to be 



Pa = 



C'.i -{[ns + \] he) \cl^ _ + 1] - ([n. + i] he) 



1/2 



(D7) 



where the positive root have been chosen since p'^ is the magnitude of the ath component of the total momentum p', therefore it 
is a positive scalar, > 0. 

By definition, the force is equal to the rate of change of momentum. 



a=l 



dt ^"^^'^ dt 



a=l 



a=l 



The explicit expression for J^'a is found to be 



Cli-{[ns + h]hc) 
\ 



Cli - (h. + 5] Ml -{[ns + ^] he)' 



2 [n, + i] ?lC„,iC„,2 



+ 



+ M' 



dCa,2 



dCa,l 

dt 



I 



+ 



[rig + ^] h?Ca,2Ca,A 



+ 



a, 3 



(D8) 



where 



a, 4 



C2 z'' 
Q.2 '~^a,l'~^a, 



- ( h« + 1] Ml + 5] M^ 



-1/2 



(D9) 



Before computing the three time derivatives dCa,i/dt, dCoc,2/dt and dCa^s/dt, we notice that k'^ (rij) = Uj/j (ij) . Hence, the 
derivative dk[/dt can be written as 



dk[ ^ dk^drn^ +n ——=f- (L ) —^h- + n —L 
dt dm dt ' ' '9Li dt ' * Srij ' ^ dLi * 



(DIO) 



The three derivatives dC^ .i/dt, dCa .2/dt and dC^ , 3 / are given by 



dk' dk' ~dt ~ ^ ^ f)y.riy. 

2—1 J — 1 J * 



3 3 C12-1 



dt 



z— 1 J — 1 J '' 



dk 



3 3 



i=i j=i 



(Dll) 
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dC, 



Q,2 



dt 



3 

E 



(1 - Sia) 



orii oLi 



(D12) 



a,3 



=2x;(i-M 



1=1 





1' 


2 










^2 





Ji [Li) T^Ui + Hi— —Li 



(D13) 



where Cq,_i, Cq.2 ™d Cq,,3 are defined in equation (D6). It is noted that the derivative dk'Jdt, and also each of dCa,i/dt, 
dCa,2/dt and dCa,3/dt, consists of two contributing parts, one is proportional to rii and the other involves Lj. The force 
expression in equation (D8) has then two contributing parts. The force contribution involving has a physical meaning that 
the boundaries are being driven to generate the extra wave modes that would otherwise be missing when such drivers were 
not present. The force contribution involving Lj is the effect of feedbacks from the moving boundaries. This feedback effect 
due to the moving boundaries tends to either cool or heat the conducting boundaries. For an isolated, non-driven conducting 
boundaries, the force contribution proportional to hi vanishes. The expression of force is then rewritten as 



2 [ris + i] ?iCa,lC„,2 



[rig + i] ^Ca,2Ca,i 



1 



1^2 

2 



2^(l-5ia) 



It can be simplified with the following definitions. 



Ca,5 — 



Cl^ - {[ns + I] hcY 

Ca.lCa.^CctA 



2 [n, + i] ftC„,iC„,2 



Lin ' 



Ca.2CaA 



[Us + i] ft 



+ 5] ^ '] - ( ["^ + 5] Cl, - ( [n, + i] ftc) 

C'a^lCc(,4 



2 ' 



C'c«,7 = 



c^-(K + i]ftc)^ 



The dynamical force can then be rewritten as 



I ^ f 



3 



dWi. 







1' 






1' 






Us - 


^2 




Us - 


'2 





' dLj dk'.dK 



(D14) 



(D15) 



(D16) 



(D17) 
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where Cq.5, Ca.e and Ca.7 are defined in equations (D14), (D15) and (D16). The force equation (D17) vanishes for the ID case, 
which is an expected result. The reason is explained as follow: Recall that equation (D4) reads 



E 



dn' 



dk' 




-1/2 



dp[ = 0. 



For the ID case, the summation runs only once and the above expression simpUfies to 



Us + 



dk' 



dp'i = 





1' 




( 











dk!: 



0. 



This is a classic situation where the problem has been over specified. For the 3D case, equation (D4) is really a combination 
of two constraints, bil^ ^^'^ '^'n,- '^^e ID case, there is only one constraint, H'^^. Hence, equation (D4) becomes 

an over specification. In order to avoid the problem caused by over specifications in this formulation, the one dimensional 
force expression can be obtained directly by differentiating equation (Dl) instead of using the above formulation for the three 
dimensional case. We have then for the force expression in ID case: 



^ c dt 
For an isolated, non-driven systems. 



c 



dk' 



dk' dt 



ldn'„ 
c 



dk' 



ndfdn'^ . 

T' = — Le 

cdL dk' ' 



(D18) 



where is the force expression in ID space. Here the subscript a of J^'a have been dropped for simphcity, since it is a one 
dimensional force. 



2. Equations of Motion for tlie Driven Parallel Plates 



Consider the one dimensional system of two parallel plates shown in Figure 14. Defining the boundary length as the 



magnitude of a vector 



Lsn ■ Si 



, where 3? denotes the region, the following relation is found from Figure 14, 

3 

i=l 



(D19) 



Hence, the velocity dL^/dt is 



dLsji _ dRrp^m. dRip^sR 
dt dt dt 

and the corresponding component magnitude is given by 



E 



dR. 



dR, 



dt 



dt 



e,; 



(D20) 



dR, 



dt 



■ei = 



7! * 



dt 



dt 



(D21) 



Substituting the result Lj^jj of equation (D21) for L„ in the one dimensional dynamical force expression of equation (D18), 



na,?R dfa^sft dH'^^ sj 



c dLa,?tt dk'„ 



'a,SR 



dR, 



dR, 



lp,St 



dt 



dt 



(mi) 



where Li.sR = and i = a. The subscript 3? denotes the corresponding quantities associated with the region 3i = 1,2, 3, e.g., 
jj denotes the field energy in region 3?. For simphcity, the following notational convention is adopted 

• _ dRa ^ ^ _ d'^Ra 
~ ^ ■ ~ ~d^ ■ 
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The force expression of equation (D22) is then rewritten as 



(D23) 



(D24) 



Before writing down equations of motion for each plates illustrated in Figure 14, the associated center of mass point relative the 
the surface point vectors -Rrp,SR,a for each plates needs to be determined. The center of mass point Rrp,cmfor plate labeled "right 
plate" in Figure 14 is related to the surface point vector Rrp,R through a relation 

Rrp,cm (t) = Rrp,U=2 {t) + Rrp,cm-n {t) , 

where Rrp,cm-^ {t) = Rrp,cm-2 {t) is a displacement between surface and the center of mass point. The ath component of the 
center of mass point Rrp,cm is then 

Rrp,cm,a (^) — ' Rrp,cm (0 — * Rrp,2 (^) ~t~ ' Rrp,cm — 2 (^) — Rrp,2,a (0 ~^ Rrp,cm — 2,a (^) • (-^^5) 

The component of the center of mass point speed is given by 

Rrp^cm,a (0 — Rrp,2,a (^) ~t~ Rrp,cm — 2.a (^) ■ (D26) 

Similarly, for the plate labeled "left plate," the center of mass point is related to the surface vector point Rip,^=2 {t) by 

Rlp,cm {t) = Rlp,U=2 {t) — Rlp,cm-^ (t) , 

and the component along the direction is 

Rip, (t) = Rlp,2,a (t) — Rip, cm— 2, a (i) , Rip (t) = Rlp,2,cx (t) — Rlp^ 

cm—2,a 

it) . (D27) 

Using the above center of mass relations, equations (D25), (D26) and (D27), along with the force equation (D24), the net force 
acting on a plate labeled "right plate" along the direction of in the configuration shown in Figure 14 is 



or 



Tn^rpRrp,cm,a — 9ci,2 



TTlrp Rrp , cm 

Rrp,t 



Rlp,< 



R 



rp,cm — 2,a 



Rip, 



cm— 2, a 



+ 9a,Z 



R. 



■dpr,cm,a 



■rp,cm^a. 



Rdi 



lpr,c'm—2,a 



R 



Tp,cm—2,(x 



(D28) 



where rrirp is the mass of the "right plate." If the plate surface is not vibrating longitudinally along the direction of e'^, the 
displacements Rrp,cm-2,a and Rdpr,cm-2,a are constants; hence, Rrp,cm-2,a = Rdpr,cm-2,a = 0. For static surfaces, the 
above net force relation simphfies to 



1^rpRrp,cm,a — ^rp,29oi,2 Rrp,cm,a Rlp,cm,' 



^rp,39oi,3 



R. 



■dpr,cm,oc 



(D29) 



where Srp.2 and Srp.s have been inserted for convenience due to the force sign convention to be set later. Similarly, for the plate 
labeled "left plate" in Figure 14, the net force relation along the direction of e'a is 

'mipRlp,cm,a {t) = ^'a,/3,K=l + J^' a, 0,^=2 

or, for the case where plate surfaces do not have longitudinal vibrations. 



mipRip 

.cm. a — Slp,l9a,l \Rlp 



R 



dpi, cm, a 



Slp,29a,2 



Rr 



R 



Ip.cm.a 



(D30) 



where mip is a mass of "left plate" and the terms sip^i and sip^2 have been inserted for convenience due to the force sign 
convention to be set later. We have now the two coupled differential equations. 



1^rpRrp,cm,a ~t~ Srp,39cK,3Rrp,c 



^rp,29a,2Rrp,cm,a ~t~ ^rp,29a,2Rlp,cm,a — ^rp,39ct,3Rdpr,cm,aj 
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mipRip Rip ,cm,a Slp,29a,2-^rp,cm,a — -Slp,iga,lRdpl ,cm,c 

Introducing the following definitions, 

' rji = [Srp,29a,2 - Srp.S^a.s] , % = -Srp,2ga,2'^rp 

'^3 = ^ip [Slp,l9a,l - Slp,29a,2] , '74 = Sip,25c«,2W,~\ 

Sffa.am^p Rdpr 6p = -sip,iga,irnip , 

-^1 ~ Rrp^cm,(xt R2 — Rlp^cm,(xi 

the coupled differential equations are rewritten as 

Ri - rjiRi - 772 7?2 = £.rp, R2 - - mRi = iip- 



(D31) 



(D32) 



The equations of motion shown in equation (D32) are a system of two linear second-order inhomogeneous differential equations. 
In order to rewrite the coupled linear inhomogeneous differential equation (D32) into a set of first-order linear inhomogeneous 
equation, a set of new variables are defined first. 



Ri = R3, R2 = Ri, 

R3 = Rl ^ Crp + rjiRi + 772i?2 = irp + ??1-R3 + ^72-^4, 

^ Ri = R2 = ^ip + ry3i?2 + V^Ri = 6p + V3R4 + mRa- 



(D33) 



Using these new variables defined in equation (D33), equation (D32) can be cast into first-order inhomogeneous equation in 
matrix form, 



R=M-R + ^ 



'Ri' 




R2 




R3 




_ R4 





10 

1 

771 ?72 

774 ?73 





[Ri] 




- - 




R2 









R3 


+ 






. Ri 




- iip - 



The above first-order inhomogeneous equation is equivalent to 



-Ri = / R^dt , R2 = / Ridt 

J to J to 



and 



i?3 








R3 


+ 


^rp 




















Mr, 


Rrj 





(D34) 



(D35) 



For the homogeneous system 



' Rs' 




' Vi 


V2 ' 




'R3' 


Ri 




Vi 






Ri 



(D36) 



R„ Mn Rr, 

the eigenvalues are found from the root of the characteristic equation 

det (a7- = - [771 773] A -I- 7?i773 - 772774 = 0. 

The two eigenvalues are 



^ m + % , / 1 r i2 , 

A3 = — ^ — + 1 4 [^1 ~ + ^2»74 



>. 1/2 
I ' 



A4 



[m - m] + mm 



1/2 



(D37) 
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And, the two corresponding eigenvectors are found to be 



1 



Ra = 



As - rii 



-1/2 



and 



i?A4 = -R3 



" 1 


, ^;3 = |l + 


'A4 - rji 




L ?72 J 




. V2 . 





-1/2 



where ^3 and ^4 are the normaUzation constants. The solutions for the matrix equation (D36) are then 



Ai^exp(A30 
exp (Aat) 



= i?A3 exp (Aat) = Ri 
The fundamental matrix solution $ (f) = (/>a3 (t) , 4>\i (t) is given by 



exp (A4t) = -R3 



exp (A4/:) 



- exp (A4i) 



^^i?4exp(A3t) i?3exp(A4t) 



i?4exp(A30 ^^ii;3exp(A4t) 



The fundamental matrix solution $ (t) has an inverse 



where 



det($(t)) 
det (i)) 



^i?3exp(A4i) -ii3exp(A4t) 
-hi exp (A3i) A^^^4 exp (Ast) 



A4 - rii 



1 



^3^4exp([A3 + A4] i). 



.A3 -r?i 

The principal matrix solution ^ {t, to) = $ (t) • (to) of equation (D35) becomes then 

1 



where 



*(i,io) 



V'li {t, to) = RzRa 



det 



(*(to)) 



V'li {t,h) i'12 {t,to) 

1p21 {t, to) V'22 {t, to) 



A4 - m 
A3 - m 



exp (As* + A4to) - exp {Xit + Xato) 



(D38) 



(D39) 



(D40) 



(D41) 



(D42) 



(D43) 



V'12 it,to) = R3R4 



V'21 {t,to) = R3R4 



V2 112 

exp {X4t + Xsto) - T exp {Xst + X^to) 



A3 -m 



A3 -m 



A4 — Vl / , , N A4 — TJi , , 

exp (Xst + Xito) exp (A4i + Xsto) 



fp22 {t,to) = R3R4 



A4 - T?l 



A3 -m 



exp {Xit + A3<o) - exp {X^t + A4<o) 



The inverse of principal matrix solution ^ {t, to) is 

^~'{t,to) = 



det (t,to)) det (to)) 



V'22 {t, to) -V'12 (t, io) 
-V'2i(i,io) V'ii(*>*o) 



(D44) 



(D45) 



(D46) 



(D47) 
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where 



det 



{t,to) 



dot 4> (to 



[V'll (*, *o) V'22 {t, to) - V'12 {t, to) V'21 (*, to)] ■ 



(D48) 



Using a variation-of-parameters technique, the solution to the inhomogeneous first-order differential equation (D35) is 

Rr, (t) = * {t, to) ■ Rr, (to) + ^ {t, to) ■ I to) ■ |* {t') dt' , 

Jto 

where it is understood the multiphcations are that of the matrix operations. Substituting into this integral equation the results for 
(i) ' ^ (*') ' * (*) and given by equations (D35), (D42) and (D47), 



Rs (t) 

RA{t) 



1 



det 



(^(io)) 





■ i?3 {to) ' 






R4 {to) 


+ 



to I det 



ipii (t, io) i/'i2 {t,to) 

V'21 {t,to) i>22 {t,to) 



^22{t\to) -^l2{t',to) 

-V'2i(i',io) V'ii(i',io) 



■0ii(i,*o) 1pl2{t,to) 

i^2l{t,to) ll'22{t,to) 



6p (i') 



dt' 



or 



det 



{to) 



V'll (t, to) {to) + 1pl2 {t, to) R4 {to) + 



' ^^22{t'M)irp{t') _ 



*o det 



{^{t'M)) 



^12 y,to)ap {t') ^^, 
to det (^*(t',to)) 



det (5 (to)) 
■012 (i,to) 



+ 



det 



^11 {t',to)^lp{t') 



to det 



dt' - /* 

Jto 



V'21 (i',to)^rp(0 



det 



-dt' 



(D49) 



R4{t) 



det 



($(io)) 



V'21 (i, to) R3 {to) + ^22 {t, to) R4 {to) + 



V'22 {f, to) ^rp {t'\^^, _ /■* V'12 {t'M)^lp {t') ^^, 



to det 



to det 



{^{t',to)) 



V'21 *o) 
det (to)) 

^ jA22 (i,to) 

det 



(*(io)) 



/ 



* V'll {f, to) ^Ip (f) /■* V2I (i', to) ^rp (*') 



/ 



'to det (f, to)) "^to det (f, to)) 

It is noted from equation (D34), Rs {to) and R4 (to) are initial speeds, 



-dt' 



{to) = Rl {to) = R3 {to) , Rlp,cm,a {to) = R2 {to) = -^4 {to) ■ 



Hence, 



Rrp,cm^a {t) — 



A4 (;to) - rii (;to) 



A3 (;<o) - Vi {;to) 



V'll (t,to) Rrp ,cm,Q 

(to) + V'12 {t,to) Rip {to) 

exp ([A3 (;to) + A4 (; to)] to) 



+ ^n{t,to) [ 

Jto V'll {f, 
Jto ^11 {f 



V'22 {f, to) ^rp {t') - V12 {f, to) ilp (f) 



, to) V'22 {f, to) - V12 (*', to) V'21 (*', to) 



dt' -h V12 {t, to) 



V'll to) 6p (*') - V'21 (f, to) ^rp (f) 



to) V'22 (f, to) - V12 (f, to) V2I (f, to) 



dt', 



(D50) 



(D51) 
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-^lp,cm,a (^) — 



A4 (;to) - m (;to) 

h {■,to) - m (;io) 



- 1 



V'21 {t, to) Rrp, 



rp.cm.a 



i'22 {t,to) Rip (to) 



exp([A3(;to) + A4(;io)]io) 



+ V'21 

rt 



^22 {t',to) irp {t') - V'12 (i', <0) W 



X 



/o (*'' *o) ^22 (i', io) - V'12 (i', to) ^21 (t', io) 

Vll (f, to) ^22 (f, to) - V12 (f, to) V2I (f, to) ' 



rft' + V22 (t,to) 



(D52) 



where substitutions have been made for the determinants det (to)^ and det {t', to)^ from equations (D41) and (D48). It 

is to be understood that the notation (; to) on r^i , A3 and A4 impHes impHcit time dependence for these terms. Finally, integration 
of both sides of equations (D5 1) and (D52) with respect to time gives the results 



R. 



■rp,cm^(x 



it) 



A4 (;to) - r/i (:to) 

A3 (;to) - vi (;to) 



(to) +V'12 (t, to) Rip 

,cm,a (to) 

exp([A3 (;to) + A4(;to)] to) 



Jtt 



i^22 (t',to)?rp (t')-V'l2 (t',to)6p (f) 



to) V22 (f, to) - V'12 (f, to) V2I (f, to) 



rft' + Vl2 {T,to) 



Vll (t^to)6p (tQ- V2I {t',to)^rp jt') 

V-ii (f, to) V22 (f, to) - V12 (f, to) V21 (f, to) 



dt' 



dr + R 



Tp,cm,a 



(to): 



(D53) 



RLp,cm,a (t) 



A4 (;to) - ?7i (;to) 

A3 (;to) - m (;to) 



V'21 (t, to) Rrp (to) +V'22 (r, to) Rip (to) 

exp([A3(;io)+A4(;to)]to) 



+ V'2i (T,to) / , 

Vii(t', 



V22 (f, to) Crp jt') - V12 jt', to) ^ip jt') 



, to) V22 (f, to) - V12 (f, to) V21 (f, to) 



rft' + V22 (T,to) 



V'll jt', to) ^ip jt') - V21 (t^ to) irp jt') 

V-ii (f, to) V22 (f, to) - V12 (f, to) V21 (f, to) 



dr + Rip (to), 



(D54) 



where the terms V'll, V'12, V'21 and V'22 are defined in equations (D43), (D44), (D45) and (D46). The remaining integrations are 
sttaightforward; hence, their explicit forms are not shown. 

As a closing remark of this section, one may argue that for the static case, Rrp,cm,a {to) and Rip^cm,a (to) niust be zero 
because the conductors seem to be fixed in position. This argument is flawed for any wall totally fixed in position upon impact 
would require an infinite amount of energy. One has to consider the conservation of momentum simultaneously. The wall has 
to have moved by the amount /SRwaii = Rwaii^t, where At is the total duration of impact, and Rwaii is calculated from the 
momentum conservation and it is non-zero. The same argument can be applied to the apparatus shown in Figure 14. For that 
system 



Rrp, cm, a (to) = \\Rlp,3 (to) + Rrp,2 (to) 



'^virtual— photonW — ^n.,SR (to) , 



Rip (to) = -Rrp,l (to) + Rlp,2 {to) 



or, for simplicity, assuming an impact along the normal direction, 

2 



■^rp,cm,a (^o) — 



TThfpC 



\'^'n^,3 (to) — '^ns,2 (to)|| , Rlp,cm,a (to) 



mipC 



|'^n,,l (to) - T^'ns,2 (to) I 



where the difference under the magnitude symbol imphes that the energies from different regions act to counteract each other. 
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